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Preface 


Equally  as  important  as  balloon  thermodynamic  problems  (and  perhaps  more 
challenging)  are  balloon  aerodynamic  drag  problems,  and  the  interrelated  prob¬ 
lems  of  dynamically  determined  balloon  bubble  shapes.  A  review  of  the  literatui  e 
reveals  that  far  from  being  solved,  these  latter  problems  have  barely  been  de¬ 
fined.  Although  we  can  draw  no  lasting  comfort  from  Icnowing  that  our  status  is 
not  unique,  #e  should  certainly  be  encouraged  by  the  fact  that  others  are  vigor¬ 
ously  (and  with  some  success)  pursuing  solutions  to  very  similar  problems.  A 
quite  concise  statement  of  our  mutual  difficulty  is  the  following;^ 

The  most  basic  problem  of  determining  the  equilibrium  figure  of  (the 
body)  also  requires  a  simultaneou.s  evaluation  of  the  flow  around  (the  body), 
which  in  turn  depends  on  the  shape.  It  is  doubtful  that  (the  body)  attains  a 
true  equilibrium  figure  under  natural  conditions,  and  the  analytical  prob¬ 
lems  arising  from  the  coupling  of  the  flow  and  shape  are  unsolved,  even 
for  the  steady  case .  * 

Quite  simply  stated,  this  paper  re-evaluates  the  efforts  that  have  gone  into 
proaucing  practical  aerodynamic -thermodynamic,  flight  performance  models.  It 
has  been  justified  primarily  by  those  programs  that  require  accurate  ascent  rate 
prediction  or  the  ability  to  account  for  the  effects  of  balloon  motion,  insofar  as 


In  the  quotation  we  have  substituted  the  term  "the  body"  wherever  Green  used 
the  term  the  raindrop. 

1.  Green,  A.  W.  (1975)  An  approximation  for  the  shapes  of  large  raindrops,  J. 
Appl.  Metcorol.  14:1578-1583. 
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they  affect  data  obtained  from  balloon-borne  sensors.  Although  our  findings  are 
significant,  they  are  not  conclusive.  However,  it  is  hoped  that  the  results  will  be 
useful  from  both  theoretical  and  practical  perspectives.  To  the  extent  that  tliey 
are,  much  of  the  credit  must  go  to  Mrs.  Catherine  Hice,  who  is  both  the  work  unit 
and  task  scientist,  without  whose  cons  tructive  criticism  the  arguments  herein 
would  not  have  been  as  well  organized  or  (we  trust)  as  convincing. 
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Factors  Affecting  the  Verticol  Motion  of  a 
Zero-Pressore,  Polyethylene,  Free  Balloon 


l.  INTRODUCTION 
1.1  Oliij(ic(iv«i> 

Knowledge  of  the  vertical  motions  of  f»'ee  balloon  systems  is  important  for 
four  reasons,  h  irst,  to  effect  a  desired  ascent  profile,  we  i>\ed  to  understand 
how  a  balloon  system  responds,  both  to  its  free*Uft  force  and  to  citanges  in  its 
operational  environment.  Second,  to  initiate  ascent  (or  dc  jccnt)  or  to  vary  the 
rate,  we  mvtst  be  able  to  predict  the  effects  of  both  deballasting  and  gas  valving. 
Thirdly,  to  design  a  balloon*bornc  experiment,  we  often  need  to  be  able  to  pi'e« 
diet  balloon  system  motions.  Fcurtli,  to  rednee  cortoin  flight  sensor  data,  we 
must  be  able  to  iiccount  for  effects  of  system  motions  on  scientific  ok-  c-  ations. 
With  these  needs  clearly  eslablishcd,  we  ic-prcscnt  ihe  relevant  .j,  and 
restate  the  problems  with  historical  comments  and  fresh  ir.sights.  In  this  way, 
we  .seek  to  bring  about  better  vinderstunding  of  the  problems,  greater  oa.se  in 
modifying  them,  and  improved  agreement  between  observed  and  predicted  system 
performance. 


(lleoeived  for  Publication  30  May  10S5) 
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1.2  liHckfruiutd 

I 

Tho  first  step  toward  achlcvinj;  our  objoctivos  was  to  review  the  rall\er  esc- 
tensive  literiuure  in  two  of  throe  probtetn  areas  Utat  deal  with  balloon  ascent  and 
float  motions.  These  fii'St  two  ureas  and  the  iiiore  significant  works  reported 
Uiercunder  are; 

(a)  Ascent  Hate  Prediction  Models 
1041  Clarke  and  Korff^ 

•> 

1082  University  of  Minnesota 
195!)  Ney  and  Winkler** 

1980  Erickson  and  Froohlich'’ 

197-1  Nelson®  _ 

1976  Kremser^ 

(b)  Comprehensive  Flight  Performance  Models 
1949  Smith  and  Murray® 

1952  Hall® 

1952  University  of  Minnesota*® 

1 06 1  Einslie** 

1063  Dingwell  ot  al*^ 

1966  Garmelos*® 

1970  Hansen*** 

1973  Fujii  etal*® 

1974  Kreith  and  Kreidcr*® 

1070  Romero  ot  al*^ 

1076  Dalis  Crema*® 

1001  Carlson  and  Horn*® 

1001  Carlson  and  Horn”® 

The  third  problem  area,  Float  Altitude  Motions,  encompasses  some  of  tl»e 
principal  ooncernr-  that  compelled  this  .study.  Becau.so  of  the  scale  of  these  mo¬ 
tions,  and  the  strictly  operational  nature  of  most  free  balloon  efforts,  it  is  not 
surprising  that  literature  in  this  ar  a  is  scarce.  However,  there  is  extensive 
reporting  on  related  natural  inotioim  of  the  atmosphere,  and  on  the  motions  of 
floating  su|)erprcssur  ?  balloons.  A  chronological  sampling  of  .such  works  in¬ 
cludes  : 


Because  of  tho  large  number  of  references  cited  above,  they  will  not  be  listed 
here.  See  Heferences,  page  41. 


(c)  Flout  Mtiluck*  Motions 

IP'U*  Smith  untl  Mvu'rwy*^ 

2 1 

1P50  FnunOiis  el  ftl 

22 

1P52  Uiuvoi'sily  of  Minnosoia 
lOfifi  lUrsch  and  Booker'^ ^ 

24 

1908  Nishlinura  and  Mlrosawa 

IPB9  Morris  and  Stofan^^ 

1P71  Nlshlmwva  ct  01“° 

*>7 

1974  Lcvanon" 

in7R  Levaiwn  and  Kushnir^® 

15)77  Julian  ot  nl^^ 

1978  Masstnan^® 

iVo  must  keep  in  nilnd  that  In  the  late  l940s  there  was  no  one  who  had  any 
experience  with  large  balloons  made  fi’om  Inextensible  plastic  fLlnis.  In  truth, 
balloons  having  volumes  in  excess  of  10  million  ft'  did  not  become  common  until 
15)55),  and  early  balloon  envelopes  were  made  from  cither  relatively  incxtoislble 
coated  fabrics  or  the  very  oxti?nslblo  radiosonde  balloon  material,  We  must  also 
boar  In  mind  that  a  large-scale  digital  computer  was  seldom  available  in  the  early 
days  of  the  plastic  balloon;  tills  both  prevented  and  discouraged  ntlompts  to  solve 
or  even  to  define  rigorously  many  of  the  problems. 


21.  Emmons,  G, ,  et  ttl  (1950)  O.ieillations  in  the  strniosphere  and  high  tropo* 

sphere.  Bull.  Am.  Motoorol.  Goc.  3UMo.  4)!l35“130. 

22,  l'mvpr.miy  of  Minnesota  (195:'.)  )-*roRress  Iteport  on  Iteseni’cji  iind  Develop- 

inoni  itt  tlie  Field  of  High  .‘XKItudc  Halloons.  Volume  iX,  Contract  Nonr- 

tToToit;; 

2:1.  liirsch,  .).  11..  and  Booker,  1)  B.  (}i)fiC)  Response  of  suporpressuro  b.al* 
loons  to  verlival  air  motions.  J.  Appl.  Metcorol.  5:226-229. 

21.  Nishimura,  J. .  and  Uirosa\v.n,  H.  (1968)  The  hunting  mcJianism  of  plastic 
balloons,  ISA.S  null.  4(1B):93-110. 

25.  Morris,  A.  1...  and  Stefan,  K.  11.  (1969)  IUrIi  Altitude  Balloons  as  Scientific 
Platforms.  National  Center  foi  Atmospheric  Be.search. 


26.  Nishimura.  J.,  et  al  (1!)7  1)  Balloon  behavior  during  level  flight,  IS.AS  Bull. 

7(10:357-268. 

27.  Levanon,  N. ,  ct  al  (1974)  On  tlte  behavior  of  superpressure  balloons  at 

150  niB,  J.  Appl.  Metcorol.  j^3:4‘J4-501. 

28.  LA>vanon,  N. ,  and  Ktishnir,  Y.  (1976)  On  the  response  of  superpressure  bal¬ 

loons  to  displacements  from  equihbrium  dcttsit.v  level,  J.  Appl.  Metcorol. 
J^5:346-:)49. 

29.  Julian.  P. ,  et  al  (1977)  The  TWERU=.  experiment,  Bull.  Am.  Mcteorol.  Soc. 

5n(No.  n):n3G-9l8, 

30.  Massman,  W.  J.  (15)78)  On  the  natvmc  .of  vertical  oscillations  of  constant 

volume  balloons,  J.  Appl.  Metcorol.  _lV:l3.5l»13r)C. 
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1.2.1  ASCENT  RATE  PREDICTION  MODELS 

Early  balloon  users  sought  (and  we  still  seek)  mathematieally  simple  models 
for  predicting  balloon  ascent  rate.  Even  now  we  anticipate  that  a  comprehensive 
aerodynamic-thermodynamic  flight  performance  model  will  enable  us  to  perform 
the  types  of  factor  sensitivity  analyses  necessary  for  the  development  of  such  a 
rate  prediction  model.  However,  i'ith  the  increased  capabilities  of  new  mici’o- 
computer.s  we  may  soon  have  both  the  capacity  and  the  speed  needed  to  solve  the  * 

comprehensive  flight  performance  model  (both  interactively  and  in  real-time)  on  ■ 

location  at  remote  balloon  launch  sites. 

A 

The  University  of  Minnesota^  developed  the  most  frequently  used  and  studied 
ascent  rate  prediction  model.  Gildenberg*  used  this  model  to  analyze  many  bal¬ 
loon  ascents;  his  objective  was  to  improve  the  accuracy  and  applicability  of  the  i 

model  by  refining  its  thermodynamic  and  aerodynamic  coefficients.  Nelson's 

efforts  were  directed  along  the  same  lines.  He  published  tlie  results  of  these  ; 

efforts  to  use  the  model  for;  (a)  balloons  with  volumes  from  1  to  30  million  ft  , 

(b)  balloons  carrying  payloads  weighing  between  100  and  10, 000  lbs,  and  (c)  bal¬ 
loons  ascending  before  and  after  sunset.  Like  Gildenberg,  he  had  limited  success. 

(ionaidering  what  might  be  learned  from  using  a  complete  flight  performance 
model  (such  as  proposed  herein),  it  is  possible  that  the  University  of  Minnesota's 
model  might  be  enhanced  to  serve  as  a  practical  ascent  rate  predictor.  Thus,  we  ■ 

reproduce  the  model  so  that  the  reader  might  appreciate  its  relative  simplicity 

(compared  to  the  proposed  comprehensive  flight  performance  model).  We  note  j 

that  the  originator.s  urged  caution  in  its  use,  t 

F  Cl  ’I'  (G  *  (  L  «  V  +  4  *  v|**3/l  T;»:'21)':'«(0.  25)  ’! 

+  C:2  ‘I’  (v!''!-2)  *  (P  «  G/T)>:“:=(l/3)  j 

where  ] 

F  is  the  free  lilt  force,  normalized  by  dividing  by  the  weight  of  the  dis-  j 

placed  air,  ; 

Cl  is  the  thermodynamic  coefficient,  7.4E-04,  j 

.i 

■  J 

’'B.  D.  Gildenberg  reti.’ed  as  a  meteorologist  and  balloon  operations  controller  at 
the  AFOL's  Holloman  AFB  balloon  facility  in  New  Mexico.  His  work  on  this 
problem  is  contained  in  unpublished  notes  and  letters  to  his  coworkers. 

t  Arithmetic  -iperations  in  this  report  are  expres.sed  in  FORTRAN  operator  sym-  ■ 
bols.  'niese  symbols  are  widely  recognized,  but  we  include  their  definitions 
here  to  ensure  against  misinterpretation:  addition  (+),  subtraction  (-),  multi¬ 
plication  (*),  division  (/),  and  exponentiation  (**). 
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G  ii  the  weight  of  the  disiplaccd  air,  lbs, 

L  is  the  atmospheric  lapse  rate,  degrees  C  per  1000  ft, 

V  is  the  ascent  rate,  ft  per  min. 

T  is  the  air  temperature,  degrees  K, 

C2  is  the  aerodynamic  coefficien'  ti.  5E-07, 

P  is  the  atmospheric  pressure,  mB. 

7 

Kremser's  ascent  rate  prediction  model  seems  to  be  the  only  model  that  is 
at  least  as  sophisticated  as  that  of  the  University  of  Minnesota,  but  it  has  not 
been  extensively  evaluated.  Therefore,  we  believe  that  a  continued  effort  with 
this  type  of  model  would  be  totally  unproductive  at  this  time. 

1.2.2  COMPREHENSIVE  FLIGHT  PERFORMANCE  MODELS 

As  with  the  ascent  prediction  models,  we  find  that  there  are  two  main  prob¬ 
lem  areas,  namely,  aerodynamic  drag  and  thermodynamic  drag,  now  designated 
as  heat  transfer.  Smith  and  Murray^  treated  both,  but  provided  limited  details 
on  flight  thermodynamics.  Their  primary  thermodynamic  concerns  were  atmos¬ 
pheric  temperature  lapse  rate  and  the  balloon  gas  superheat. 

They  treated  aerodynamic  drag  in  the  conventional  manner;  they  used  a 
spherical  shape  -  the  only  comparable  shape  for  which  drag  data  was  readily 
available.  On  the  other  hand,  they  did  note  thsi,  during  the  early  part  of  the  as¬ 
cent,  balloon  shape  is  characterized  by  a  "flabby,  unfilled  portion",  which  would 
uffect  the  aerodynamic  drag.  Although  this  was  a  significant  observation,  they 
failed  to  capitalize  on  it  because,  perhaps,  balloon  float  altitudes  at  that  time 
were  relatively  low.  Thus,  the  flabby  portion  would  exist  for  an  insignificant 
time  during  the  ascent. 

Smith  and  Murray  made  a  further  observation.  They  noted  the  inability  to 
solve  the  general  equations  of  motion  using  mettiods  then  available. 

g 

Hall  recognized  aerodynamic  and  thermodynamic  drag  problems  as  two  of 
the  pi-incipal  difficulties  of  his  day: 

The  greatest  uncertainties  in  the  analytical  premises  are  in  the  magni¬ 
tude  and  variation  in  the  drag  coefficient,  . . .  and  the  mechanism  and 

rate  of  heat  transfer  between  the  balloon  and  its  environment. 

Apparently  for  the  same  reasons  as  Smith  and  Murray,  he  used  the  sphere  as 
his  model  for  balloon  shape.  However,  he  noted  that  the  balloon's  "flexibility" 
would  "certainly"  be  a  governing  factor  with  respect  to  the  drag  coefficient. 

Hall  appears  to  have  been  the  first  to  include  in  a  flight  performance  model 
the  heat  transfer  processes;  (a)  free  and  forced  convection  between  air  and  bal¬ 
loon  film,  (b)  free  convection  between  inflation  gas  and  balloon  film,  (c)  solar 
energy  input,  and  (d)  infrared  heat  exchange  between  balloon  film  and  environment. 


His  work  was  organized  for  analysis  rather  than  for  solution  by  numerical  meth- 
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ods;  his  style  is  that  of  University  of  Minnesota  researchers. 

One  cannot  overemphasize  the  efforts  and  contributions  of  the  University  of 
Minnesota  research  team  to  both  the  aerodynamics  and  thermodynamics  of  free- 
flight  plastic  balloons.  Their  work,  which  also  included  design  of  both  balloons 

and  flight  instrumentation,  is  recorded  in  16  volumes  (with  flight  data).  In  these 

3  10 

volumes  there  are  four  specific  comments  '  on  aerodynamic  drag  problems  -  ■ 

that  are  most  relevant  to  our  effort  to  provide  an  improved,  reformatted,  com¬ 
prehensive,  flight  performance  model;  ' 

Unfortunately  in  the  case  of  an  ascending  nonextensible  type  balloon,  the  ^ 

shape  is  not  constant  with  altitude;  varying  from  the  shape  of  a  small  " 

sphere  with  long  depending  folds  of  fabric  at  take-off  to  roughly  a  spher-  \ 

ical  shape  at  altitude.  It  is  not  possible  therefore  to  use  a  single  func¬ 
tion  CD(Rn)  to  predict  the  drag  at  all  altitudes.  ] 

The  difficulty  in  predicting  the  value  of  CD  in  advance  lies  in  the  fact  J 

that  the  shape  of  the  balloon  is  not  constant  with  altitude,  and  thus  one 
cannot  carry  over  the  results  of  wind-tunnel  experiments  on  any  par¬ 
ticular  shape  of  model. 

The  described  dynamic  pressure  loading  on  the  balloon  is  such  as  to 
make  the  balloo'i  more  oblate  than  its  original  natural  shape. 

It  is  not  difficuii  to  imagine  that  a  balloon  free  to  change  shape  with 
dynamic  forces  will  hav  lone  so  appreciably  before  the  velocity  of 
dimpling  has  been  l  ea c  id.* 

Again,  the  difficulty  in  solving  the  drag  problem  see>ns  to  have  caused  investiga-  j 

tors  to  ignore  it  at  least  to  the  extent  that  they  did  not  attempt  to  synthesize  | 

their  observations  mco  a  formal  statement  .  he  problem,  t  j 

EmsUe'  ^  was  the  first  person  to  p  i  an  aerodynamic -thermodynamic  j 

flight  performance  model  in  a  system  ations  for  solution  on  a  large-  le 

digital  computer.  Tliis  system  of  ecv  s-  s  was  his  means  of  investigating  ...al- 
loon  dynamics.  It  included;  (a’’  perfci  oas  law,  (b)  an  equation  of  vertical 

motion,  (c)  a  gas  energy  equation  and  f'l;  ,  equation  for  fabric  or  film  energy  j 

(film  is  now  the  preferred  no  .  j 

Emslie  made  two  e  ‘  •.  iy  nt  points  that  relate  to  our  current  ef- 

fort.  First,  although  he  u.sea  a  w.iswint  drag  coefficient  (again  a  tecnnological 
expedient),  he  noted  the  gross  asymmetries  in  the  folds  of  the  film  below  the  gas 


-  i 

^Dimpling  velocity  is  defined  as  that  velocity  at  which  the  dynamic  pressure  in¬ 
duced  by  the  ascent  rate  equals  the  internal  pressure  on  the  crown  of  the  balloon. 

t". . .  a»  ther  characteristic  which  unifies  science  is  the  ability  to  ignore  prob¬ 
lems  h  are  not  yet  capable  of  solution  . . 

31.  Rivi  t,  P.  J.  (1983)  A  world  in  which  nothing  evei‘ happens  tv/ice,  J.  ?er. 

Hsch.  Soc.  34(No,  8):6B1. 
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bubble  when  the  balloon  is  below  its  float  altitude  (Figure  1).  Second,  he  ex¬ 
pressed  the  mass  of  the  film  surrounding  the  gas  bubble  as  a  function  of  the  en¬ 
closed  volume.  Although  he  understated  the  mass,  basing  it  on  the  surface  area 
of  a  sphere  witii  a  volume  equal  to  the  volume  of  the  enclosed  gas  bubble,  he 
nevertheless  recognized  that  the  mass  of  the  film  iiivolved  in  the  heat  transfer 
process  was  not  constant  (see  Figure  2). 


Figure  1.  Representative  Configuration  of  an  Ascending 
Tailored  Natural  Shape  Balloon,  2.01  Million  ft^,  Flight 
No.  H8 1-012.  Note  the  extreme  asymmeti’y  of  both  the 
gas  bubble  and  folds  of  undeployed  film  below  the  bubble 
shortly  after  launch 

Prior  to  Emslie's  work,  performance  models  were  used  for  relatively  low 
altitude  balloon  flights  and  relatively  small  balloons  (with  a  few  exceptions,  of 
course).  His  work  came  at  a  time  when  we  were  beginning  to  fly  routinely  at 
altitudes  well  above  100, 000  ft  and  on  balloons  with  volumes  of  10, 000, 000  ft^  and 
larger.  For  such  high  performance  systems,  ascent  ballast  capacity  was  at  a 
premium  and  ascent  times  were  on  the  order  of  hours  —  especially  at  night. 

Thus,  models  developed  to  forecast  accurately  the  ascent  profile  required  even 
more  accurate  formulation.  Therefore,  it  is  unfortunate  that,  in  expanding 
Emslie's  work,  and  in  translating  it  into  computer  code,  Dingwell  et  al  restated 
the  film  energy  equation,  using  constant  film  mass.  Instead,  they  should  have 
redefined  Emslie's  model  of  the  relationship  between  shell  mass  and  instantaneous 
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Figure  2.  A  Large  Tailored  Natural  Shape  Balloon,  21.  77 
Million  ft^.  Note  that  the  balloon  envelope  material  Is  not 
concentrated  around  the  gas  bubble,  but  is  distributed  be-> 
tween  the  gas  bubble  and  the  "rope"  of  film  yet  to  be  de¬ 
ployed.  Balloon  is  shown  shortly  after  launch  on  Plight  No. 

H78-052 

volume,  expressing  the  fact  that  the  actual  area  of  the  balloon  surface  that  en¬ 
closes  the  gas  is  greater  than  the  surface  area  of  the  enclosing  shape,  but  that  It 
is  far  less  than  the  constant  area  of  the  whole  envelope  (see  Figure  3),  except  in 
the  vicinity  of  the  natural  float  altitude. 

In  continuing  Emslie's  work,  Dingwell  et  al  developed  a  system  of  nine  simul¬ 
taneous  equations  to  solve  for  the  following  dependent  variables  as  functions  of 
time: 

(1)  altitude 

(2)  vertical  velocity 

(3)  gas  temperature 

(4)  film  temperature 

(5)  gas  weight 

(G)  instantaneous  balloon  volume 

(7)  atmospheric  pressure 

(6)  atmospheric  density 

(0)  payload  weight. 
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FRACTIONAL  SURFACE  AREA 

Figure  3.  Approximate  Mass  of  the  Shell 
Surrounding  the  Gas  Bubble  as  a  Fraction 
of  Total  Film  Mass,  in  Belatlon  to  the 
Ratio  of  Instantaneous  Specific  Lift  to  Spe¬ 
cific  Lift  at  Natural  Ceiling  Altitude,  Ihe 
approxinnation  based  on  the  simplified  par¬ 
achute  shape  model  is  given  in  subroutine 
MYBLN.  The  spherical  approximation  is 
based  on  the  function; 

R  =  0.  5  *  {1  -  cos  I  r->  (V/Vniax)**(l/3)) }, 

where  R  is  the  fractional  surface  area, 

Vmax  is  the  maximum  volume,  and  V  is 
the  instantaneous  volume.  Note  that  for  a 
balloon  with  a  ceiling  altitude  of  100, 000 
ft,  the  fractional  surface  area  at  launch  is 
about  10  percent  of  the  entire  envelope 
area 

To  solve  these  equations,  one  had  to  determine  certain  parameters  as  functions  of 

time  and  the  above  dependent  variables.  These  were  the  relevant  heat  transfer 

coefficients,  heat  transfer  areas,  and  optical  properties.  For  these  they  provided 

tables,  graphs,  and  mathematical  models,  including  those  for  heat  transfer,  both 

by  forced  and  free  convection,  and  by  radiation. 
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Germeles  reported  extensions  and  improvements  of  analyses  and  computer 
codes  reported  by  Dingwell  et  al.  Furthermore,  he  continued  the  constant  film 
mass  error.  Hansen^^  made  use  of  this  work  by  Germeles.  but  apparently  did  not 
modify  it. 


1  ^ 

In  Japan  Pujii  et  al*^  developed  a  far  less  sophisticated  routine;  it  included  a 

constant  drag  coefficient  (0.  2)  and  a  he  .'ansfer  model  that  did  not  explicitly 

contain  the  film  mass  term. 

1  fi 

Kreith  and  Kreider  ”  further  refined  the  work  of  Germeles  and  Dingwell  et  al, 
but  made  no  reference  to  Emslie.  Although  they  added  a  routine  to  compute  CD  as 
a  function  of  Reynolds  number,  and  made  significant  changes  in  some  of  the  heat 
transfer  models,  they  left  the  computer  codes  substantially  the  same  as  those  re¬ 
ported  by  Germeles.  Their  work  is  now  the  generally  accepted  standard. 

17 

Romero  et  al  referred  to  the  works  of  both  Germeles,  and  Kreith  and 

18 

Kreider,  but  Balls  Crema  et  al  referred  to  only  Kreith  and  Kreider.  Neither 
work  cited  Emslie  and  both  continued  to  use  constant  film  mass,  Romero  et  al, 
however,  did  define  separate  values  for  drag  coefficients  in  turbulent  and  laminar 

flows,  0.45  and  1.35  respectively. 

19 

Carlson  and  Horn  followed  the  lead  of  Kreith  and  Kreider,  using  a  system 
of  eight  roughly  "equivalent"  equations  —  including  the  film  mass  term,  still  rep¬ 
resented  as  a  constant.  However,  they  modified  the  assumptions  to  allow  the  in- 
flatant  to  absorb  and  emit  energy.  This  is  a  significant  change  that  reqw'res 
further  study  before  (.eneral  acceptance,  especially  in  light  of  some  of  the  changes 
in  this  proposed  model. 

With  regard  to  balloon  shape,  Carlson  and  Horn  assumed  tiiat 

. . .  over  much  of  the  flight  profile  the  balloon  shape  is  close  to  a  sphere. . 

and,  accordingly,  also  used  a  drag  coefficient  model  based  on  Rejmolds  number 
only,  In  addition,  they  observed  that  computed 

.  .  balloon  ascent  velocities  between  launch  and  tlie  tropopause  are  very 
sensitive  to  the  values  of  CD, 

and  suggested  that  the  balloon  could 

. . .  experience  significant  skin  friction  drag  in  addition  to  the  pressure 
drag  normally  found  on  a  sphere. 

20 

Carlson  and  Horn  added  significantly  to  their  previous  observations  when 
they  commented  that  the  apparent  virtual  mass  coefficient,  used  in  tlie  first  equa¬ 
tion,  might  be  inappropriate  for  "the  balloon  configuration.  " 
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1.2.3  FLXDAT  ALTITUDE  MOTIONS 

The  vertical  motions  of  a  zero-pressure  balloon  at  float  altitude  are  complex. 
If  we  are  to  understand  them  at  all.  we  must  also  understand  the  zero-pressure 
balloon's  interaction  with  its  use  environment.  Although  we  found  that  very  little 
exists  on  this  subject,  we  did  find  the  following  to  be  both  pertinent  and  interest¬ 
ing. 
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(a)  Smith  and  Murray  noted  the  influence  of  vertical  winds  on  vertical  bal¬ 
loon  motions.  They  even  included  these  wind  models  in  their  flight  performance 
model,  but  the  effects  at  today's  relatively  "high"  altitudes  may  be  insignificant. 

(b)  Nishimura  and  Hirosawa^^  treat  a  "hunting"  motion  that  relates  to  the 

balloon  construction  (a  subject  area  not  to  be  ignored  with  regard  to  its  effects  on 

balloon  performance  at  natural  float  altitude). 

30 

(c)  Massman  comments  that  . .  the  Brunt-Valsala  oscillations  of  the  bal¬ 
loon's  EDS  (equilibrium  density  surface)  can  have  a  period  as  short  as  5  min  in 
the  stratosphere."* 

1.3  Concliuions  Ihicd  on  the  Literature  Search 

Based  on  findings  and  observations  that  appear  in  the  cited  literature,  we 
believe  that  four  areas  deserve  particular  attention  in  the  formulation  of  any  flight 
performance  model;  (a)  balloon  shape  during  ascent,  (b)  mschanism  for  heat 
transfer  between  the  gas  and  the  film,  (c)  aerodynamic  drag  coefticient;  and  (d) 
added  mass. 

1. 3.  1  BALLOON  SHAPE  DURING  ASCENT 

The  shape  taken  by  the  partially  full  balloon  (gas  bubble),  while  either  as¬ 
cending  or  floating,  governs  the  effective  envelope  mass,  the  gas  bubble  surface 
area,  and  the  areas  Involved  in  heat  transfer  processes  (Figure  4).  Also  during 
ascent  and  descent  the  shape  affects  the  drag  area  (Figure  5)  and  the  air  flow 
around  the  balloon,  hence,  the  aerodynamic  drag  coefficient.  Possibly,  as  we 
shall  see,  it  also  affects  tlte  added  mass. 

Clearly  ascending  (or  descending)  balloon  shapes  are  far  fiom  spherical  — 
even  though  the  leading  surface  of  the  gas  bubbles  in  Figures  6  and  7  appear  to  be 
hemispherical.  Factually,  a  partially  full  balloon  is  asymmetrical  in  every  plane, 
and  this  asymmetry  is  further  exaggerated  by  the  gore  deployment,  which  is 
governed  in  turn  by  the  gore  pattern,  t  Figures  I,  2,  and  4  show  bubble  shapes 
quite  typical  of  today's  large,  fully  tailored,  natural  shape  balloons;  the  maximum 

*The  Brunt-Vaisala  period  is  defined  to  be:  2  *  ff/SQRT  [g  (Beta  +  dT/dh)/Tl 
seconds,  with  the  terms  defined  In  Brunt. 

tl'or  a  short  commentary  on  the  development  of  balloon  gore  patterns  see 
Dwyer.  33 

32.  Brunt,  D.  (1927)  The  period  of  simple  vertical  oscillations  in  the  atmosphere. 

Quart.  J.  Roy.  Meteorol.  Soc.  53:30-32. 

33.  Dwyer,  J.  F.  (1978)  Zero  pressure  balloon  shapes,  past,  present,  and 

future.  Scientific  Ballooning  (COSPAR),  W.  Riedler,  Ed.,  Pergamon 
Press,  pp.  9-19. 


Figure  4.  Balloon  Ascent  Shape.  Flight  No. 

H84<-003,  'Hie  early  ascent  shape  and  relative 
gas  bubble  surface  area  of  this  tailored  natu^ 
ral  shape  balloon,  model  no.  SV  -0178,  differs 
considerably  from  those  of  the  fully  tailored 
natural  shape  balloons  shown  in  Figures  1  and 
2,  It  is  an  intermediate  size  heavyload  balloon 
having  a  maximum  volume  of  5.  142  million  ft^ 

horizontal  cross  section  of  each  of  these  balloons  is  far  less  circular  than  those  of 
earlier  balloons  made  with  either  rectangular  or  semi'tailored  gores  {compare  the 
bubble  shapes  in  Figures  1  and  7). 

Dut  iug  inflation,  the  balloons  in  Figures  6  and  7  assumed  the  characteristic 
shape  shown  in  Figure  8.  Gore  pattern  types  that  produce  such  shapes  ai'e  shown 
in  Figure  9  as  patterns,  numbers  2  and  3.  and  (to  a  lesser  degree)  pattern  num¬ 
ber  4.  On  the  other  hand.  Figures  10  through  12  are  representative  of  the  pre- 
launch  shapes  of  balloons,  such  as  those  shown  during  ascent  in  Figures  1,  2.  4. 
and  5.  These  latter  shapes  arc  characteristic  of  balloons  made  with  fully  tailored 
gores,  pattern  number  1,  Figure  9. 

In  addition  to  the  gas  bubble  proper,  we  should  also  consider  the  shape  and 
effects  of  the  trailing  undeployed  balloon  shell.  Tltis  is  the  mass  of  film  that  (as 
we  noted)  Smith  and  Murray  called  the  . .  flabby  unfilled  portion. . .  "  and  rt  - 
searchers  at  the  University  of  Minnesota  described  as  ". . .  long  depending  folds 
of  fabric. . "  Emslie  also  noted  the  gross  asymmetries  in  the  folds  of  undeployed 
material,  and  Carlson  and  Horn  called  attention  to  tlie  fact  that  large  balloons,  in 
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Figure  S.  Ascent  Cenfigur&tion  of  a  itelatively 
SmaH.  Fully  Tailored  Balloon.  Model  No.  LTV- 
019,  Having  a  Maxlmuin  Volume  of  628,000  ft3. 

Flight  No.  H81-006.  The  highly  assymetric 
horizontal  cross-section  is  common  in  a  balloon 
of  this  size  at  liftoff 

the  early  stages  of  ascent  (when  the  existence  of  this  surplus  is  most  evident), 

have  shapes  that  are  .  significantly  different  from  that  of  a  sphere. .  "  It  Is 

probable  tliat  this  shape  feature  plays  an  Importont  role  in  determining  the  drag 

34 

coefficient,  much  as  in  the  case  of  a  sphere  with  a  splitter  plate. 

To  better  represent  overall  balloon  shape  in  a  computable  configuration,  we 
selected  the  existing  balloon  shape  model  shown  in  Figure  13.  This  contrived 
shape  has  at  least  two  distinct  advantages.  First,  it  provides  smooth  transition 
from  the  modeled,  partially  full  state  to  the  Ml,  natural  shape  state,  a  most  im¬ 
portant  consideration  in  the  analysis  of  vertical  motions  that  occur  at  or  near  jdte 
natural  ceiling  altitude.  Second,  it  permits  reasonably  accurate  computation  of 
the  instantaneous  mass  of  the  balloon  film  involved  in  the  heat  transfer  process  — 
including  load  cap  film,  if  such  is  present. 


34.  Hoerncr,  S.  F,  (19G5)  Fluid-Dynamic  Drag  (published  by  author). 

35.  Dwyer,  J.  F.  (1980)  The  Problem ;  Instantaneously  Effecting  Controlled  Bal¬ 

loon-System  Descent  from  High  Altitude,  AFGL-TH-80-0277,  AD  A100255. 


Figure  6.  Ascent  Configuration  of  an 
Early  Moby  Dick  Balloon  Having  Semi- 
Tailored  or  Rectangular  Gores  and  a 
Maximum  Volume  Less  Ilian  100,000 
fH,  Flight  No.  S-140,  27  November*. 
10r»3.  llie  hemispherical  crown  .and 
large  air  pockel  are  characteristic  of 
balloons  constructed  with  the  afore¬ 
mentioned  gore  patterns 


1.3.2  MECHANISM  FOR  HEAT  TRANSFER  BETWEEN  THE 
GAS  AND  THE  FILM 

Tr.nditionally,  one  has  assumed  that  the  mechanism  for  heat  transfer  between 
the  gas  and  balloon  wall  is  free  convection.  Wo  have  no  direct  evidence  to  support 
this  assumption;  neither  do  we  Imve  knowledge  of  the  sensitivity  of  the  models 
based  on  this  assumed  mechanism.  However,  we  do  have  evidence  that  an  ascend¬ 
ing  balloon  is  quite  asymmetric  (see  Figures  1,  2,  5,  and  7),  and  upward  looking 
cameras  have  shown  how  balloons  rotate  considerably  during  ascent.  These  two 
facts  suggest  that  the  gas  should  be  in  constant  motion,  agitated  by  large  internal 
v.'ines  of  the  envelope  material  that  are  surplus  to  the  instantaneous  Shape.  Tliere- 
forc.  one  might  ask  whether  this  heat  transfer  is  due  to  forced,  rather  than  to  free 
convection.  The  differences  in  computed  flight  pcrfoi’mances  based  on  these  oppos¬ 
ing  assumptions  arc  not  known;  should  tliey  not  be  significant,  then  the  more  easily 
computed  model  should  be  used.  Wc  do  not  yet  have  sufficient  reason  to  break 
with  convention  on  this  issue  and  consequently  continue  to  use  free  convection. 


'-i-V*'. V.,\ 
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Figure  7.  Ascent  Coiifiguration  of  ft  Seini*TailoJ'ed  Bullootw 
Model  No.  TrV-OOlF,  Flight  No.  H73'Olfi.  This  tapcless, 
semi-cylinder  balloon,  having  a  volume  of  804,000  ft^,  lias 
tlie  same  configuration  characteristic  of  the  Moby  Dick  bal¬ 
loon  shown  in  Figure  6 

1.3.3  AERODYNAMIC  DRAG  COEFFICIENT 

Drag  coefficient  models  found  in  the  literature  have  r.'ingetl  from  a  single 
value,  through  a  five-part,  piecewise  continuous  function  of  Reynolds  number, 
to  the  categorical  conclusion  that  a  single  ftmetion  relatioiuship  between  drag  coef¬ 
ficient  and  Reynolds  number  is  not  possible.  In  the  literature,  three  bases  for 
arguments  support  the  latter  conclusion:  (a)  inconstant  shape,  (b)  sluipc  deform - 
ability,  and  (c)  dimensional  reasoning.  Based  on  the  arguments  that  follow,  we 
conclude  that  any  valid  model  tliat  is  to  determine  acrodjmamic  drag  coefficients 
must  consider  at  least  two  dimensionless  variables:  Reynolds  luimbor  and  Fronde 
number.  Further,  because  the  shape  has  no  single  characteristic  length,  it  is 
reasonable  to  expect  that  we  will  need  a  third  dimensionless  variable,  fractional 
volume.  Tills  latter  variable  is  defined  us  the  ratio  of  the  instantaneous  volume 
to  the  maximum  volume;  it  has  tlie  effect  of  normalizing  the  shape  .so  that  a  single 
dimension  of  the  balloon  might  serve  as  a  characteristic  length  to  determine  the 
effective  drag  area. 


IS 


Figin'«  3.  tun&tion  Cootiguration  of  a 
Scmi'Tallored  Balloon,  Model  No.  TTV“ 
OOIF,  Flight  No.  H81-014,  During  in¬ 
flation  the  tapeless,  semi-cylinder  bal¬ 
loon  is  charactcri7.ed  by  highly  irnegulat' 
dcployii'ient  of  inatorial  excess  to  Uie 
amount  needed  to  enclose  the  gas  bubble. 
What  appear  to  be  load  tapes  are  alumin¬ 
um-backed  polyethylene  seam  reinforce¬ 
ments;  tltcse  make  the  system  more 
radar  reflective 


1 . 3 . 3 . 1  Inconstant  Shape 

The  University  of  Minnesota  study  summarises  well  the  argument  based  on 
inconstant  shape: 

Unfortunately  in  the  case  of  .an  ascending  non-oxtensible  type  balloon,  the 
shape  is  not  constant  with  altitude,  varying  from  U»e  siwpc  of  a  small 
sphere  with  long  depending  folds  o'  fabric  at  »  ike-of.  to  xvdghly  a  .spher¬ 
ical  shape  at  altitude.  It  is  not  possible  thevefore  to  use  a  single  function 
CDiRn)  to  predict  the  drag  at  all  altitudes. 

( emphasis  added) 

Schlichting,  also  in  this  regard,  notes  that  the  use  of  Reynolds  nuiiiber 
alone  presupposes  both  the  same  sliapc  and  orientation.  When  one  considers  the 


36.  Schlichting,  H.  (i068)  Doundary-layor  Tlteory,  McOraw-Hill  Lo>k  Co. ,  New 
York.  6th  Edition,  p.  16. 


l^lgure  ii,  Gore  PoUern  Types  for  Natural  Shape  Balloons.  Tlie  patterns 
shown  are:  (1)  fully  tailored,  (2)  rectangular.  (3)  seini*tallored,  and 
(4)  tapered  tatjgent.  These  respectivo  patterns  are  used  to  construct  the 
following  balloon  types;  (a)  fully  tailored  balloons,  |l>)  cylinder  balloons, 

(c)  seini'cyllncior  (sometimes  called  tailored  tapelcss)  balloons,  and 

(d)  tailored  balloons 

docmnenteci  balloon  .shapes  and  the  differences  between  the  forward  surfaces  pre¬ 
sented  by  rising  and  descending  balloons,  it  is  clear  tliat  these  two  prcsvippositions 
aru  invalid.  Haiiuc,  use  of  Iteyuolds  number  alone  is  insufficient. 

1 . 3 .  3 . 2  Shape  Dof ormal)  ility 

We  have  recognized  for  a  long  time  that  there  is  a  significant  difference  be¬ 
tween  the  .static  shapes  of  ft'cc  balloons  just  prior  to  launch  and  the  djmanUc 
shapes  taken  by  the  same  balloons  during  ascent.  Wo  concluded  from  tins  that: 

(a)  the  ascent  sitapes  of  a  balloon  represented  defoi'mntions  of  static  shapes, 

(b)  the  aerodyuamic  $hnpe.s  of  the  balloons  were  dependent  on  norodynamio  drag, 

(c)  .sohitions  to  the  drag  problem  involved  free  surface  phenomena,  and  (cl)  the 
Froude  number  should  play  a  role  equally  ns  important  as  the  Reynolds  number. 
We  are  still  surprised  tliat  tlcc  connection  between  shape  deformation,  free  sur¬ 
face  phenomena,  Frovede  number,  and  aerodynamic  drag  was  not  made  long  ago, 
in  the  earlier  motion  studies. 

!n  the  litci'ature  on  balloons,  there  arc  numerous  obseiwations  that  suggest 
the  applienbilit  •  of  free  surface  phenomena  to  the  aerodynamic  drag  problem. 

Hall  implied  .such  in  his  use  of  the  term  "flexibility",  and  in  his  certainty  that 
flexibility  would  bo  a  governing  factor  in  the  determination  of  drag.  Tice  Univer¬ 
sity  of  Minnesota  study  likewise  implied  the  relevance  of  such  phenomena  when  it 
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Figure  10.  Inflation  Configuration  of  a 
Fully  Tailored  Balloon,  Model  No.  LTV- 
OIB,  Havir^  a  Maxinoum  Volume  of 
355.000  ft%  Flight  No.  H80-029.  From 
the  top  end -fitting  outward,  the  fully 
tailored  gores  are  uniformly  deployed 

described  the  dynamic  pressure  acting  downward  on  the  top  of  the  balloon,  making 
it  .  more  oblate  than  its  original  natural  shape.  "  It  (U.  of  M. )  further  rein¬ 
forced  our  unique  and  unusual  interpretation  by  the  comment,  "It  is  not  difficult  to 
imagine  that  a  balloon  free  to  change  shape  with  dynamic  forces  will  have  done  so 
appreciably  before  the  velocity  of  dimpling  has  been  reached.  "  [emphasis  added) 
Interpretation  of  this  deformation  as  a  free  surface  phenomenon  is  also  strength¬ 
ened  by  the  Carlson  and  Horn  reference  to  the  possibility  that  the  balloon  apex 
region,  during  the  early  stages  of  ascent,  might  be  . .  more  oblate  than  a  sphere 
due  to  pressure  differences  across  the  film. " 

Where  such  deformations  do  occur,  Schlichting  notes  that  drag  based  only  on 
Reynolds  number  is  invalid  and  that  the  F'rcude  number  must  be  considered.  In¬ 
deed,  when  we  re-examined  the  University  of  Minnesota's  treatment  of  dimpling 
velocity,  we  found  that  the  ratio  of  dynamic  pressure  tc  static  internal  pressure 
could  be  reduced  to  the  Froude  number. 

We  considered  three  other  aerodynamic  problem  areas  to  be  potentially  en¬ 
lightening  with  respect  to  free  surface  phenomena  and  Froude  number,  insofar 
as  they  relate  to  shape  dcformability:  they  were  raindrops,  parachutes,  and 


Figure  11.  Inflation  Configuration  of  a 
Fully  Tailored  Balloon.  Model  No.  LTV  - 
013A,  Having  a  Maximum  Volume  of 
2.  9  Million  ft3.  Flight  No.  H80-039. 
From  the  top  end -fitting  outward,  the 
fully  tailored  gores  are  uniformly  de¬ 
ployed 


air -supported  structures.  Although  none  of  these  areas  yielded  anything  directly 

applicable  to  our  problem,  the  review  did  provide  Some  rewarding  insights. 
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Gillaspy  commented  on  the  raindrop  problem,  one  which  is  quite  analogous 
to  ours; 

A  sphere  falling  in  a  fluid  medium  will  attain  a  constant  or  terminal  ve¬ 
locity.  When  falling  at  terminal  velocity,  all  of  the  forces  on  the  sphere 
are  in  equilibrium.  If  the  sphere  is  composed  of  solid  material,  this 
equilibrium  is  the  balance  between  the  weight  and  the  aerodynamic  drag 
forces  on  the  sphere.  However,  in  a  liquid  drop  the  balance  is  much 
more  complicated.  Other  forces  arise  from  the  fact  that  the  drop  is 
liquid  and  deformable.  ( emphasis  added  ] 


37.  Gillaspy,  P.  H.  (1981)  Experimental  Petermin^ion  of  the  Effect  of  Physical 
Properties  on  the  Drag  of  Liqul3~!Bropr..  l^h.l!).  Dissertation,  tlnxverslty  o 
Nevada  (funded  under  U.  S.  Army  Research  Office  Contract  No.  DAAB29- 
77-G-1072). 


Figure  12.  Inflation  Configuration  of  a 
Fully  Tailored  Balloon.  Model  No.  SV- 
017B,  Having  a  Maximum  Volume  of 
5. 142  Million  ft^.  Flight  No.  H84-003. 
From  the  top  end -fitting  outward,  the 
fully  tailored  gores  are  uniformly  de¬ 
ployed 


He  developed  a  model  that  accounted  for  influences  of  Reynolds  number  and  Bond 
number.  L.iUe  the  Froude  number  in  our  proposed  approach  to  balloon  aerody¬ 
namic  problems,  the  Bond  number  accounts  for  the  effects  of  gravity. 

Perhaps  balloon  problems  are  generically  closer  to  parachute  problems  than 
they  are  to  raindrop  problems,  primarily  because  the  .stresses  in  a  balloon  shell 
are  closer  to  stresses  in  a  parachute  canopy  than  they  are  to  surface  tensions  in 

a  raindrop.  On  this  premise.  Von  Karnian's  introduction  of  the  Froude  number 
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into  the  analysis  of  parachute  opening  shock  provides  some  additional  encour¬ 
agement  that  the  Froude  number  might  indeed  be  one  key  to  the  solution  of  aero¬ 
dynamic  drag  problems  involving  balloons. 


Figure  13.  Simplified  Parachute -Shape  Model, 
In  this  model,  V  is  the  volume  of  the  inflation 
gas,  A  is  the  apex,  Bfjg  is  the  base  of  the  nat¬ 
ural  shape  generator.  Bps  is  the  base  of  the 
model,  T  is  the  point  of  tangency  between  the 
lower  portion  of  the  model  and  generator  shape, 
and  the  actual  balloon  gorelength  is  the  distance 
between  A  and  Bpg 


We  found  in  the  recent  literature  on  air-supported  structures  only  limited 

references  to  Froude  number;  the.<ic  were  with  respect  to  large  tensioned,  pneu- 
39  40 

malic  structures.  '  We  cite  them  only  to  indicate  that  application  of  free  sur¬ 
face  phenomena  to  problems  dealing  with  deformable  barrier  surfaces  has  further, 
and  more  recent,  precedent. 

Clearly  balloon  deformability  justifies  our  interpretation  of  balloon  aerody¬ 
namic  drag  as  a  free  surface  problem.  Furthermore,  it  supports  the  conclusion 
that  we  cannot  determine  aerodynamie  drag  coefficients  by  Hcynolds  number  alone; 
that  we  must  also  consider  another  dimensionless  variable  —  the  Froude  number. 


39.  Tryggvason,  B,  V,,  and  Isyumov,  N,  (1978)  Similarity  requirements  for  in¬ 

flatable  structures.  Proceedings  of  the  'Fhird  U.  S.  National  Conference 
(on  Wind  Engineering  lle.seardb).  University  of  Florida,  Gainesville,  Floi'- 
ida,  pp.  335-330. 

40,  Tryggvd.son,  B.  V ,  (1979)  Aeroclastic  modelling  of  pneumatic  and  ten.sioned 

fabric  structures.  Proceedings  of  the  Fifth  International  Conference  (on 
Wind  Engineering),  Fort  Collins,  Colorado,  pp.  1081 -1072. 
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1.  3.  3.  3  Dimensional  Reasoning 
41 

Finally,  Landau,  on  the  basis  of  dimensional  reasoning  argues: 

If  the  force  of  gravity  has  an  important  effect  upon  the  flow,  then  the 
latter  (drag  force)  is  deternjined  not  by  three  but  by  four  parameters 
. .  (including  the  acceleration  of  gravity) —  From  these  parameters 
we  can  construct  not  one  but  two  independent  dimensionless  quantities. 

These  can  be,  for  instance,  the  Reynolds  number  and  the  Froude 
number. . .  (and)  . . .  two  flows  will  be  similar  only  if  both  these  num¬ 
bers  have  the  same  values. 

1.3.4  ADDED  MASS 

The  conditions  under  which  the  added  niass  term  applies  are  not  well  defined, 
20 

and  Carlson  and  Horn  questioned  (as  did  we)  whether  the  term  is  applied  prop¬ 
erly  to  the  configuration  (of  the  partially  full  balloon  in  vertical  motion).  We  as¬ 
sume  that  their  doubt  (like  ours)  applies  only  under  the  stated  conditions,  when 
that  film,  excess  to  the  instantaneous  bubble  shape,  can  deploy  asymmetrically 
to  form  pockets  of  ambient  air  (see  Figures  6  and  7).  This  type  of  deployment 
would  have  the  effect  of  increasing  the  volume  in  the  added  mass  product  repre¬ 
sented  as  CM^BUOY  [see  Eq.  (1)  in  Section  2.1],  This  seems  quite  probable  only 
at  low  and  intermediate  altitudes,  before  the  nadir  cone  angle  becomes  great 
enough  to  prevent  the  formation  of  such  air  pockets.  If  adjustment  of  the  added 
mass  term  is  ever  deemed  necessary,  the  added  mass  coefficient  CM  can  prob¬ 
ably  be  developed  as  a  function  of  fractional  volume  for  each  specific  balloon  con¬ 
struction  type.  This  is  not  thought  to  be  required  for  the  large,  fully  tailored 
balloons  that  predominate  today.  Consequently,  such  a  function  will  not  be  intro¬ 
duced  at  this  stage  in  the  development  of  a  comprehensive  flight  performance 
model. 


2.  A  NEW  FLIGHT  PERFORMANCE  MODEL 

In  the  following  model,  tiie  operator  D(  )  represents  the  first  derivative  with 
respect  to  time.  The  definitions  of  terms  in  the  equations  are  given  once  in  the 
subsection  on  each  equation,  and  again  in  the  glossary,  where  (along  with  their 
dimensions)  their  FORTRAN  names  are  given.  Complex  model  components  are 
explained  separately;  otherwise  they  are  clarified  by  explicit  comments  included 
in  the  program  FORTRAN  codes  or  implicitly  by  references  to  specific  equations 
and/or  figures  appearing  in  cited  documents. 


41.  Landau,  L.  D. ,  and  Lifshitz,  £.  M.  (1959)  Fluid  Mechanics,  Addison-Wesley 
Publishing  Co. ,  Inc. ,  p.  63. 
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2.1  Hie  Differential  Equations 

The  following  eight  differential  equations  define  the  model  when  the  balloon  is 
partially  full;  either  when  it  is  floating  or  when  it  is  moving  vertically,  upward  or 
downward.  In  these  cases:  (a)  the  gas  volume  V  is  free  to  expand  or  contract, 

<b)  tlie  gas  pressure  Pg  is  a.ssumed  to  follow  the  ambient  pressure,  (c)  the  tem¬ 
perature  of  the  gas  Tg  reacts  to  the  gas  expansion  or  contraction,  and  (d)  the  gas 


weight  Wg  is  constant  if  the  apex  valve  is  closed. 

D[v)  =  (BUOY  -  DRAG  -  WS)/l  (WS  +  CM  *  BUOY)/Gl  (1) 

D[WJ  =  -  DB  (2) 

D(  Z\  =  V  (3) 

DI  Tf|  =  (Q1  +  Q2  +  Q3  +  Q4  -  Q5  +  Q6  +  Q7)/(CF  WE)  (4) 

D1  Wg|  =  -  VV  «  SWG  (5) 

D[  Tg)  =  (-  01  -  07  +  Q8  +  Qf)  -  QlO  +  on  -  SW  V  «  v)/[  Wg  *  (CV  +  RG)] 

(6) 

DIVJ  =  V  >:=  {D(Wg|/Wg  +  DlTgl/Tg +v/(RA  *  Ta))  (7) 

Dl  Pg|  =  -  SW  V  (8) 


When,  however,  the  volume  of  the  gas  bubble  equals  the  maximum  volume  of 
the  balloon,  and  the  gas  in  the  balloon  is  still  expanding,  either  gas  must  be  ex¬ 
pelled  through  the  ducts  or  the  balloon  will  eventually  burst.  To  model  this  vent¬ 
ing,  new  relationships  are  required;  Eqs.  (5a)  through  (8a)  [which  replace  Eqs. 

(5)  through  (8)1  provide  just  such  a  model.  In  execution:  (a)  the  gas  pressure 
changes  as  the  balloon  vents  gas,  (b>  the  bubble  volume  remains  constant,  (c)  the 
gas  temperature  reacts  to  restricted  expansion,  and  (d)  the  weight  of  the  gas  is 
reduced  due  to  the  venting  process.  Equations  for  this  alternate  balloon  state  are: 

D[  Wg|  =  -  (VV  +  VD)  ■!'  SWG  (5a) 

DlTgl  «  (-Q1  -Q7+Q8  +  C)9  -QlO  +  Qll  -  Pg  *  Dj  Wg) /SWG)/(Wg  ■!'  CV) 

(6a) 

{7a) 


DIV]  =  0. 

DlPgl  =  Pg«  (DlWgl/Wg  +  D(Tgl/Tg)  . 


(8a) 
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We  introduced  a  simple  version  of  this  duct  venting  model  in  a  paper  on  bal- 
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loon  design,  and  modified  it  herein  to  work  with  the  dynamic  case.  It  is  a  major 
change  from  all  previous  flight  performance  models.  Those  models  used  a  routine 
called  burping;  the  name  (origin  unknown)  is  somewhat  inelegant,  but  the  routine  is 
mathematically  effective . 

However,  for  the  analysis  of  vertical  motions  of  full  balloons  at  float,  it  is 
important  that  performance  model  outputs  represent  (as  nearly  as  possible)  actual 
flight  performance;  use  of  the  burping  model  precludes  this.  Yet.  even  our  more 
sophisticated  process  does  not  account  for  volume  increases  when  balloons  change 
shape  duo  to  venting  backpressure.  For  a  venting  balloon  carrying  a  payload  less 
than  or  equal  to  its  design  payload,  tile  balloon  shell  is  fully  deployed  and  taut;  the 
volume  cannot  change  perceptibly.  I'he  same  balloon,  carrying  a  significantly 
heavier-than-design  payload,  has  excess  envelope  material  and  excess  potential 
volume  wlieii  it  begins  to  vent.  As  backpressure  is  created  and  rises  due  to  vent* 
ing  excess  gas,  gas  expands  into  the  potential  volume  and  observed  performance 
may  differ  considerably  from  the  model's  output.  The  described  performance  is 
typical  of  fully  tailored,  natural  shape,  free  balloons.  Cylinder  balloon  perform¬ 
ance,  on  the  other  liand,  is  considerably  different  and  more  complex,  but  cylinder 
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balloons  arc  no  longer  used  routinely.  (See  Nishimura  and  Hirosawa.  ) 

2.1.1  THE  EQUATION  OF  MOTION 

D1  v  1  =  (BUOY  -  DRAG  -  WS)/1  (WS  +  CM  *  BUOY)/Gl  , 

where 


V  is  the  vertical  velocity  of  the  system, 

BUOY  is  the  weight  of  the  displaced  ambient  air, 

DR.AG  is  the  aerodynamic  force  resisting  balloon  vertical  motion, 

WS  is  the  system  weight,  including  the  gas, 

CM  is  the  coefficient  of  added  or  virtual  mass, 

G  is  the  gravitational  constant. 
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This  equation  is  essentially  equivalent  to  Eq.  (1)  of  Horn  and  Carlson. 
However,  we  have  adjusted  it  to  enaole  the  use  of  weight  in  place  of  m.ass  as  .a 
primary  dimension.  It  contains  references  to  two  previously  discussed  problem 


42.  Dwyer,  J.  F.  (li)t!2)  Polyethylene  Free  Balloon  Design  From  the  Perspective 

of  User  and  Designer,  AFGL-TR-82-0350.  AD  A1275S3. 

43.  Horn,  W.  J. ,  and  Carlson,  L.  A.  (1983)  THERMTRAJ ;  A  Fortran  Program 

to  Compute  the  Trajectory  and  Gas  Film  Temperatures  of  Zero  PressurF 
Balloons ,  NASA  Contractor  Report  168342”. 
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areas:  added  mass,  expressed  by  the  term  CM  *  BUOY,  and  aerodynamic  drag, 
expressed  by  the  term  DRAG.  As  already  stated,  we  continue  with  the  tradltlojial 
approach  to  added  mass,  but  not,  however,  to  aerodynamic  drag.  The  term  DRAG 
can  be  expanded  as ; 

DRAG  =  0.  5  •-!'  RO  'i'  CD  «  HC  *  v  #  ABS(v)  , 
where 

RO  is  the  density  of  tlie  ambient  air  computed  by  the  subroutine  VIRON, 

CD  is  the  aerodynamic  drag  coefficient, 

HC  is  the  horizontal  gas  bubble  cross  section,  as  defined  in  subi'outine 
MYBLN. 

We  reserve  ti'catment  of  ihe  term  CD  to  Section  3,  where  we  will  provide  specific 
comments  on  the  problems  of  modeling  it. 

2. 1.2  DEBALLASTING  EQUATION 

Tile  deballasting  equation  accounts  for  reductions  in  the  dead  weight  payload 
on  the  balloon.  It  is  included  as  a  differential  equation  (ratlier  than  a  simple  func¬ 
tion  of  time)  as  a  matter  of  choice,  and  because  a  prior  version  of  this  perform¬ 
ance  model  included  provisions  for  a  cryogenic  ballast  system  that  will  be  rein¬ 
stituted  when  sufficient  theoretical  or  practical  interest  arises. 

D(W|  s  -DD  , 

where  W  is  the  weight  suspended  bcMieath  the  balloon  and  DB  is  the  ballast  pouring 
rate . 

2.1.3  A  LTITUDE  EQUATION 

4S 

This  is  identical  to  E'q.  (2)  of  Horn  and  Carlson. 

D[Z|  =  V  , 

where  Z  is  the  altitude  (above  msl)  of  the  system. 

2.  1.4  FILM  TEMPERATURE  MODEL 

D1  TfJ  =  (Q1  +  Q2  +  Q3  +  Q4  -  Q5  +  Q6  +  Q7)/(t:F  >■'  WE)  , 

where 

Tf  is  the  balloon  film  temperature, 

Ql  is  the  rate  of  free  convective  heat  transfer  between  the  gas  and 
balloon  wall, 
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Q2  is  the  rate  of  direct  solar  energy  absorption  by  the  balloon  wall, 

Q3  is  the  rate  of  IR  energy  absorption  by  the  balloon  wall, 

Q4  is  the  rate  of  free  (or  forced)  convective  heat  transfer  between  the 
balloon  wall  and  the  air, 

Q5  is  the  rate  of  IR  energy  emission  by  the  balloon  wall, 

Q6  is  the  rate  of  absorption  of  reflected  solar  energy, 

Q7  is  the  rate  of  radiative  exchange  between  the  gas  and  balloon  wall, 

CF  is  the  specific  heat  of  tl\e  balloon  film, 

WB  is  the  gas  bubble  envelope  weight  computed  from  the  balloon  shape 
parameters  as  determined  by  subroutines  MYBLN  and  NELSON, 

This  temperature  model  and  Eq.  (4)  of  Horn  and  Carlson"*^  are  comparable, 
but  a  significant  difference  exists  between  our  term  WE  and  their  term  MASSP. 

The  latter  term  is  a  constant;  it  accounts  for  the  mass  of  the  entire  balloon  enve¬ 
lope.  On  the  other  hand,  the  term  WE  refers  to  only  part  of  the  fabricated  salloon 
envelope:  that  part  which,  at  a  given  instant,  surrounds  the  bas  bubble.  WE  de¬ 
pends  upon  both  the  balloon  envelope  construction  and  the  degree  of  inflation. 
Conceptually,  WE  was  developed  independently  of,  and  without  recourse  to,  the 
work  of  Emslie.  However,  it  can  be  viewed  and  should  be  viewed,  as  a  more 
accurate  version  of  his  model,  even  though  it  is  a  unique  development. 

As  a  factorable  term  in  the  denominator,  WE  (like  MASSF)  significantly  af¬ 
fects  the  values  compvitod  by  tliis  function,  most  particularly  for  very  large  bal- 
loon.s  and  for  small  fractional  volumes.  One  should  also  note  that  we  have  dropped 
the  product  (Tf  ^  D[WE|/WE),  a  second-order  term  that  results  from  the  fact  that 
WE  is  not  constant. 

Terms  Ql  through  Q7  are  treated  in  more  detail  in  Section  2.2,  the  section 
on  heat  transfer  models. 

2.  1.  n  GAS  WEIGHT  MODEL  (SLACK  BALLOON) 

D(  Wg|  =  -VV  SWG  , 
where 

V\  g  is  the  balloon  ga.s  weight, 

VV  is  the  apex  gas  valve  discharge  rale, 

SWG  is  the  specific  weight  of  the  balloon  gas. 

This  gas  weight  model  is  comparable  with  Eq.  (3)  of  Horn  and  Carlson,^®  but 

it  incorporates  a  mathematical  model  of  the  EV-13  apex  valve  discharge  rates 

3  ^ 

based  on  data  in  an  earlier  report. 
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2.  l.G  GAS  TISMPEHATUKE  MODEL  (SLACK  BALLOON) 

D(  Tg)  =  (-  Ql  -  Q7  +  QO  +  Q9  -  QIO  +  QU  *  SW  «  V*  v)/|  \Vg  *  (CV  +  BCi)]  , 


where 


Q8  is  the  rate  of  absorption  of  direct  solar  energy  by  the  gas, 

Q9  is  the  rate  of  absorption  of  reflected  solar  energy  by  the  gas, 

QIO  is  the  rate  of  emission  of  IR  energy  by  the  gas, 

Qll  is  the  rale  of  absoi'ption  of  IR  energy  by  the  gas, 

SVV  is  the  specific  weight  of  the  ambient  air  computed  by  subroutine  VIRON, 

V  is  the  instantaneous  volume  of  the  gas  bubble, 

CV  is  the  specific  heat  of  the  gas  (at  constant  volume), 

RG  is  the  specific  gas  constant  for  the  Inflatant. 

This  gas  temperature  model  is  comparable  with  Eq.  (5)  of  Horn  and  Carlson^^ 
and  with  Kreith  and  Kreider.  Again,  terms  Q1  and  Q7  through  Qll  are  treated 
in  Section  2.2. 

2.1.7  GAS  VOLUME  MODEL  (SLACK  BALLOON) 

D(  VI  =  V  *  1D(  VVgl  /Wg  t  D[  Tgl  /Tg  +  v/(RA  »  Ta) )  . 
wliere 
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RA  is  the  specific  gas  constant  for  air, 

Ta  is  the  ambient  air  temperature,  compvited  by  subroutine  VIRON. 

Tlu'  gas  volume  model  is  included  as  a  differential  equation,  rather  than  as  a 
definite  func  tion  of  temperature,  pressure,  and  mass,  for  previously  cited  rea- 
soi^s. 

2.1.8  BALLOON  GAS  PRESSURE  MODEL  (SLACK  BALLOON) 

The  nature  of  our  duct  venting  model  requires  that  we  include  balloon  gas 
pies.sure  as  a  differential  equation  for  the  case  of  t)ie  full  balloon;  the  differential 
equation  is  included  in  the  slack  balloon  case  only  because  it  is  required  by  the 
symmetry  of  the  solution  pi'ocess. 

Dl  Pg)  =  ••  SW  V  , 
where  Pg  is  the  balloon  gas  pressure. 

2.  1.  9  GAS  WEIGHT  MODEl,  (FULL  BALl-OON) 

This  gas  weight  model  differs  from  Eq.  (6)  for  slack  balloons  in  that  it  in¬ 
cludes  the  duct  venting  model  previously  discussed  in  Section  2.1.  In  this  I’cgard 
it  differs  also  from  Eq.  (3)  of  Horn  and  Carlson.**^ 


DlWg)  =  -  (VV  +  vn)  f  SVVG  . 
whei't?  VO  is  the  duet  discharge  rate. 


2.  1.  10  GAS  TEMPERATUnE  MODEL  (FULL  BALLOON) 

DI  Tg)  =  (-  Q1  -  Q7  +  QU  +  QO  -  QIO  +  Qll  -  Pg  «  D(  Wg) /SWG)/(Wg  ^  CV)  . 

3 

This  fjas  tcmporaturc  model  is  comparable  to  Eq.  (5)  of  Horn  and  Carlson 
and  with  Krcith  and  Xreidor,  except  that  it  has  been  modifiod  to  reflect  the 
effect  of  venting  excess  gas.  Again,  terms  Ql  and  Q7  through  Qll  are  treated  in 
Section  2.2. 

2.1.11  GAS  VOLUME  MODEL  (FULL  BALLOON) 

■Phe  gas  volume  model  is  included  as  a  differential  equation,  rather  than  as  a 
definite  function  of  leinpcratiiro,  pressure,  and  mass,  for  previotisly  cited  rea¬ 
sons. 

n|Vl  =0  . 

2.  1.  12  DALIOON  GAS  PRESSURE  MODEL  (FULL  BALLOON) 

Ralloon  gas  pressure  is  computed  as  an  essential  factor  in  the  model  of  ducts 
venting  cxc'css  lifting  gas  when  the  volume  of  the  inflatant  tends  to  exceed  the  max¬ 
imum  volume  of  the  balloon, 

DlPgl  r  Pg  -  (D|Wgl/Wg+ D|Tg]/Tg)  . 

2.2  Ileal  Transfer  Moilels 

lloat  energy  added  to  or  lost  by  the  inflatant  is  a  major  factor  in  overall  sys¬ 
tem  performance.  .Added  heat  energy  decreases  the  density  of  the  inflatant,  there¬ 
by  increasing  the  buoyancy.  When  heat  energy  is  lost  the  effect  is  opposite.  To  re- 
il^'ialc,  heal  Iran.sfer  models  (relative  to  polyethylene  balloon  fliglit  analysis)  were 
introduced  by  liull^  pi’ioi'  lo  the  deveJopiuciU  of  flight  performance  models  for  solu¬ 
tion  by  large  scale  digital  computers.  Over  the  last  32  years  these  models  have 

44 

evolved  into  11  elements,  the  last  five  of  which  result  from  efforts  by  Carijon 
to  coi’i'clate  theoretical  and  aclual  fliglit  performance.  Assumptions  underlying 
these  la.st  five  elemeuls  must  still  be  validated,  for  he  notes  that  the  absorption 
and  emission  values  deduced  for  the  inflatant  m  e  not  independent  of  t!ie  balloon 
envelope  materials  used  on  the  flights  from  which  tiie  data  were  obtained. 


44.  (tarlson,  L.  A.  (IhTO)  A  new  thermal  analysis  model  for  high  altitude  bal¬ 
loons,  Proceedings,  Tenth  AFGl.  Scientific  Balloon  ‘Symposium,  Catherine 
L.  Rice.  Ed.  ,  pp.  107 -20(:. 
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Given  our  use  of  FORTRAN  notation  and  nonsubscriptcd  symbols,  and  in  order 
eep  oui'  nonvsnclatnrc  and  symbols  reasonably  close  to  those  of  Carlson  and 
Horn,  wc  express  convective  heat  transfer  between  a  gas  and  some  object  ns  fol¬ 
lows: 

Q  =  CH  dT  *  Area  =(!?■>  Nu/L]  *  dT  *  Area  , 
where 

h  is  the  thermal  conductivity  of  the  gas, 

Nu  is  a  Nussclt  number, 

I.  is  ft  length,  characteristic  of  the  object's  shape  (in  this  application, 
the  maximum  horizontal  diameter  of  tlie  assumed  shape), 

dT  is  the  difference  in  their  temperatures, 

Aren  is  the  involved  surface  area  of  the  object, 

CH  is  the  convective  heftt  transfer  coefficient,  generally  correlated 
with  the  Nusselt  number  by  the  relationship;  ( k  Nu/  . 

2.2. 1  FREE  CONVECTIVE  HEAT  TRANSFER  BETWEEN 
THE  GAS  AND  THE  BALLOON  WALL 

Ql  =  CQl  #  CG  =!■  GNl  *  Nu  *  (Tg  -  TO  «  SA/DM  , 

where 

CQl  is  a  correction  coefficient  to  be  established  on  the  basis  of  exper¬ 
ience  (initially,  CQl  »  1), 

CG  is  the  gas  thermal  conductivity  cucfflciont, 

GNl  is  tlie  Nusselt  number  correction, 

Nu  is  the  Nu.ssolt  number, 

SA  is  the  surface  area  of  the  gas  bubble  shape  comt  uted  by  subroutine 
MYBLN. 

DM  is  the  diameter  of  the  bubble  model  computed  by  subroutine  MYBLN. 

Since  we  still  consider  this  process  as  free,  rather  than  forced  convection, 
we  represent  the  Nusselt  nvimber  in  the  general  form ; 

Nu  =  a  (b  +  c  (Pr  Gr)"*(d)|  , 
where 

a  is  an  arbitrary  constant  (as  are  b,  c,  and  d), 

Pr  is  the  Prandll  number  for  Helium,  0.67, 

Or  is  tl»e  Grashof  number  expressed  as 

G  «  I  (DM  *  SWG/VSG)>!'*2j  ABSd'f  -  Tg)/Tg  , 

SWG  is  the  specific  weight  of  the  gas, 

VSG  is  the  viscosity  of  the  gas, 

ABS(  )  is  the  symbol  for  absolute  value. 
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Ths  viilues  used  for  Hie  constants  ti,  b*  is  and  d  have  not  been  consistent 
during  the  evolution  of  the  flight  porformanve  models  (See  Table  1).  VVe  use 
t-aiTson's  values  In  tixis  proposed  modeli  except  that  we  Include  tl>c  value  of  "a" 
in  *GN1  I  a  torrection  factor.  VVe  believe  that  our  reeoinnu't\dotl  et'nmriij;  will 
require  that  other  eivpirieal  coustnnts  be  eliangeil  alsoi  but  only  f  if  ter  there  has 
been  an  opportunity  tt  correlate  the  results  with  actual  flight  data. 

Table  l.  Nusscll  Nuinbei-  Model  Constants  for  Free  Heat  Ti'nnsfer 
Between  the  Gas  and  tlic  Balloon  Film,  Si.'ov  ing  Similai'itics  and 
Differences,  and  Arranged  Clironologically.  Notes  refleet  valid 
ranges  fori  (1)  lOLT  s  (Pr  >  Gr)  s  lOLTil,  (2)  (Pr  «  Gr)  «  lOilP 
(31  lOi  ;P  <  (Pri'Gr),  (1)  (Pr Gr)  't  l.r**  101:«,  and  (5) 'l.  5* 

10K«  -?  (Pr  Gr) 


b 

C 

d 

Notes 

Germeles*^ 

1 

G 

0.  13 

0.  33 

1 

Fujii  et  al*'' 

1 

0 

0.  35 

0.  25 

Kreith  and  Kreider*® 

3 

0 

0.5P 

0.25 

2 

3 

0 

0. 13 

0.33 

3 

Balls  Creina  et  al^® 

1 

0 

0.  12 

0.3.3 

t  n 

tiarhson  arid  5-Ionv 

2.5 

♦J 

0»  (lO 

A  nr 

•1 

1 

2.. 5 

0 

0.  13 

0.  33 

5 

2.2.2  DIIIECT  SOLAH  KNKBGY  Ai:;J0HPT10N  BY 
the  da i. I. DON  VVAl.l.. 

G2  «  CXJ2  AV  V'V  V  lu:  , 

where 

CQ2  is  .a  correction  coefficient  to  be  established  on  the  basis  of 
experience  (initially,  CQ2  =  1). 

AV  is  the  effective  UV  absorptanee  of  iho  film, 

I'Y  is  tile  effective  UV  flux. 

This  model  o.''  absorption  of  UV  flux  is  e.s.sentially  equivalent  to  that  devel¬ 
oped  by  Germ e lee  .md  refined  by  Kreith  and  Kreider  (see  program  comment  card.s 
for  more  detailed  references).  We  believe  that  the  use  of  a  constant  cross  sec¬ 
tional  shape  is  adequate  for  two  reason.s:  (a)  due  to  uncertaintie.s  in  the  actual 
shape,  and  (b)  due  to  availability  of  the  factor  CQ2  far  making  necossarv,  small 


luljvtBtmoiils.  h»  any  evenl,  tho  sliftpo  of  Iho  cross  sccUon  will  bo  reasonably  con¬ 
stant  for  any  uivon  short  ijcriocl  of  timo. 

2.9.3  III  JiNJiilKIY  A».SOI!PTI()M  BY  TtlK  BALLOON  WALL 
Q3  s  CQ3  0  AH  <<  UZ  <'  SH  *  . 

whore 

CQ3  Is  a  coi'j'ooUon  coefficient  to  be  established  on  the  bii.sis  of 
(initially,  CQ3  •<  I), 

All  is  tlie  offocti\»e  IH  ob.'Jor()tance  of  the  film, 

H?.  1.':  the  Stchn-Tlolt'/mann  constant, 

SR  IB  the  effective  111  receptor  surface  area,  namely 
HC  +  (HC  .  SA)  <•  LOG<O.OOG2C/RE)/LGTKN, 

TI  Is  the  equilibrium  I'adintlon  temperature, 

BE  is  tiie  specific  lift  of  the  inflntant  (a  function  of  altitude), 

LOG(  )  is  Iho  natural  logarithm, 

LGTEN  la  tlie  natural  logarithm  of  10  (converted  to  real  in  program). 

•13 

This  model  is  that  of  Horn  and  Carlson,  '  except  for  the  definition  of  the  cf- 
foctive  .'Surface  area  SH.  This  area  vnric.s  from  the  entire  surface  area  stirround- 
ing  the  gas  bubble  (at  launch)  to  only  the  area  of  the  horizontal  cross -sectior.  of 
the  gas  bubble  (at  altitudes  greater  than  (»0,000  ft.  approximately);  this  Is  con¬ 
sistent  With  Kreith  and  Kreidor.  I'or  altitudes  up  to  tiO.OOO  ft,  the  foregoing 
function  defining  .SR  i.s  ipiUc  etfeelivc. 

2.2.«1  CONVlsCTlVK  HEAT  THANSh'KH  HKTWKKN  THE  BAI..!.,OON 
WALI.  AND  THK  AIR,  KHEK  AND  EOHCED 

q-l  s  <  -Q-l  CA  CiN((  •:  Nu  (Ta  -  Tf)  SA/DM  , 

where 

<.:Ql  is  a  cori'cclion  coefficient  to  be  established  on  the  basis  of 
experience  (initially,  CQ<1  =  I), 

t:A  is  the  air  thermal  cendudivity  coefficient,  computed  by  sub¬ 
routine  VIHON, 

GNk  is  the  Nussell  number  correction,  c  =  2  (free),  f  =  3  (forced). 

2.?...1.  I  .I'h'cc  Convettive  Heat  Transfer  Between  the 
Balloon  Wall  and  Uic  Air 

a  balloon  at  rc.m  (if  sticli  a  i-a.-^e  ever  truly  cxi.sts)  the  Nu.sscll  number  is, 
.I.S  in  .Seenon  9,2.  1.  of  ihc  gcnciMl  form! 

Nu  *  a  I  b  +  e  >  (Pr  '■  Gr)’''’==(d)  1  . 

where  the  Prandtl  number  for  nil'  i.s  0. 72  and  the  Cra.shof  number  is  expre.ssed  in 
this  easf  aS! 
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Gr  »  G  *  [  (DM  *  SW/VS)**2  ]  *  ABS  (Tf  -  Ta)/Ta  , 


where 

SW  is  the  specific  weight  of  air  computed  by  subroutine  VIBON, 

VS  is  the  viscosity  of  air  computed  by  subroutine  VIBON. 

The  values  used  for  the  constants  a,  b,  c,  and  d,  once  again  have  not  been 
consistent  during  the  evolution  of  the  performance  models  (see  Table  2).  As  be¬ 
fore,  we  use  Carlson's  constants  in  the  proposed  model;  except  that  the  value  of 
"a"  is  included  in  the  correction  factor  ’'GN2". 


Table  2.  Nusselt  Number  Model  Constants  for  Free  Heat  Transfer 
Between  the  Balloon  Film  and  the  Air,  Showing  Similarities  and 
Differences,  and  Arranged  Chronologically 


a 

b 

C 

d 

Hall® 

1 

0 

0.  50 

0.25 

Germeles^^ 

1 

0 

0.  13 

0.33 

Fuji!  et  al^^ 

1 

C 

0.65 

0.25 

Kreith  and  Kreider^® 

1 

2 

0.60 

0.25 

1  ft 

Balls  Crema  et  al 

1 

0 

0,56 

0.25 

19 

Carlson  and  Horn 

1 

2 

L. 

0.60 

0.25 

2. 2. 4. 2  Forced  Convective  Heat  Transfer  Between  the  Balloon 
Wall  and  the  Air 

For  balloons  in  motion  the  Nusselt  number  is  usually  expressed  in  the  gen¬ 
eral  form : 

Nu  =  b  +  c  -  (Re)«*(d) 

where  b  is  an  arbitrary  constant  (as  are  c  and  d)  and  Re  is  the  Reynolds  number. 
Here  again,  the  values  of  the  constants  b,  c,  and  d,  have  not  been  consistent  dur¬ 
ing  the  evolution  of  the  performance  models  (see  Table  3).  In  the  proposed  model, 
we  again  use  Carlson's  constants,  except  that  we  ma'<te  no  exception  for  balloons 

Q 

with  maximum  volumes  greater  than  If!  million  ft  .  We  believe  that  any  such 
correction  should  await  evaluation  of  the  effects  on  the  model  output  of  other  non- 
arbitrary  changes. 


32 


i-diW  ■‘W 


• 


HW*»in-iC.n  art  rf«-V^’^'W-^V-^^  TW-TV  ^ry-TVvXVTUYU  YX/VI/'Vy  vuvirvin/TfVi^rv^TVV^CMris?^^ 


Table  3 .  Nusselt  Number  Model  Constants  for  Forced  Heat  Transfer 
Between  the  Balloon  Film  and  the  Air,  Showing  Similarities  and  Dif- 
fei'ences,  and  Arranged  Chronologically.  Notes  reflect  valid  ranges 
for:  (1)  laminar  flow,  (2)  turbulent  flow,  (3)  1.8*  10E3  £  (Re)  ^ 
1.4*10E5,  (4)0.4*  10E5  s  (Re)  £  1.4«10E5.  and  (5)  volumes 
greater  than  19  *  10E6  ft^ 


C 

d 

Notes 

Germeles^^ 

H 

0.  37 

0.60 

Fujii  et.al^^ 

H 

0,  52 

0,50 

1 

0.03 

0.  BO 

2 

1  fi 

Kreith  and  Kreider 

2 

0.30 

0.57 

3 

2 

0.41 

0.55 

4 
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Carlson  and  Horn 

0 

0.  J7 

0.60 

0 

0.74 

0.  60 

5 

_ 1 

2,2.  5  IR  ENERGY  EMISSION  BY  THE  BALLOON  WALL 
Q5  =  CQ5  *  ER  *  BZ  *  SA  *  Tf»*4  , 


CQf)  is  a  correction  coefficient  to  be  established  on  the  basis  of 
experience  (initially,  CQS  »=  1). 

KR  is  the  effective  IR  emissivity  of  the  balloon  wall  film. 

Except  for  the  fact  that  SA  has  been  redefined  in  accordance  with  our  new 

shape  model,  this  is  identical  to  the  respective  heat  transfer  model  of  Horn  and 
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Clarlson. 

2.  2.  (i  IlEFLEt;TE1)  SOLAR  ENERGY  ABSORPTION  BY  THE 
BALLOON  WALL 

1  r 

This  equation  is  based  on  Eq.  (47),  Kreith  and  Kreider.  ’ 

Qfi  =  i;Q(i  *AV  *  (2*HC)  *  GS  *  FK  *  RL*  QA  *  f  (AL,RE)  , 


where 

CQK  is  a  correction  coefficient  to  be  established  on  the  basis  of  ex¬ 
perience  (initially,  CQ6  =  1), 

AV  i.s  the  effective  UV  absorptance  of  the  film  according  to  Carlson, 
GS  is  the  solar  con.stant, 

FF  is  the  directional  reflectivity  factor  according  to  Figure  15  of 
Krcith  and  Kreider, 
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RL  is  the  reflectance  modeled  after  Figure  10,  Kreith  and 
Kreider,  16 

QA  is  the  cosine  of  the  solar  zenith  angle. 

AL  is  the  balloon  altitude, 

KE  is  the  radius  of  the  earth, 
f(  )  is  the  function:  [1.  -  SQRT(AL./RE/2. )) , 

SQRT(  )  is  the  FORTRAN  notation  for  square  root. 

2.2.7  RADIATIVE  EXCHANGE  BETWEEN  THE  INFLATANT  AND 

THE  13AL.LOON  W'ALL 

Q7  =  CQ7  *  El  *  BZ  *  SA  »  (Tg**4  -  Tf**4)  , 
where 

CQ7  is  a  correction  coefficient  to  be  established  on  the  basis  of  ex¬ 
perience  (initially,  CQ7  =  1), 

E!  is  the  coefficient  of  radiative  exchange  between  the  inflatant  and 
the  balloon  wall  film. 

The  relationships  expressed  in  the  models  of  Q7  through  Qll  are  based  on 
Carlson,  but  in  the  default  mode  of  our  propesed  model  we  reject  Carlson's 
hypothesis;  thus,  terms  Q7  through  Qll  are  set  equal  to  0. 

2.2.8  DIRECT  SOI^AR  ENERGY  ABSORPTION  BY  THE  INFLATANT 
Q«  CQC  -  AG  FV  *  HC  , 

where 

CQ8  is  a  <;orre<  tion  coefficient  to  be  established  on  the  basis  of  ex¬ 
perience  (initially,  CQ8  =  1), 

AG  is  the  effective  coefficient  of  absorptivity  of  the  inflatant  in  the  UV. 

2.2.9  REFLE(.;TED  SOLAR  ENERGY  ABSORPTION  BY  THE  INFLATANT 
Q9  =  Cqg  =?  AG  (2  *  HC)  «=  (IS  *  FF  KL»  QA  *  f  (AL,RE)  . 

where 

CQ9  is  a  correction  coefficient  to  be  established  on  the  basis  of  ex¬ 
perience  (initially,  CQ9  =  1). 

2.2.10  IR  ENERGY  EMISSION  BY  THE  INFLATANT 
OlO  =  CQIO  ♦  EG  <=  BZ  *  SA  «  Tg*»4  . 

where 

CQIO  is  a  correction  coefficient  to  be  established  on  the  basis  of  ex¬ 
perience  (initially,  CQIO  >=  1), 

EG  is  the  effective  IR  emissivity  of  the  inflatant. 
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2. 2. 11  IR  ENERGY  ABSORPTION  BY  THE  INPLATANT 

Qil  =  CQll  ■»  EG  *  BZ  1'  SR  *  TI**4  , 
where 

CQll  is  a  correction  coefficient  to  be  established  on  the  basis  of  ex* 
perieuce  (initially,  CQll  *  1). 

19 

This  relationship  is  taken  directly  from  Carlson,  but  it  is  necessary  that  we 
change  his  effective  surface  area  term  to  be  consistent  with  the  assuniption  that 
the  IR  energy  absorbed  is  dependent  upon  altitude  (see  Section  2. 2.  3). 


3.  MODELING  THE  AERODYNAMIC  DRAG  COEFFICIENT 


3.1  General  Conairferationt 

It  is  improbable  that  one  can  develop  an  adequate  model  of  balloon  aerody¬ 
namic  drag  coefficients  by  statistical  means  alone.  Any  reasonably  approximate, 
mathematical  model  must  account  for  a  number  of  hard-to-quantify  phenomena 
and,  therefore,  might  become  quite  complex.  For  example,  the  degree  to  which 
balloon  envelopes  are  deformed  by  dynamic  pressure  due  to  vertical  motion  cer¬ 
tainly  depends  on  shell  stresses  relative  to  film  yield  stresses  (which,  in  turn, 
depend  on  balloon  film  temperature).  However,  one  might  develop  adequate  ap¬ 
proximations  by  analyzing  separately  the  flight  data  for  heavily  loaded,  moderately 
loaded,  and  lightly  loaded  balloon  shells.  On  the  other  hand,  one  also  might  gain 
some  Important  theoretical  Insights  by  studying  relevant  works  on  raindrops  and 
parachutes  (for  example  see  Figure  14).  In  any  case,  one  must  always  temper 
judgement  with  experience  —  consider  the  influence  of  gore  pattern  on  balloon 
ascent  configuration.* 

We  concluded  earlier  that  the  drag  coefficient  must  depend  on  Reynolds  num¬ 
ber,  Froude  number,  and,  most  probably,  fractional  volume  -  three  dimension¬ 
less  parameters.  Given  an  unambiguous  definition  of  the  characteristic  length, 
and  accurate  flight  data  (elapsed  time,  altitude,  atmospheric  temperature,  and  all 
initial  flight  conditions),  we  can  compute  the  average  value  of  each  of  these  pa¬ 
rameters  for  each  increment  of  altitude;  leaving  only  the  corresponding  drag  co¬ 
efficient  to  be  determined. 


’"considerations  such  as  thesy  miglit  have  made  the  statistical  analyses  of  rise 
rates  by  Nolan,  and  Keeney’^  more  valuable  to  balloon  users.  Further,  if  they 
had  used  dimensionless  terms,  they  could  have  reduced  to  a  minimum  the  num¬ 
ber  of  multiple  regression  analysis  (MHA)  terras  -  this  would  have  improved 
the  quality  of  tnclr  predictors. 


Nolan,  G.  F, ,  end  Keeney,  P.  I,.  (1973)  Analysis  of  Factors  Influencint;  Rate 
of  Rise  of  Large  Scientific  Balloons,  AF(;RL-TR-73-07o3,  Ad  778070. 
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REYNOLDS  NUMBER 

Figure  14.  A  Comparison  of  Reynolds  Number  vs  Drag 
Coefficient  for  a  Water  Drop  and  a  Sphere,  According  to 
Gillaspy.  37  By  analogy  the  depicted  relationship  sug¬ 
gests  that,  for  a  given  range  of  relative  envelope  stresses 
and  for  a  given  F'roude  numher,  the  drag  coefficient  for 
an  ascending  balloon  might  be  expected  to  increase  with 
increasing  Reynolds  number  -  in  some  unspecified  range. 

As  in  all  of  Section  3,  wo  assume  that  the  balloons  are 
constructed  from  fully  tailored  gores  (at  least  from  gores 
that  are  fully  tailored  in  the  apex  region) 

Given  the  proposed  comprehensive  flight  performance  model  (including  the 
required  sub-models,  such  as  the  atmosphere  model,  balloon  shape  model,  and 
others),  the  problem  of  determining  a  drag  coefficient  is  straightforward.  We 
need  only  to  establish  an  acceptable  closure  accuracy  for  the  altitude  computation, 
assume  a  drag  coefficient,  and  then  iteratively  solve  the  model  over  each  corre¬ 
sponding  time  interval  until  we  find  a  drag  coefficient  value  for  which  the  altitude 
closure  accuracy  is  satisfied.  Otlier  things  being  equal,  we  will  then  have  a 
reasonably  equally  weighted  set  of  four  dimensionless  variables  for  each  altitude 
increment:  a  drag  coefficient,  Reynolds  number,  r'roude  number,  and  fractional 
volume.  Usually,  however,  this  will  not  be  the  case;  unless,  of  course,  we  are 
unusually  fortunate,  or  we  have  carefully  selected  flight  data  to  account  for  qual¬ 
itatively  or  imprecisely  defined  phenomena  that  significantly  affect  free  balloon 
ascent  rates. 


3>2  S«^eg«tion  of  Datx  for  Analytii 

We  have  already  presented  in  sufficient  detail  the  effects  of  gore  patterns  on 
balloon  ascent  configurations.  We  are  forced  to  conclude  therefrom  that  fully 
tailored  balloons  and  cylinder  type  balloons  will  have  vastly  different  drag  coef¬ 
ficient  models.  Thus,  because  cylinder  balloon  types  (full  cylinder  and  semi- 
cylinder)  are  infrequently  used  and  even  more  infrequently  manufactured  for  new 
programs,  only  fully  tailored  balloons  should  be  considered  in  our  present  effort 
to  model  free-balloon  drag  coefficients. 

Thermodynamic  phenomena  are  covered  by  sub-models  included  in  the  pro¬ 
posed  flight  performance  model.  However,  there  is  some  concern  about  the  accu¬ 
racies  of  the  assumptions  on  which  pre.sent  (and  prior)  heat  transfer  processes 
are  based.  Because  of  this,  we  seek  to  make  our  drag  coefficient  solutions  aS 
independent  as  is  possible  of  thermodynamic  considerations.  The  first  and  obvi¬ 
ous  choice  is  to  use  only  data  from  flights  launched  at  night.  In  this  way,  thermo¬ 
dynamic  inaccuracies  arise  only  from  assumptions  about  the  infrared,  and  con¬ 
ductive  heat  transfer  nnodels.  Errors  in  the  models  dealing  with  direct  and  re¬ 
flected  solar  radiation  are  eliminated. 

A  natural  transition  point  both  in  the  vertical  motion  of  the  balloon  and  in  the 
dynamic,  mechanical  responses  of  the  balloon  envelope  m.aterial,  occurs  near  the 
tropopause.  Both  of  these  results  are  due  to  the  reversal  of  the  ambient  temper¬ 
ature  gradient.  The  former  phenomenon,  the  slowing  of  the  ascent  rate,  is  well 
understood,  and  quite  thoroughly  documented.  The  latter  phenomenon  has  gener¬ 
ally  been  associated  with  balloon  bursts  due  to  cold  brittleness  of  the  polyethylene 
film;  its  implications  with  respect  to  both  subsequent  balloon  failures,  and  altered 
resistance  to  ascent  shape  deformation  are  not  well  understood.  With  respect  to 
shape  deformability,  we  believe  that  the  noted  relaxation  of  the  strain  in  the  en¬ 
velope  material  in  the  crown  of  the  balloon^®  supports  the  contention  that,  for  a 
balloon  ascending  above  the  tropopause,  relative  stress  (actual  stress  divided  by 
yield  stress)  changes.  Thus,  the  shape  can  deform  more  easily.  Therefore,  we 
suggest  that  flight  data  above  and  below  the  tropopause  be  segregated  for  the  pur¬ 
pose  of  drag  coefficient  modeling.  The  least  advantage  of  this  approach  will  be 
the  existence  of  a  logically  distinct  set  of  data  by  which  a  developed  model  may  be 
tested.  (When  the  tropopause  temperature  is  different  than  the  minimum  temper¬ 
ature,  prudent  judgement  is  required. ) 

Finally,  we  recommend  that  flights  be  separ.oted  into  two  other  classes;  heav¬ 
ily  loaded  and  light  or  moderately  loaded  balloons.  To  accomplish  this,  stress 
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indices  such  as  those  used  by  the  NSBF  or  AFGL  (see  Dwyer  )  should  be  adequate. 

4G.  Rand,  .T.  L.  (1902)  Balloon  PHlm  Strain  Measurements,  Workshop  on  Instru¬ 
mentation  and  Technology  for  Scientific  Ballooning,  XXIV  COSPAR  Plenary 
Meeting,  Ottawa,  Canada,  IG  May  -  2  June  1982. 


3.3  Modeling  Using  Multiple  Rcgrcuion  Analysis 

Those  familiar  with  modeling  with  MRA  concede  that  it  is  more  an  art  than  a 
science,  hut,  having  decided  to  use  dimensionless  variables,  we  have  taken  the 
first  and  biggest  step  in  our  analysis  in  a  sound  scientific  manner.  If  we  then 
segregate  the  input  flight  data  as  suggested,  we  will  be  taking  the  second  step  in 
a  sound  scientific  manner;  this  will  help  to  minimize  variations  due  to  other 
factors,  as  previously  noted.  Beyond  this  point  (if  one  is  to  develop  a  practical 
model  using  MRA)  it  appears  that  one  must  rely  on  both  mathematical  art  and  sci¬ 
entific  insight  and,  what  is  equally  as  important,  care  in  the  specification  and 
collection  of  the  input  data. 

4.  REFINING  THE  THERMODYNAMIC  COEFFICIENTS 

It  does  not  appear  that  one  can  easily,  if  at  all,  find  conditions  wherein  the 
aerodynamic  drag  may  be  ignored.  Consequently,  the  aerodynamic  drag  coeffi¬ 
cient  model  must  be  developed  before  it  is  possible  to  introduce  refinements  of 
the  models  for  direct  atid  reflected  solar  radiation.  As  with  development  of  the 
model  for  the  aerodynamic  drag  coefficient,  segregation  of  the  input  data  is  de¬ 
sirable  if  not  imperative.  However,  the  solution  process  does  not  appear  to  be 
so  complicated;  the  general  format  of  the  models  of  the  heat  transfer  processes 
are  fairly  well  established  and  ail  that  remains  is  (hopefully)  to  correct  the  coef¬ 
ficients  . 

5.  CONCLUSIONS 

We  have  laid  a  foundation  for  the  development  of  a  comprehensive  flight  per¬ 
formance  model  based  on  practical  and  theoretical  considerations. 

We  have  proposed  that  the  aerodynamic  drag  coefficient  model  be  based  on 
th. 'ee  dimensionless  variables;  Reynolds  number,  Froude  number,  and  fractional 
volume. 

We  have  shown  that: 

(a)  there  has  been  a  longstanding  and  widespread  error  in  the  definition  of 
the  instantaneous  mass  of  the  balloon  film  involved  in  the  heat  transfer  processes, 

(b)  the  gas  bubble  cannot  be  modeled  realistically  as  a  sphere, 

(c)  the  gas  bubble  is  asymmetrical  except  when  it  is  at  or  very  near  its 
natural  ceiling  altitude, 

(d)  the  actual  gas  bubble  shape,  and,  most  probably  the  added  mass,  is 
directly  related  to  the  type  of  gore  pattern. 


(e)  theory  does  not  support  models  of  drag  coefficient  based  on  Reynolds 
number  only, 

(f )  theory  does  support  the  use  of  the  Froude  number  as  one  of  the  variables 
that  effects  tlie  drag  coefficient  of  a  free  balloon,  and 

(g)  that  fractional  volume  is  a  reasonable  way  to  accommodate  variations  in 
overall  balloon  shape,  consistent  with  the  need  to  specify  a  characteristic  length 
for  use  in  establishing  a  reference  drag  area. 

Finally,  we  have  proposed  a  procedure  for  the  analysis  of  actual  flight  data 
to  enable  the  development  of  a  practical,  but  also  theoretically  sound,  model  of 
the  aerodynamic  drag  coefficient  of  a  zero-pressure,  free  balloon,  and  subse¬ 
quent  refinement  of  the  heat  transfer  models  for  direct  and  reflected  solar  energy. 
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Appendix  A 

Software  for  Devolopinj,  Verifying,  and  Using 
Aerodynamte  and  Tharmodynamie  Constants  and  Models 


Al.  I'liOCUAM  UNDCI) 

Pi'ogi'am  FINDCD  is  a  FORTRAN  coded  set  of  routines  and  models  discussed 
in  the  main  test;  U  is  w ,  I**"’!  to  run  on  an  IBM  PC  ftnd  has  a  compiled  executable 
version.  Tins  progi  nin,  working  on  the  assumptions  tiuit  all  of  tiie 

thermodynamic  models  are  suffioicittly  aeonrate,  collects  for  each  point  in  any 
chosen  flight  profile  the  values  of  drag  eoeffioient  ('U,  lieynolds  Number  RN, 
Fronde  Number  FRN,  and  fractional  volume  VB. 


AM  IViigraiii  logic 

The  logic  by  which  the  values  of  the  terms  CD,  HN,  FRN,  and  V13  arc  deter¬ 
mined  is  shown  in  Figures  Al  luid  A’2.  Fundamentally,  it  is  an  iterative  method  of 
adjusting  the  value  of  CD  between  oaeh  sueoessivc  set  of  points  until  the  actual  and 
computed  altitudes  are  satisfactorily  close  for  the  related  elapsed  time. 

AI.2  Slratpgy 

It  appears  that  the  .accuracy  of  related  aerodynamic  terms  can  be  enhanced  by 
initially  restricting  analyjes  to  flights  launched  and  ascending  in  darUnes.s  -  this 
eliminates  potential  errors  due  to  solar  energy  input  models-  Fventvially  one  will 
have  to  alter  this  program  to  accommodate  deballasting  seo,iiences  «  not  at  all  a 
difficult  task.  This  will  be  required  because  on  most  bigh*altitude  night  flights 
deballasting  is  required  to  maintain  ascent  I'atcs  compatible  with  mission  profiles. 
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(z) 


I - , - J - time  (T) 

Tj  Tj*, 

Figure  Al.  Representative  Segment  of  Vertical  Flight  Profile 
A1.3  All  Aerodynamic  Drag  Coefficient  Model 

The  definition  of  an  adequate  drag  coefficient  model  is  still  to  be  found.  If 
physical  modeling  is  an  art  then  mathematical  modeling,  being  one  step  beyond, 
might  be  considered  a  black  art.  There  is  some  discussion  -  both  in  the  main 
text  and  in  the  notes  imbedded  in  the  program  comments  —  regarding  the  use  of 
MRA  as  a  method;  doubtlessly  the  data  developed  by  this  program  will  suggest 
more  explicit  approaches. 

A2.  PKOGRAM  KROUDK 

Program  FROUDE  is  also  a  set  of  FORTRAN  coded  routines  and  models  dis¬ 
cussed  in  the  main  text;  it  too  is  written  to  run  on  an  IBM  PC  and  has  a  compiled 
executable  version,  FLITE.EXE.  It  assumes  the  existence  of  a  drag  coefficient 
model  in  the  following  format: 

i=20 

CD  =  SUM  lAi-  (FR*=  Cli)  *  {RN**C2i>  *  (VB**C3i>j  . 
i‘l 

This  model  format  was  selected  by  the  author  as  one  easily  adapted  to  MRA  mod¬ 
eling,  one  suitable  for  expressing  simple  series  models,  and  one  with  the  inherent 
capacity  to  e.xpress  quite  complex  relationships.  Like  FINDCD,  this  program  per¬ 
mits  interactive  alteration  of  most  of  the  various  model  coeffi'-ients.  Consequently, 


J- .  r  .•  . 


SET  ALLOWABLE 
ALTITUDE  ERROR,  DZ 


it  can  be  used  for  either  flight  performance  prediction  or  for  analysis  of  perform¬ 
ance  sensitivity  to  changes  in  particular  coefficients. 


A3.  INPUT  FORMATTING 

There  are  two  additional  FORTRAN  programs,  CD.  EXE  and  QCDDATA.EXE, 
written  to  format  and  store  the  required  input  files  for  FLITE.  EXE  and  QCD.EXE. 

A3.1  CU.RXii: 

This  program  supports  only  FLITE.EXE;  it  stores  the  CD  model  coefficients 
and  exponents  in  the  required  format. 


A3.2  0C1)I)ATA,EXB 

This  input  file  formatting  program  supports  both  FLITE.EXE  and  QCD.EXE. 
It  has  one  particularly  interesting  feature;  it  distinguishes  between  radar  flight 
data  and  altitude  translated  from  a  standard  altitude  table.  In  the  latter  case  it 
provides  the  altitude  corrected  for  the  local  atmospheric  temperature  profile  and 
launch  site  pressure. 
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PRDGRCiM  FINDCD 


K  -X-lt « 41  ««  'X  II  «»«-»«««  II 

*  # 

»  PROGRAM:  FINDCD  20  FEB  1985  * 

x  » 

»  THE  EXECUTABLE  VERSION  OF  THIS  PROGRAM  IS  DESIGNATED  AS  'GCD"  # 

41  WHICH  HAS  BEEN  COMPILED  UNDER  MICROSOFT  FORTRAN??  TO  DE  RUN  ON  ♦ 

*  AN  IBM  PC.  » 

*  » 

*  THIS  PROGRAM  IS  USED  TO  DETERMINE  THE  REYNOLDS -NO. ,  FROUDE  NO. .  * 

*  FRACTIONAL  VOLUME  AND  RELATED  AERODYNAMIC  DRAG  COEFFICIENT  FOR  * 

»  WHICH  THE  MODEL -PREDICTED  ASCENT  RATE  AND  ACTUAL  FLIGHT  ASCENT  * 

*  RATE  AGREE  WITHIN  A  GIVEN  TOLERANCE,  OVER  A  REASONABLY  LARGE  # 

*  ALTITUDE  SPAN.  THE  RESULTS  ARE  INTENDED  TO  BE  USED  AS  INPUT  IN  * 

*  A  MULTIPLE  REGRESSION  ANALYSIS  TO  PROVIDE  A  MODEL  OF  THE  DPtAG  * 

*  COEFFICIENT  AS  A  FUNCTION  OF  REYNOLDS  NO.,  FROUDE  NO.,  AND  THE  * 

*  FRACTIONAL  VOLUME.  THE  DF.SIRED  OUTPUT  VARIABLE  VALUES  ARE  FOUND  * 

*  BY  ITERATIVELY  ADJUSTING  THE  ASSUMED  DRAG  COEFFICIENT  UNTIL  * 

*  THE  GIVEN  TOLERANCE  IS  ACHIEVED.  ENHANCED  ACCURACY  IS  ACHIEVED  * 

*  BY  REDUCING  THE  VARIABILITY  OF  FACTORS  INFLUENCING  THE  ASCENT  » 

*  rate;  primarily,  this  IS  EFFECTED  BY  USING  NIGHTTIME  ASCENT  DATA  * 

*  UNAFFECTED  BY  SOLAR  INPUT.  » 

»  « 

»  THE  USE  OF  MRA  TO  ACHIEVE  THE  MODEL  18  ONLY  ONE  APPROACH,  AND  » 

*  ONE  SHOULD  NOTE:  1>  THAT  THE  USE  OF  MRA  IS  VERY  MUCH  AN  ART,  AND  * 

»  2)  THAT  THE  RESULTING  MATHEMATICAL  MODEL  MAY  BE  SIGNIFICANTLY  * 

X  DIFFERENT  THAN  THE  TRUE  PHYSICAL  MODEL.  CAUTION  IS  URGED  IF  THE  * 

*  RESLuTING  FORM  OF  THE  MRA  MODEL  IS  TO  BE  USED  TO  PREDICT  OTHER  » 

»  PHYSICAL  RELATIONSHIPS.  * 

«  » 

*  t-IlK  COMMENTS  ON  SPECIFIC  LINES  OF  CODE,  SEE  PRI3BRAM  FROUDE.  * 

*  » 

*  THIS  PROGRAM  WAS  DEVEL.OPED  AT  THE  AIR  FORCE  GEOPHYSICS  LABORATORY  * 

*  AS  PART  OF  IN-HOUSE  WORK  UNIT  NO.  76591114  * 

*  •» 

COMMON  C A , DT I ,  E ( 30 , 2 ) , PE , PO , RO , SW , T I , T I R I ,  T I RO , TR ,  TROP ,  VS 
DIMENSION  A (5) , B <5) , C <5) , D <8, 5) , Y (8,6> ,0(0,6) , FLY ( 1 00, 2) , FR ( 100, 4) 
CHARACTER*82  FINN, FOUT.FNAME, HEADl ,HEAD2 

DATA  A/5..  14159265-1,  .5,  .2928932,  1.7071068,.  1666666666/, GS/96./ 

DATA  B/0.  ,  2. ,  1 . ,  1 .  ,  2,  / ,C/  .01745329252, .5, .  2928932,  1 . 7071068, .  5/ 
DATA  BZ /3. 699SE-1 0/ , G/-32. 1741/, RE/20855278. / , RA/53. 35/, RG/386. 076/ 
DATA  DTM/ 20 .  / ,  l)TV/0 . 5 / , DT/ 8 . / , JF , LL2, LL3/3* 1 / , LAUNCH , LEAP/2«0/ 

DATA  AC3,AC4,DB,DD0,DDl,T,7T.tTT,VD,W,VV, 07, OB, (29, 010.011/16*0./ 
DATA  ALF/ 1 . 83E-07/ , BET/ . 682/, GAM/1443. /,CV/5B6. 73/ , CF/42B. / 

DATA  AYV/ . 001 / ,  RYV/ .114/, TYV/ . 885/ , AYR/ . 03 1 / , EYR/ . 03 1 / , RYR/ .127/ 
DATA  TYR/ , 042/ , AYRG/ . 0028/ , WOW/ . 0048/ , CM/ .  5/ , VL/ . 0 1 / 

D ATA  CO 1 , CQ2 , CQ3 , CO 4 , COS , Ca6 , CQ7 , CDS , C09 , CO 1 0 , CQ 1 1  / 1 1 * 1 .  / 

DATA  GNl , GN2, GN3/3* 1 .  / 
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HEflDl='  RN  FRD  VB  CD' 

HEAD2='  TIME  ACT  SPEED  ERROR’ 

1  FORMAT (A) 

2  FORMATdX, ’ENTER  NAME  OF  INPUT  FILE;  B:-fil«Bpec. FLY  ’,\) 

3  FORMAT (tX, ’AND  NAME  OF  OUTPUT  FILE;  B:f ilespec. JFD 

4  F0RMAT<13(ElS.a/>,ElS.8) 

5  FORMAT <2F9. 3) 

6  FORMAT <2Fa.0> 

7  F0RMAT(IX,4E13.7) 

8  FORMAT(4(ElS.a/>,ElS.8) 

9  FORMAT (13) 

10  F0RMAT<10(E15.8/),E1S.8) 

11  FORMAT  (A,  13,  A) 

12  FORMAT <6E15. 8) 

13  F0RMAT(1X,2E13.7) 

14  FORMAT (2E13. 7) 

15  FORMAT (3E13, 7) 

17  F0RMAT(1X.3E13.7) 

18  FORMAT <2X, A) 


CALL  INFORM 


OPEN ( 5 , F I LE= ' LPT 1 ’ , FORM- ’ FORMATTED ’ ) 


50  write;*,!)’  DEFAULT  GAS  i,  FILM  VALUES  ?  C  0/1  =  N/Y  3’ 

read;*, 9) I 

IF  ;i.EG!,0)  THEN 
lal 

write;*, 1)’  INSERT  APPROPRIATE  DATA  DISK  IN  B-DRIVE  AND’ 
write;*,  1)’  ENTER  GAS  t’.  FILM  FILE  NAME.  Bt  f  ilespec.GAF’ 
read;*, i)FNAME 

OPEN  <  3 ,  F I  LE“FNAME ,  FOr<M=  ’  FORMATTED  ’  ) 

READ ;3, 4) ALF , BET, GAM, CV, CF.  WOW, AYV. RW. TYV, AYR, EYR, RYR, TYR, AYRG 
close;  3) 

WRITECfi,  DFNAME 

write;*, 1)’  ’ 

END  IF 

IF  ;i.NF..l)  GOTO  50 


60  write;*, D’  DEFAULT  HEAT  XFER  COEFs. ,  VIRTUAL  MASS  COEF.  AND' 
write;*, 1)’  EFFECTIVE  ZERO  ASCENT  RATE  7  C  0/1  =  N/V  3’ 
read;*, 9) I 
IF  lI.EGI.O)  THEN 
1  =  1 

write;*, D’  INSERT  APPROPRIATE  DATA  DISK  IN  B-DRIVE  AND’ 
WRITE.;*,  D’  ENTER  COEFFICIENT  FILE  NAME,  B:  f  i  1  espec .  CMV’ 
RE.AD  ;*,  DFNAME 

fJF'EN  (.3,FILE=FNAME,F0RM=’F0RMATTED’  ) 

READ  < 3 . 1 0 ) CQ 1 , CQ2. COS, C04, CDS , C06 . CD7 , COS , C09 , CO 1 0 . CD 1 1 
READ  ( 3 ,  e )  GN 1 ,  GNr ,  GN3 ,  CM ,  vi. 

CLOSE ;3) 

WRITE  (5,  DFNAME 

write;*. D’  ’ 

END  IF 

IF  ;i.NE.D  GOTO  60 
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WRITE(*,1)’  INSERT  THE  FLIGHT  DATA  DISK  IN  B-DRIVE  AND’ 
WRITE(*,2> 

READ (*, 1) FINN 
WRITE<*, l>’  ' 

WRITE<#,3) 

READ(#,  DFQUT 
WRITE (*,!)’  ' 

OPEN ( 3 , F 1 LE=F 1 NN , FORM= ’FORMATTED’ ) 

READ(3,1)FNAME 

READ (3, 4) SIG, GL, CLEN, AD, VTM. THK, TCP, WB, WP, WT, FL, DBB. SPD. PO 
R_E AD ( 3 , 4 > CD.  ERR , DU3 ,  DCD .  DUS ,  AGAS, ELL, GH,XD,XG,XL,TS,XIN.X!NN 
DO  100  1=1,30 

100  READ(3,5)Ea,  1),E(I,2) 

C  option:  SELECTION  OF  ANALYSIS  CUT-OFF  POINT 
110  INN=INT(X1NN) 

WRITE  <*, ID’  THERE  ARE’ , INN, ’  FLIGHT  DATA  POINTS.  ' 

WRITE(*, D’  ENTER  NUMBER  TO  BE  ANALYZED.’ 

READ (#.9> INN 

IF  (INN.GT.INT(XINN) )  GOTO  110 
DO  120  1=1, INN 

120  READ(3,6)FLY<I,1),FLY<1,2) 

CLOSE (3) 

WRITE<*, D’  ENTER  DIFFERENTIAL  GROWTH  LIMITS;  ALT.  C  3,’ 
WRITE!#, D’  ASCENT  RATE  C  ft/sec  3  AND  TEMP.  C  deg  R  3* 
READ!*. iS)ALlM,VLIM. TLIM 
WRITE!*,  n’  ’ 


WRITE!*,  U'  TO  CHANGE  DT  tma>!3  OR  DI  Ivent]  ENTER  NEW  NON-ZERO’ 
WRITE!*,!)’  VALUE,  OTHERWISE  ZERO.  ENTER  ZERO  OR  ONE  FDR  AGAB. ’ 
WRITE!*.  !>’  DTM,  DTV,  AGAS  I  g,-.s  absorbs  in  IR.  0/1  =  N/Y  3’ 
WRITE!*, M)DTM,bTV,  AGAS 
RE AD ( * , i S ) GG 1 , 6B2 ,  ABAS 
IF  IGGl.GT.O.)  DTM“GG1 
IF  !GG2.GT.O,)  DTV=GQ2 
WRITE!*,!)’  ’ 

C  OPTION:  ORIGINAL  OR  REVISED  DRAG  COEFFICIENT  ESTIMATE  AND  DRAG 
C  COEFFICIENT  INCREMENT  FOR  ITERATIVE  ADJUSTMENT 

WRITE!*,!)’  TO  CHANGE  CD  OR  DCD  ENTER  NEW  VALUE,  OTHERWISE  ZERO 

WRITE!*, 13)CD,DCD 

READ!*, 14)661. GG2 

IF  (GOi.GT.O. )  CD=GGl 

IF  (GG2.GT.0.  )  DCD=GG2 

WRITE!*, 1)'  ' 

WRITE (5. ISIFOUT 
WRITE!5,  D'  ’ 

WRITE!5, !)’  DTM  DTV  AGAS  ALiN  VLIM  TLIM’ 

WRITE <5, 7)DTM. DTV. AGAS. AL:M,VLIM, TLIM 
WRITE!!!,  1)-  ’ 

WRITE  (!?..  1)  ’  CD  bCO- 
WRITEIS, 13)CD.DCD 
WRITEIH,!)'  ’ 


H 


Hniiiiiii 


i.GTEN=LQe<10.  ) 
IN=INT<XtN) 
LL1=1N 
LLO=IN 


EGAS=A8AS 
AYRG=AGASttAyRG 
AYRGlsl . -AYRG 

AV=AVV#  M . +TYV*AYRG1 / ( 1 . -RYV*AYRB1 ) > 
AG=AYRG#TYV/ ( 1 . -RYV»AYRG1 > 

WTX=0. 

SPD=1.&9#SRD 

Dn4=P0 

DD5=P0 


DO  190  1=2,30 
TROP  =  E( 1-1,1) 

IF  <E<I,2).GE.E(1-1,2))  GOTO  200 
190  CONTINUE 
200  TIR0=1.8*E<l,2>-5.55 

T1R1=.74*TIR0 

DTI=-.2&*TIR0/(TRaP-E(l, 1) ) 

IF  <XL.GE.20. )  GOTO  300 
RL=-.0025*XL+. 15 
GOTO  600 

300  IF  <XL.GE.30. )  GOTO  400 
RL=.  1 
GOTO  600 

400  IF  (XL.GE.40.)  GOTO  500 
RL=.005*XL-.OS 
GOTO  600 

500  RL=.0075*XL-.,  15 

600  XD=C(1>*XD 
)(L=C(1  >*XL 

CX=COS ( XL ) *COS ( X») 
3X=SIN(XL)*5IN(XD) 

CALL  IJPSONd  ,S16.ABL) 

GB=CLEN/GL 

CALL  NELSON<UO,GK',&iIG) 
WC=UO*ABL*WOW*TCF-*GL.*I5U 

DO  700  J=1,S 
DO  700  1=1,8 
Y  <l,-.n-o. 
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Y(l.  1)='0. 

Y(2. 1 )-WP+WT 
/(3,1)=ELL 

CALL  VI RON (ELL, BP, DDl .  I)D3. DD5,LL1 , LL3, DD7, DD9,  TK,  1  > 
BE=<1.-,13B135)*SW 

vo=<i.-4Pl>*<wb+wp+wt>  /be 

Y(4, 1  )--=TR 

Y(5. 1)=V0*PE/TR/R6 
Y(A,  n-TR 

Y(7,  1>=:V0 
Y(a,  1)“=PE 

+  f+++4-*++++4+4-ff+++++++  BE(3IN  RUN6E-KUTTA  ITERATION  ++4+++++44.++4+++444+ 

800  DO  1600  KJ=2,5 
KK=LJ--1 

CALL  V I  RON  <  Y  ( 3 ,  KK ) ,  BP,  DD 1 ,  DD3 ,  DD5,  LLl ,  LL3 ,  DD7 ,  DD9 ,  Tk' ,  1  > 

BE-(1.-. 13B185#TR/Y(6,Kk) )*SW 

IF  (WT.NE.WTX)  THEN 
X=LCJG(1.+WB/(WT*WP)) 

CALL  UPS0N(2,X,SXRX) 

CALL  UPS0N(3,sisX,VU> 

VT«VU#QL»»3 

WTX=Wr 

IF  (VT.BT.VTM)  VT=VTM 
ENDIF 

IF  (KK.ea.U  THEN 
VB=Y<7,kK>/VT 

CALL  rlYBLN  ( VB ,  GL ,  GB .  GN ,  DM ,  SB ,  DM ,  RM,  HC ) 

U=l.)0 

IF  (GN.LT.CLEN)  CALL  NELSON <U,GB, SIG> 

WE“SB»  (WB-WC )  +WDW»U#ABL*TCP«Bi.  »GL 

SA=SB*ABL*GL*GL 

SR'-HC 

IF  (BE.GE.0.00626»  SR-HC-^ <HC-SA) *l.OG( .  00626/BE> /L8TEN 
ENDIF 

PG-DH*Bt: 

SWG=Y(5,l<k)/Y(7,K»a 
Vi>G=ALF*V  (6,  KK>  ♦#BET 
Ce=GAM*VSG 
SPEED=Y(l,Kk) 

IF  (LAUNCH. EG. 0)  SPEED=SFD 

SPDSO-SPEED*SPEED 
RN=DN*SW*ABB (SPEED) /VS 
EU0Y=SW*Y(7,KK) 

WB=Y (2, KK) +Y <5. KK) +WB 
F'JRCE=BUav-WS 

DRAG". 5*R0*CD*HC*SPEED#AB5 (SPEED) 


W  -  W  -  W  «  y  n -  w  *■  w 


0(1.  >  =6+  (FOI^GE-DRAG  >  /  (WS-*GH*BUOY » 

n-  (LAUNCH. EQ.O)  D(l,K.n=0. 

D(2,KJ)==‘-DB 

D(3,K0)=V(),I<K> 

»•■♦•++++♦++■*•++•♦+■♦••♦++++++■♦■  f+-  BEGIN  HEAT  TRANSFER  RATES  ++t- «■+•♦•++++++++++++++ 
EYRG=EGAB*3. 42E-  Oi*  ( Y  (6,  KK>  / 1 . a>  **. SIS'.' 

EYRG1  =  1.-EYRG 

eR=E YR* ( 1 . +TYR*E VRG 1/(1, -RYR#EVRQ 1 ) ) 

AR=ER 

EI=EYRG*EYR/ ( 1 . -RYR*EYRG1 ) 

EG=EYRG»TYR/ ( 1 . +RYR*EYRB1 » 

CALL.  PR6R  ( GR ,  .  A7 ,  DM ,  SUG ,  Y  ( 4 ,  KK ) ,  Y  ( 6 ,  KK ) .  VSG ) 

IF  (GP.LT.  15E+07)  GNU=2. 4-. A*GR**0. 25 
IF  (GP.GT. 15E+07)  GNU=. 13*6P»* (1. /3. ) 

Q 1 =CQ 1 4CN 1 »S A* ( Y ( 6 , KK ) - Y  <  4 , KtO  ) *CG«GNU / DM 

DA=SX+CX^^C0a  (C  ( 1 )  *  (GH-XG+T/240. ) ) 

ARMS=  (GP/PO)  iKBGIRT  ( 122B.  6+376750. 44*DA#OA)  -613.a#ABS  (DA) ) 

TRM 1 = . 5*  <  E  X  P ( - . 65* ARMS  >  +E  XP ( - . 075* ARMS ) > 


FV='GB*Ci3*TRNl 
IF  (QA.LT.O.)  THEN 
FV=0. 

OB=-SG>RT  ( 1 .  (RE/  <RE+Y  (3,  KK)  > )  **2> 

IF  ((3B.LT.DA)  THEN 

ZZ=SQKV(1 .  •-rJA*PA)*(RE+Y(3.KK)  l~Rfe 

CALL  V I RON ( Z  Z ,  P AM , DDO ,  DD2 , DD4 , LLO . LL2 , DD6 , DDB , TX , 2 ) 

AM=35. 1*PAM/P0 

TRM2= . 5* ( E  XP ( - . 65*AM) +EXP  < - . 09S*AM ) > 
F(/=>GS*CS*TRM2*TRM2/TRM1 
ENDIF 
END  IF 


02=C(>i2*AV*FV*HC 

Q3=C(J3*AR*eZ*Sr<*TI**4 

IF  (ADS (SPEED) -LT.VL)  THEN 

CALL  PRQR(6P,  .67,DM.SW.Y(4,KK>,TR,VS) 
6N  ,)=GN2*  ( 2.  + ,  6*GP**6 . 25 ) 

ELSE 

6NU=. 37*GN3*RN**0. 6 
ENDIF 

Q4=CQ4*5A*(TR-Y  (4,K'K)  )  *CA*GNU/DM 
n5=^(:G5»ER*DZ*SA*Y  (4  ,  KK)  **4 


.QTifV'.V.V".,  <!>: 


ZS=S7.2957e*ATAN<SG!RT  ( <  t  .~QA»OA>  /  (QA«OA>  >  > 

FFsl. 

IF  <ZR.GT.77.2)  FF*=.09375*ZS-S.437S 
IF  i(FF.EQ.I.».AND.(ZS.6T.25. >)  FF=.0l5?*ZS+.61fe9 

0A=CDA*AV*(2.*HC)»BS»FF#RI.*(l.-S0Rt<Y(3,Kk'>/RE/2. )  >*QA 

IF  lEI.NE.O.)  THEN 

Q7=CG'7»E  1  »B2*  <  Y  <  6 .  k'K )  **4-Y  <  4 ,  KK)  ##4  >  #SA 
Q9=CGB*AQ*FV»HC 

e9«=CG9»AG*  <2.  *Hr ) *GS#FF*RL# ( 1 . -S0KT< Y  <3, KK) /RE/2. > ) #eA 
Q10=CD10*EG»BZ*SA#Y  (<i),KK)*»4 
0 1 i =CQ  n  #EQ*BZ  »SR#T 1 #*4 
END  IF 

++++++4-+++++++-»"*-4-+++-i'+++-»-  END  HEAT  TRANSFER  RATES  +++++++++++++'fH+++++++ 

D ( 4 , K J )  =  <  0 1 +02+03+04-05+06+07 ) /CF/WE 

IF  ( <VV.QT.O. )  .AND. IKK.EQ. 1) >  THEN 
IF  (PG.GE. 0.936)  W=.559-PG/8 
IF  <FG.LT. 0.936)  VY=.72222#S0RT( <1.872-PG>»PG) 

W=6. 127*YV*S0RT  (SW/SWG-t.  > 

END  IF 

IFIVD.LE.O.)  1VENT=0 

IF  <IVENT.EG!,0)  THEN 
D(5,Ka)=-VV»SWS 

D ( 6 , K J )  =  ( 08+09-0 1-07-010+011 -SW# Y  <  7 , kK  >  *D (3 . K J ) ) / Y ( 5 , KK' ) / <CV+RB  > 
D<7,KJ)-Y(7,KK)+(D<5,kJ)/Y<5,KK)+D(6,ka> /V(6,k:K)+D<3,KJ)/RA/TR) 
D<a,Ka)=-SW*D<3,KJ) 

ELSE 

IF  (LEAF'.EO.O)  DT=DTV 

VD= AI)*B0RT  ( 2 .  *G+ APS  ( Y  ( B ,  Kk )  -PE )  /  SWG ) 

D <S, KO > =- <VV+VD ) +SWG 

D<6.KJ>=(08+R9-Gl-07-Q10+Gll-PE*D<5.KJ> /SWG) /Y<5,Kk> /CV 
D(7.K0)=0. 

D  (S.KJ  >  =Y  <B,  kk)  *  <D  <5,  KJ )  /Y  <5,kk)  +D<6,KJ  >  /Y  (6,I<K)  > 

END  IF 

DO  1500  ki=i,a 

00= A  < K J  >  *  <  D  <k I ,  k.) )  -B  < kvl )  #0  (k I ,  kk ) ) 

Y  ( K I ,  K J )  =  Y  ( K I ,  kk:  >  +DT+00 
1 500  G ( K i , KJ ) =Q ( K I , kk ) +3. *00 -C ( k J ) *D  <k I , KO  > 

1600  CONTINUE 

++.^+++++++++++++++-*-+4-++  end  PUNGE-KUTTA  ITERATION  +++++++++.-+++++++++++ 

IF  (V<7,5) .GT.VT)  THEN 
VT=Y(7,  1) 

IVENT=:1 
GOTO  800 
END  IF 


IF  (IVENT.EO.O)  THEN 

IF  <ABS(Y(l,S)-V<l,i>>.QT.VLIM)  GOTO  17oO 
IF  <ABS<Y<3,5)-Y(3,l>>.GT.ALIM>  GOTO  1700 
IF  <ABS<Y<4,S)~Y<<;.1))  .GT.TLIM)  GOTO  1700 
IF  (ABS(Y(6,5)~Y(6,t)).LE.TLIH>  GOTO  1750 
1700  DT*DT/2. 

IF  (DT.LT.O.S)  THEN 
WR1TE<*,1>>  BAD  EXIT  t  DT  <  0.5  3’ 

STOP 
END  IF 
GOTO  SOO 
1750  DTX=DTM 

DO  1800  1=1,7 

IF  <Y(l,5)#bTX.LT.ALIM)  GOTO  1900 
1800  DTX=DTX-2. 

1900  CONTINUE 
END  IF 

C  PRE-LAUNCH  STABILIZATION  CLOCK  TIMER 
IF  (LAUNCH. EO.O)  THEN 
TT=TT+DT 
DT=DTX 

IF  (TT.OE.TS)  THEN 
LAUNCH=1 

TT=AN1NT<FLY(2, 1)-FLY<1.1> ) 

DAL=ERR#<FLY(2,2)-FLY(1.2> ) 

IF  (DAL.GT. 100. )  DAL=10b. 

COLL  VIR0N(V(3,5).BP,DD1,DD3.DD5,LL1,LL3,DD7,DD9,TK, 1> 
VB=Y(7,KK)/VT 

CALL  MYBLN<VB.6L,GB,6N,DH,SB,DM.RM,HC> 

WV=SW/VS 
GOTO  2100 
END  IF 
GOTO  2300 
ENDIF 

C  ELAPSED  PLIGHT  TIMER 
T=T+DT 

C  INTEGRATION  INCREMENT  CONTROL  TO  ENSURE  THAT  DT  IS  LESS  THAN  OR 

C  EDUAL  TO  THE  ACTUAL  TIME  IN  CORRESPONDING  FLIGHT  INCREMENT 
IF  <LEAP.E(3.0)  THEN 
TT=AN1NT(TT-DT) 

IF  (TT.LE.DT)  THEN 
DT=TT 
LEAP=1 
GCno  2300 
ENDIF 
GOTO  2300 
ENDIF 


*  VIST *■*«««  •  rwt«rwvn«  i*n*i  rw-T«>»*  •u~i*nv»-ivvjMnrv^v^«>i 


CC  BINARY  CHOP  ROUTINE  TO  CHECK  ALTITUDE  CONVERGENCE  AND  TO  CONTROL 
CC  NEW  DRAB  COEFFICIENT  ESTIMATE 

C  COMPUTES  ALTITUDE  CLOSURE 

X«Y<3,S)-FLY<JF+1,2) 

C  ALTITUDE  NOT  WITHIN  LIMITS 

IF  <ABS<X>.GT.DAL)  THEN 
IF  IX.BT.O.)  THEN 
AC3al. 

SGN-l . 

GOTO  19S0 
END  IF 
ACA«1. 

SBN=-1. 

1950  DCD=DCD/ (AC3+AC4) 

CD«CD-)-&GN«DCD 

*  RESETS  T.  TT.  DT,  LEAP  AND  RUNGE-KUTTA  VARIABLES  FOR  INTERVAL 

»  RECOMPUTATION 

TT=AN1NT  <FLY< JF+1 , 1 ) ~FLY< JF, 1 ) ) 

DT=DTM 

LEAP=0 

TsTTT 

DO  2000  1*1,9 
VIIjDsYlIjib) 

2000  (J(I,1)*G<I,6) 

GOTO  800 
END  IF 

C  ALTITUDE  WITHIN  LIMITS 
IFIT=X 
AC3«0. 

AC4=0. 

0CD».O96*CD 

CC  END  OF  ROUTINE 

CC  ROUTINE  TO  COMPUTE  AND  STORE  OUTPUT 

SPEED=ABS ( ( Y (3. 5) -XXAL) / (T-TTT) ) 

VB=Y(7,KK> /VT 

CALL  MYBLN ( VB , GL , QB , GN , DH . SB , DM , RM,  HC  > 

CAUL  V I  RON  ( Y  ( 3 , 5  > ,  BP  ,  DD 1 , 003  ,  DD5  ,  LL I ,  LL3 ,  DD7  ,  DD9  ,  TK. ,  1 ) 

WV=SW/VS 

C  DETERMINES  AVERAGE  V'ALUES  IN  INTERVAL 

XXVB*<XXVB+VB)/2. 

XXHC=(XXHC'*-HC)/2. 

XXRM=(XXRM'*-RM)/2. 

XXWV*<XXWV+WV>/2. 

RN=2. *XXRM*XXWV*SPEED 
FRD=  <SPEED**2>  t  (2. »XXRM*B) 


rt  o«>«>rN-i:rv  tru  n  isnw^  TV•;T^^•J^^vTW'i■^V^TV^^TV▼l^'fVyV^< 


C  SIOREin  AVERftGE  VALUES  IN  INTERVAL 
FR(aF-,  1)=RN 
FRfjr,2)=FRD 
FR(^^I^3)=XXVB 
FR<JF,4)==CD 

C  WRITES  OUTPUT 
WRITE(5,  DHEADl 
WRITE<S,7)RN,FRD,XXVB,CD 
WRITEIS, I)HEAD2 
WRITE<S,7)T.Y(3.5).SPEED.2F1T 
WRITE(5,  U’  ’ 

CC  ENp  OK  ROUTINE 

C  TERMINATION  CHECK 
JF“JK+1 

IF  (aF.ED. INN)  GOTO  2500 

C  ROUTINE  TO  SET  INITIAL  VALUES  FOR  NEXT  ITERATION  INTERVAL 
TT=ANINT<FLY<JF+1, l)-FLY<aF, 1) ) 

UAL^ERR*  (FLY  <>)F+ 1 , 2)  -FLY  <JF,  2) ) 

IF  (DAL.QT. 100. )  DAL=100. 

DT=--DTN 

LEAP=--0 

TTT=T 

2100  XXAI=Y<‘1.,5) 

XXHC=HC 

XXRM«RM 

XXVB-VB 

XXWVaWV 

DO  '.’200  1-1,0 

Y(I,6>==V(I.!i) 

2200  0(1,6)  «U( 1, 5) 

2300  DO  2<1(T0  1  =  1,8 
Yd,  1  )=V(I,S) 

2400  Q(I,  1)=(?(I,5) 

GOTO  800 

n  TERMINATION  SEOUENCE;  WRITES  OlFRuT  FILES  AND  CLOSES  FILES 

:?r.oo  CLOSE  (5) 

OPEN ( 4 , F ILE=FOUT , STATUS^' NEW' , FORM= ' FORMATTED ' ) 

WRITE(4,  DFNAME 
INN=INN~1 
WRITE (4,9) INN 
DO  2550  1=1, INN 

2550  WRITr£<4,  12)FR(I,  1),FR(I,2>,FR(I,3>,FR(I,4) 

CLOSE (4) 

STOP 

END 

===:===-.:=;i:==:=:r^======-..======-.=:=u=r=PR06RAM  END  ====«=j=T,=:=u:=====!==r-.=^«s==;t,ss:== 

FOR  INFORMATION  ON  SU8R0UTINES,  SEE  PROGRAM  FROUDE, 
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s  •-  .V 


••.'Wn-T^.y’oVs- 


•JiUlliKOUl  IME  PRGR(X0,XKX:;,X3,X<1,X'5,X6) 
XO--!:M7<1xX1^*X2!'  t  (X2«*X'5/X6)^»2>  «AB5  (X^i-XS) /xn 
EMI' 


SUDROUTt NE  MYDUN (VB. GL, C5&.  6N. DM,  SB,  DM,  f^M,HC> 
CALL  BRALIN(5,VP,'2D> 

CALL  BRAUN (I,  ZD. TH) 

CALL  BRAIJN(2,2E,RB) 

CALL  D'RAUN(3,ZB,aO) 

Tl  l"0 . 0 1  V'»532925£»TH 
GB-l . Z  < 1 . -UO+RB/SIN  <TH) ) 

GM-Ol  «GB 

DM''..'.03l2»BN 

RN=l)M/2, 

VC=-3.  I<1i3927*<  <GN*Rt))i*»3)/3./TAN(TH) 

ZM'-.ZilOO^tON 

DH“ZM-ZB*GN 

D> l.-GB 

CALL  BRAUN<'»,0,00) 

HC*=3. 1415?27»RM*»2 

RB“l.-00 

END 


SUBROUTINE  UPSON <J.X,Y> 

DIMENSION  IIP(B,3) 

DATA  UP/  1 , 23S««47BS.  .  3SaPS?42S4,  -  .049794241 ,  -  .  374M?a6S, 
•  .  1431693335,.  l*»/..lU’0O2,-.  1170171749,4.  ,0... 375469014941, 
.  O4140447394-32, . 0 1 79374S40922,  - 1 . Q3972293Q64I--03, 

*2. 2592421 7923C-Ct4, -I .444170710145-05,6. , . 12605506, 

» .  04 1 374 1 0905 1,-0.  .43907397 454E -03 .  - .  090740353 . 

«.  000007007 233  4 . 04002202 1 7E-03 ,  - .  0 1 0  U  07457943 . 4 .  / 

Y-UP  Cl , vl ) 

Ill-lN(<UP(Q,vl)) 

DO  100  i«i,in 
100  V’Y+UP(I 
END 


SUBROUTINE  NELSON (Y, R. X ) 

DIMENSION  W<7,4) 

DATA  W/-. 195390893, -,9529395, .819023839,  .581211894, 

*-l .633868387,1. 118915175, -.248042045,.  19664587,4.64442895. 

♦-3. 2484 1125,-3. 472893627, 7. 232507341, -4. 25933462, . 7859002, 

»3 . 383282482 , -6 . 630204025 , 4 . 076038773 ,5.56431 556 ,-10. 0539 1 7 . 

*5 . 40409 12,-. 859 1 52769 , -2 . 3999589 . 2 . 9527 1 P  <  , - K  645782348 , 

70639 1858, 4 . 499006036, -2. 28485846, . 325235255/ 

9=0. 

XX==1.-R 

DO  200  J=l,4  . 

ZX=W(1,J) 

DO  lO0'l==2,7 

ZX  =  ZX+-W< »,  J)*X**a-l) 

Y=Y+ZX»XX**<J-1 ) 

Y=1  .-Y 
END 

SUBROUTINE  BRAUN<J,X,Y) 

DIMENSION  PR(10.15) 

DATA  8R/56. 151674,36.869415,-963.384153.4340.60078.-23751. 1652. 
*60663. 6084 , -55452. 3401 , 0. , 2.  , 7. ,0. 00012929. 1.491118, 

*0. 89926756885, -12. 4916446722, 15.759250787, 10.1738213142, 

*-26. 202606089,0. , 2. ,  7. , .000128, 1.791098667,0.901115015.  ’ 

*-12, 302923389, 23. 7:' 1897323, -8.8057025141 , -8.2295252039, 

*0. ,-0,00037994265363,0.019166414526, 1.87893200657, 

*1'.  I  06906  70257,  1 .59554475729,-5.  19''400V4837. 2.5901953799,0.  ,2.  , 

47. .  1. , - .00533990 13856, -.246383 177,0. .0. ,0. ,0. , . 999685,8^ ,0. , 

*13. 475326 39 , -13. 453150437, 0. , O. , 0. . 0.  , 0. . . 998704865, 1 . , 2 '  , 

*5. 701 ! 6145,-5. 668903443,0. ,0. ,0. ,0. .0. , .996376, 1. . 2. , 

*3. 17  795699766.-3. 13250803076,0. .0. ,0. ,0. ,0. , .99216, 1 . ,2 
♦  15.  7”5‘v;415,  -67.659164535(,90.6791946''9,-3I:!.  706776  ;.93, 0  .  ,0.  , 

*,96  ,4.  .  .924773778V6,  1.83075490637,-13,8588229434,0.  ,0., 

,  .83035916.0.  ,0.  ,  .31091805100.  .27740682536.-2.63309861357, 

*7.  4<.'6l  ■•2484  ,-11. 064485699,8. 48105265334  .  -2.  69143362923, 

4. 408.2. .  7. . .340858423189  18492521608,0. .0. ,0. ,0. ,0. ,.27218,1. . 

.  362984275675.  -.437916794366, 1.02'561964267.  -i  .44464774487, 

'.,0,0... 0500769.2. , 4. ,5. 3392134074,-28. 1122836054 , 

■'37.0145051  127,0.  ,0.  .O.  .0. , . OO5240.69,  9.  ,0.  ,  .  383179458027, 

*-2.  '.^07  204088 1 1 .0.  ,0.  ,0.  ,0.  .O.  ,  .00015.  1.  ,  2.  / 

IF  UI.LI.5)  GOTO  200 

DO  I'M.i  .1=5,15 

IF  (X.8T  .&R(a,.J)  )  goto  200 

OaMTINUE 

LI  -INT  (Brv(9.  J)  ) 

MM-INT (BRf 10,3) ) 

IF  (l.L.ED.3)  THEN 

Y=^  (“08(2,  T  )  -SORT  (PR  (2,  J)  *BR  (2,  J)  -4  .  *BR  <3,  ,1  >  *  (BR  ( 1 , 0  )  -  X  )  )  > 

Y=Y/ (2.*&R(3, J) ) 

RETURN 
END  IF 
V=BR(1 , J  > 

DO  300  1=2. MM 


»/v-  ■4f«  I  l  V/K$HH*\t>M\PtH\hS\. 


300  Y=Y+BR(I, J>*X*»(I-1) 

END 

SUBROUT  I  N’t  V I  RON  ( X  A ,  PO,  EF .  EG  ,  PX ,  J  X ,  Jf  V ,  RR .  DY ,  TK ,  NN  ? 

COMMON  CA . DT I ,  E ( 30, 2 ) , PE , PO,  RD, SW, T I , T IR 1 , T I RO, TR . TROP, VS 
DATA  CC/34. 1A3195/’ 

X=. 304S037*XA 

X=6356. 766*X/ (635676A. +X > 

IF  (<X.LT.E(JX,1>).0R. (X.GE.EF))  THEN 
IF  <X.LT.E(3X, 1) )  THEN 
JV=1 
PX=PO 
GOTO  200 
END  IF 

100  IF  (ABSIDY) .3E.0.01)  PX=PX« < (E< JX , 2> /EG) ** (CC/RR) ) 

IF  (ABS(DY) .LT.0.01)  PX=PX*EXP(-CC*DX/E ( JX, 2) > 

200  JX-JV 

JV=JVtl 
EF-E(JV.  IJ 
EG=E(JV,2) 

DY=EG-E<JX,2) 

DX=EF-E<JXj 1) 

RR='-DY/DX 

IF  (X.GE.EF)  GOTO  100 
ENDIF 

DX=X-E(JX,1) 

TK=E(0X,2)+Rr<*DX 

IF  (ABG(DY)  .GE.0.01)  PQ^^PX*  ( (E  <  JX,2> /TK)  *#  <CC/RR)  > 

IF  (ABS(I)Y)  .LT.O.Ol)  PD=PX*EXP (-CC*DX/E (OX ,  2) )  . 

IF  (NN.EG.l)  THEN 
TR=I.0*rK 
PE-2.08B58«PQ 
SW’-.0217<1Q/|!(-P0/TK 
R0=SW/32.  1741 

VS=7.  3023327!:-07*  (TR**I.S1  /  <TR+1V8.  72» 
CA=3.30517E-04<KTK«*1.5)/(TK+245.4*(10**  (--12./TK>  )  ) 
n-TlRl 

IF  (X.LT.TROP)  TI«7IR0+DTI*(X-E(1, 1>) 

ENDIF 

END 


SUBROUTINE  INFORM 
C  3  SEP  19S4 

DIMENSION  LINE (23) 

CHARACTEf?*72  LINE 
CHARACTER* 15  FINN. FNAME 

1  FORMAT (13) 

2  FORMAT (A) 

3  FORMAT ( 4 X,  A) 

4  FORMAT  (A.  \.) 

FINN=^'QCD.FAX’ 

L=0 

OPEN ( 1 , F ILE=FINN, FORM-' FORMATTED' ) 
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READ(1,2>FNAME 
READ (1,1)1 

IF  <I.NE.678)  60  TO  9 

l:  NUMBER  OF  LINES  OF  TEXT  UNIT  TO  BE  DISPLAYED 
READCl, 1) I 
IF  (I.EQ.O)  GO  TO  9 
N=0 

IF  <I.LE.20)  N=<24-I)/2 
CALL  CLEAR (25) 

DO  7  K=1,I 

7  READ (1.2) LINE (K) 

DO  8  F=1,I 

8  WRrTE(*,3).LINE<K) 

IF  (N.GT.O)  CALL  CLEAR (N) 

PAUSE 

IF  (L.E(3.  1  )  GOTO  6 

WR1TE(*,4)'  DISPLAY  PROGRAM  NOTES  ?  E  0/1  =  N/Y  1 

READ(*,1)L 

IF  (L.EQ. 1 )  GOTO  6 

9  CLOSE (1) 

CALL  CLEAR (25) 

EMD 

SUBROUTINE  CLEAR (J) 

1  FORMAT (A) 

DO  2  1  =  1,  J 


PROC-iRAM  FROUDE 

CCDMMON  f:  A ,  DT I ,  E  <  30 ,  5  >  ,  PE ,  PO ,  RO,  SW ,  T I ,  T I  ft  1 ,  T I  F!0 ,  TR ,  TROP ,  VS 
DIMENSION  Adi)  ,  El (5)  „ C  (5) , D  (8,5) ,  Y (B, 5) , G!  (0, 5)  ,  RFV <4,  20> 

CHP'lRAC:TEft*32  FINN,  FOOT, FNAME, FMOD,  FXMOD,  HEADER 

*  * 

*  PROGRAM  NAME:  FROUDE  .  .  20  FEB  1905  * 

*  * 

*  IMF  EXEOUTABLE  VERSiON  OF  THIS  PROGRAM  IS  DESIGNATED  AS  ’FLITE'  * 

»  WHICH  HAS  BFCN  COMPILED  UNDER  MICROSOFT  FORTRAN??  TO  BE  RUN  ON  # 

*  AW  IBM  PC.  * 

*  * 

*  THE  DRAG  COEFFICIENT  AG  A  FUNCTION  OF  REYNOLDS  NO.,  FROUDE  NO.,  # 

*  AND  FRACTIONAL  VOLUME  IS  AS  FOLLOWS*.  * 

*  CD  =  SUM  CRFVd.  J)*(  (RN**RFV<2,,7)  )*(FR**RFV(3,J)  )*(VEnn*RFV(4,  J) )  3  * 

*  OVER  THE  RANGE  OF  TO  JSEG,  WHERE  1 . LE.  USED.  LE.  20,  AND  * 

«  WHERE:  RN  IS  THE  r<EYNDLD3  NO.  ■* 

»  FR  IS  THE  FROUDE  NO.  * 

n  VB  IS  THE  INSTANTANEOUS  FRACTIONAL  VOLUME  OF  THE  BALLOON  * 

*  RTVd.J)  ARE  CONSTANTS  DETERMINED  BY  MULTIPLE  REGRESSION  * 

n  ANALYSIS  OF  ACTUAI.  FLIGHT  DATA.  * 

*  « 

«  I  HE  SHAPE  OF  THE  BALLOON  IS  ASSUMED  TO  BE  THAT  OF  TIE  SIMPLIFIED  * 

*■  PARAC;HUIE- SHAPE  MODEL  DESCRIBED  IN  REPORT  NO.  AFGL-Tft-80-0277.  * 

*■  # 

*  Tilts  FTiOGRAM  WAG  DEVELOPED  AT  THE  AIR  FORCE  GEOPHYSICS  LABORATORY  * 

*  AS  F'ART  OF  IN  HOUSE  WORK  UNIT  NO.  76591114.  * 


**»«#»«««'«.'»«»•*«'»*»*»«'«***•««*«*«'*«*«»*»«»«•»***«**»'»*«'***»*•****#»•*»*« 
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LV-VTfi  IV'Z.  i  ■11S92i)S4,  .  5, .  2928932,  1 . 70?  lO.'.H,  .  1 666<>6i)6A6/ ,  GS/96.  / 

Di-.TA  B/0 .  . ,  K  .  1 .  ,  2 -  / „ C/ .  01 745329252,  .  5,  .  292B932,  1 . 7071 06B,  .3/ 
D7VrA  EiZ/3.6995Er-10/,G/32.  1741/,RE/208ij527a.  / , RA/53. 35/ ,RG/3a6. 07Z>/ 
DATA  DTM/20.  /  ,DTO/a.  5/,I)T/8.  /  ,LL2,l.L3/2*l  /  ,  LAUNCH/0/ 

DATA  AC;',  AC4,i)D,DOO,DDl,T,TT,TTT,VD.  VT,W,  07, 08.  (39,010,,  01 1/16*0.  / 
DATA  AI.F/  1 . 03tf--07/ ,  BET/.682/,bAM/  1443.  /,  00/586.  73/,  CF/428.  / 

DA TA  AV V/ . 00 1 / , RYV/ .114/, TYV/ . 885/ , AYR/ . 03 1 / , EYR/ .031/, RYR/ , 127/ 
DATA  T YR / . 84  2 / , A YRB / . 0028 / . WOW/ . 004 0 / , CM / . 5 / , VL / . 0 1 / 

DATA  CO  1 , C(a2 , C03 , CQ4 , COS. CQ6 , C07 , COB , CQ9 ,C0l6, CO  11/11*1./ 

DATA  CD,GNl,GN2,f3N3/4*l./ 

HEADER^'  TIME  ALT.  SPEED  RAD  CD  FR  REN 

+  YB’ 


1  FORMAT (A) 

2  FORMAT!  n, 'ENTER  NAME  OF  INP’UT  FILE!  B: lespec.FLY 

3  FORMAT (IX. 'AND  NAME  OF  OUTPUT  FILE;  B: ti 1 eapec . FLT  ’,V) 

4  F0RMAT(13(ElS.a/).F.15.8> 

5  FORMAT  (2FS'.  3) 

6  FORMAT (2Fa.O) 

7  rdRMAT(lX,6E13.7) 

8  FaRMAT(4(E15.8/),E15.8) 

9  FORMAT ( I 3> 

1 0  FORMAT ( 1 0 (E 1 5 . 8/ )  ,  E 1 5. 8) 

13  FORMAT (E 15.0) 

14  FORMAT (1X.3E13.7) 

15  FORMAT (3E 13,  7) 

16  FaRMAT(lX,F9.  1 , F8.0, F6. 0, F6. 1,E9.3,F7.4, ,2E9.3) 

10  FORMAT ( 2 X, A) 

C  CAll  S  TO  THE  SCREEN  A  SERIES  OF  FACTS  AND  NOTES  REGARDING  THIS 

C;  PROGRAM.  THE  FILE  ADDRESSED  CAN  BE  USER  AUGMENTED  OR  UPDATED. 

CALI.  INFORM 

C  OPENS  THE  PRINTER  AS  THE  OUTPUT  FILE 

OPEN  ( 5 ,  F 1 LE'’ '  Lf^T  1  ’ ,  FORM=  ’  FORMAT  TED’  > 

C  OPTIONS:  VALUES  OF  SAS  it  FILM  CONSTANTS 

50  WRITE!*. 1)'  DEFAULT  GAS  ?<  FILM  VALUES  ?  C  0/1  =  N/Y  3’ 

READ!*, 9) I 
IF  (I.EQ.O)  THEN 
1  =  1 

WRITE!*. 1) '  INSERT  APPROPRIATE  DATA  DISK  IN  B-DRIVE  AND’ 
WRITE!*, 1)’  ENTER  FILM/GAS  FILE  NAME,  B: -f i 1 espec.BAF’ 

READ!*,  UFNAME 

OPEN !3,FILE=FNAME, FORM=’ FORMATTED’ ) 

READ ( 3 , 4 ) ALF , BET  ,  G AM, CV , CF , WOW , AYV,  R YV , T YV , A YR , E YR , RYR , TYR, AYRQ 
CLOSE (3) 

WRITE <5,  UFNAME 
END  IF 

IF  (l.NE.U  GOTO  50 
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OPTIONS:  HEAT  TRANSFER  COEFFFICIENTS,  VIRTUAL  MASS  COEFFICIENT, 
AND  EFFECTIVE  ZERO  ASCENT  RATE 

WRITE!#, 1)'  DEFAULT  HEAT  XFER  COEFb, ,  VIRTUAL  MASS  COEF.  AND’ 
WRITE <#,!)’  EFFECTIVE  ZERO  ASCENT  RATE  ?  C  0/1  =  N/Y  3’ 

READ!*, 9) I 
IF  il.EQ.O)  THEN 
1  =  1 

WRITE!#,!)’  INSERT  APPROPRIATE  DATA  DISK  IN  B-DRIVE  AND' 
WRITE!#,!)’  ENTER  COEFFICIENT  FILE  NAME,  B: f i 1 eepec . CMV’ 
READ!*,  DFNAME 

QPEN!3,FILE=FNAME.F0RM=’F0RMATTED' ) 

READ !3, 10) C«1 , C02. CD3, C04. CQ5, COS,  CP?, COB, CQ9, COlO,  COl 1 
RE AD ! 3 , 8 ) GN 1 , GN2 , SN3 , CM , VL 
CLOSE! 3) 

WRITE  (5.  DFNAME 
END  IF 

IF  II.NE.l)  GOTO  60 
INPUT:  FLIGHT  DATA 

WRITE!#,!)’  INSERT  THE  FLIGHT  DATA  DISK  IN  B-DRlVE  AMD’ 
WRITE!*, 2) 

READ!*, i)FINN 
WRITE!*.!)’  ’ 

WRITE!*, 3) 

READ!*,  DFOUT 
WRITE!*, !>’  ’ 

OPEN ! 3 , FT  LE =F I NN , FaRM= ’ FORMATTED ’ ) 

READ! 3,  DFNAME 

READ ! 3 , 4 ) S I G , GL , CLEN , AD , VTM , THK, TCP ,  WB , WP , WT . FL , DBB , SPD , PO 
READ (3, 4) DUt , DU2. 0U3, DU4 , DUS, AGAS.ELL, GH,  XD, XB, XL, TS,  X  IN, X INN 
DO  100  1=1,30 
READ(3,S)E!I, 1),E!I.2) 

CLOSE !3) 

DRAG  MODEL  COEFFICIENTS 

WRITE!*,!)’  INSERT  APPROPRIATE  DATA  DISK  IN  B-DRIVE  AND’ 
WRITE!*,!)’  ENTER  FILE  NAME  OF  DRAG  MODEL?  D: -f  i  I  espec .  RFV’ 
READ!*. ! )FMOD 
WRITE!*,!)’  ’ 

OPEN  ( 3 ,  F I  LF=FMI)D ,  F0RI1=  ’  F  ORMATTE  D  ’ ) 

READ(3,9)FXM013 
WRITE!*,  DFXMOD 
WRITE !♦,!)’  ' 

REAP(3„9)JSEG 
DO  110  J  =  !,,ISEG 
DO  DO  1  =  1,4 
REAP (3, 1-)RFV!I.J) 

L.DSE  !3! 


INPUT:  INTERVAL  GROWTH  LIMITS  FOR  RUNGE-KUTTA  VARIABLES 
WRITE<*,n'  ENTER  DIFFERENTIAL  GROWTH  LIMITS:  ALT.  C  ft  3,’ 
WRITE(*,1>'  ASCENT  RATE  C  ft/sec  3  AND  TEMP.  C  deg  R  3' 

READ<*, i5)ALlM,VLIM,TLIM 
WRITE(*,1>’  ' 

OPTIONS:  MAXIMUM  INTEGRATION  TIME  INCREMENT.  INTEGRATION  TIME 

INCREMENT  FOR  VENTING  FROM  DUCTS,  ENABLEMENT  OF  ENERGY 
ABSORPTION  BY  INFLATANT. 

WRITE<»,1)'  TO  CHANGE  DT  t«a«3  OR  DT  Event!  ENTER  NEW  NON-ZERO’ 
WRITE!*, D'  VALUE,  OTHERWISE  ZERO.  ENTER  ZERO  OR  ONE  FOR  AGAS. ’ 
WRITE!*, 1)’  DTM,  DTV,  AGAS  ^  inflatant  absorbs,  0/l=N/Y  3’ 
WRITE!*, U) DTM. DTV. AGAS 
READ!*, 15)GG1,GG2,AGAS 
IF  !QGi.6T.O.)  DTM=GG1 
IF  ! G62 .  GT .  0 .  )  I)TV=662 
WRITE!*,!)’  ’ 

option:  PARTIAl.  FLIGHT  PROFILE  COMPUTATION 
WRITE!*, U’  TO  COMPUTE  PARTIAL  PROFILE,’ 

WRITE!*,  1)’.  ENTER:  TIME  tsec3  «(  ALT.  Cft3’ 

READ  !*,  <b)  TSTOP,  ASTQP 

OUTPUT  TO  PRINTER 
WRITE!5, 18)FDUT 
WRITE!5,1)’  ’ 

WRITE !5, ia)FMOD 
WRITE!5,1)’  ’ 

WRITE!5,1)’  DTM  DTV  AGAS  ALIM  VLIM  TLIM’ 

WRITE !5, 7) DTM, DTV, AGAS, ALIM. VLIM.TLIM 
WRITEO,  1)’  ’ 

PRINTING  OF  PROFILE  HEADINGS 
WRITE ! 5, 1) HEADER 
WRITE!S,1)'  ’ 

INITIALIZATIONS  AND  NON-RECURRING  COMPUTATIONS 
LBTEN=LOGI10.  ) 

IN=INTIXIN> 

LH  =  IN 
LLO=IN 
EGAS==AGAS 
AYR6=AGAS*AYR.; 

AYRGl=t.-AYRG 

AV=AYV* ! 1 . ♦TYV*AYRGi / ( t . -RYV»AYR6l ) ) 

AG--AYRG#TVV/  <  1 .  -RYV*AYRG1 ) 

Dl'B=DBB/60, 

WTX=0. 

SPD=1.6‘('*SP0 

DD4=P0 


n  n 


rf'7I^TVm^F.M'n^KITKIT12-RW■UI'nJRAJl'nJhSnll-^lJrnJXnJUnl-^'^r^tAnMVk^’«V^W'^  V" 

j 

I 


C  IK  MODEL  [  SIMILAR  TO  REF.  15,  PAGE  57  1 
DO  190  1=2,30 
TROP  =  E<r-1,  1) 

IF  <E(I,2) .GE,E<1-1,2) >  GOTO  200 
190  CONTINUE 
200  TIR0=l.a*E(l.2)-5.55 

TIR1=.74*T1r6 

DT I =- . 26#T I RO/ < TROP-E < 1 , 1 ) > 

j 

ALBEDO  MODEL  C  SEE  REF.  15.  FIG.  10  ] 

NOTE;  XL  IS  IN  DEGREES  IN  THIS  MODEL. 

IF  (XL.GE.20.)  GOTO  300  • 

RL=-.0025*XL4-.  15  = 

GOTO  600  ' 

300  IF  (XL.6E.30.)  GOTO  400 

RL=. 1  3 

GOTO  600  ■ 

400  IF  <XL.GE.40.)  GOTO  500  I 

RL=.005*Xu.-,03  I 

GOTO  600  'i 

500  RL=.0O‘/S*XL-.  15  i 

i 

C  CONVERBIONS  TO  RADIANS  | 

600  XD=C<1)*XD  i 

XL=C(ll»XL  i| 

I 

CX=COS<XI..)*£:OS(XD)  j 

SX»S1N<XL)»'3IN<XD)  1 

C  APPROXIMATE  CAP  WEIGHT  ROUTINE  j 

CALL  UPSON (1, SI 6, ADD  i 

GB=CLEN/GL  '  ' 

CALL  NELSON  <U0,6Ei.SlG>  I 

WCsUOx  Att.  *  WOWkTCP«GL»GL  j 

C  RUNGE-KUTTA  ZEROING  ROUTINE  J 

DO  700  J=l,5  i 

DO  700  1  =  1,  a 

Y  < 1 , 0 ) =0 . 

VO  '  (P  (  I,  J >=!.>. 

C  RO'JIINE  TO  SET  INITIAL  VALUES  OF  F\-(..  VARIABLES 

Y  (1,1X1. 

Y(2,  1)-Agp+WT 
Vf3,l)=ELL 

CALL  VlR(}NvEU..BP,DDl,DD~,DC5,LLl,Ll3,DD7,DD9,TK.  n 

BE^  (1  .  •  .  .»3Bia5>*SW 

V0--  ( 1  .  *  .-  )  *  (WB-*  WP+WT)  /'BE 

Y(-'t,l'MR 

V'5.1.'---‘’'i*FE/7R'RG 

W.v,t'  = 

Y  I  1  '  v-O 

Y  (B.  I 


69 


!  ‘ 


1 

< 


f  *  +++  ♦  •♦  •f-»  -*‘4-  +  -*  <  •l••f+•^^  +•♦*'f•♦••^•»•4•+•f  •♦-♦  •♦  •♦■+■♦••♦-+•♦  •♦  +  ^•^■f++•f+•♦•4•+  +  •^•+•f  4- ♦•4-l-4-++4-++-f 


♦  BEGIN  I^UNGE-KUTTA  ITERATION  + 

+  44+4+4  1  44+4-t  4  •♦  14  4-+++ V4-4  +4-+44*f4  +  »>+4-f4-f4»f+‘f+4*4-4  f4  +4-++++4-4*++++4-+++++-f++4*+ 


{ 


I 

>< 


i 


I 


! 

■4 
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BOO  DO  1600  K.J-2,5 

h::»;>K,1~l 

CALL  VIRCIM(Y<3,KK)  ,  DP.  DDl ,  DD3,  DDS.LLI ,  LL3,  DD7.  DD9,  TK.  1 ) 
BE  =  <1.  --.  138ia5<^rR/Y, 6,(00  )«5W 


C  (^OLITIME  (0  ADJUST  MAXIMUM  BALLOON  VOLUME  DUE  TO  DEBALLASTING 
IK  (WT.NE.W1X)  THEN 
X=LOe  <  1 ,  +WB/  <WT+t«(!-  >  > 

CALL  UPSON ( 2, X . S IGX  > 

CALL  UPS0N(3,SIGX,VU) 

Vr==VIJ-*CiL**3 

WTX=WT 

IK  (VT.GT.VIM)  VT=VTM 
EWDK- 

IK  ((kK.EO.  n  THEN 

c  eractional  volume 

UP=Y  <7,KK) /VT 

C  POUriNE  TO  DETERMINE  AREA  OK  CAP  ENCLOSING  THE  GAS  BUBBLE 

r AI.L  MYBLN  <  VB ,  til. .  GB , GN,  DH,  SB,  DM, RM,  HC > 

U«UO 

IF  <GN.I.T,CLEN)  CALL  NELSON (U,  GB,  SIG) 

C  EFFECTIVE  GAS  ENv'ELOPE  FILM  WEIGHT  t  INCLUDUJR  THE  CAP  PORTION  3 

WE=£3D*  (WB-WC)  •t  W0W*U*ABL*TCP*6L*GL 

C;  FFFECTTVE  GAS  ENVELOPE  SURI-ACE  AREA 

SA  •-SB'''ABI.*Gl..-*GL 

C  EFFECTIVE  SURFACE  AREA  ADSORBING  IR  C  SEE  REF.  15,  PAGE  55  3 

SfT=HC 

I F  ( BE .  GE .  0. 0<'626 )  SR=HC  ♦  <  HC-SA )  *L0t5  ( .  00626  /BE  >  /  LGTEN 
END  IF 
P(v=l)mil;iE: 

SWt-i=Y  .-'Y(7,(CK) 

V5CL-A1.  r<i  Y  C'K )  **DET 

CG-GAM-'VSG 

SPEED- y<  1,(0') 

C  PPL  LAUNCH  WIND  SPEED  DURING  THEl^MODVNAMIC  STABILIZATION 
IF  LAUNCH. ED. 0)  SPEED--=SPO 

SPDS«=SPEED#SPEED 
RN“DN*SW»ADS (SPEED) /VS 
BU0y-SW*Y(7,)''K) 

WS-y  ( 2 ,  >  +Y  ( 5, PI-:)  +WB 
FORCE -^BUOY-WS 


t 

■i 

-ji 

! 


'i 

1 


3 


1 

1 

j 

I 

i 

j 
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r  AIV^ODYHAm r  DKt^a  R'OirfTNp  rtIR  RCSINCi  BAUl.CION 

II-  (I  niiNCH.  Lu.  t )  nu-N 

i>n  vr«<.  .1=^'.  jyKb' 

9('o  n>=-i::D  t  t  ,d  >  »  <RN*i*RrV(L',J>  )»(l‘R**RFV<3i,0)  >*(VB*#RFV(4,a> ) 

tNDJI- 

DR:Aki=.  5*RO*C!>*HC»aPEi;i;D*AEiS  (SiREED) 

n  ( 1 , 1  J )  (FORCE-DRAG;  /  (WGhCM»BLIOY) 

C  PRE-LAUNCH  VLI?T:[L'AL  HOIIOM  INHIBITOR 

IF  t  L  AUNCH .  EO  -  0  >  D  (1 ,  i  .J )  “U . 

D(  J.,1  J)  =  Y(1,K»-  ) 

*  «  +  ♦+  +  +  I  4  ^  +  +  +  ++-*-H++Ff++++*"*-+l-f+++-f-l-+4+-F+'++  +  -F+ 

4  BEGIN  HEAT  TI^ANSFER  RATES  + 

C  SEE  REF.  19,  EO.  2a 

E YRB-EGAS*3 -  A 2E -06»  ( Y  < i ,  k'l.: )  / 1 . 8)  **• .  8 1 52 

EYRGI-1 . -EYRG 

ER=EYRii  ( 1 ,  4  TVR*EYRG1  /  ( 1 .  -RYR*EYRB1 ) ) 
mR^^ER 

E I ^EVR6»EYR/ ( 1 . -RYR*EYRG1 ) 

EG-LYRG*TYR/  ( 1 .  •(RYh'*EVIiGl  » 

C  MUSSELT  NUMBER  ROUIINE  t  SEE  THIS  REF-ORT  AND  REF.  18  3 
CAI., L  F'RGR  < BP, .  67,  DM,  BWB,  Y  < 4 , K»0  ,  Y  ( A,  KK) ,  VSB ) 

IF  (OP.LE.  15E+0  7)  6NU=2.-»-.6»6P*a-0.25 
IF  ( lip . (3T  .  1 5E4'07 )  F.NLI= .  I T.»GP** ( 1 .  /3. ) 

C  CQNVfXTIVE  HEAT  TRANSFER  BETWEEN  BAS  AND  BALLOON  ENVELOPE  FILM 
Cl  1  r0 1  *BN  1  *SA*  ( Y  ( 6 ,  kk  >  “Y  < 4,  kk > )  »CG*GNU/DM 

r,!A-PXH:X*COS  (C  C  l )  iKCiH-XB+-T/240.  )  > 

CC  ROUTINE  TO  DETERMINE  EFFECTIVE  SOLAR  ENERGY 

C  SEE  REF.  15,  EGI,  38 

ARMG;=(BP/PO)*  (SORT <  1228.64-376750.  44*OA*OA) -613.  Q*ABS(QA)  > 

C  SEE  REF.  12,  E(?.  48 

TRM1  =  , 5* (EXP 6S*ARMS>4EXP (-.095*ARM5>  > 

FV=GS*CS*TRm 
IF  (QA.LT.O. i  THEN 
FV-0. 

0P=-S(3RT  ( 1 .  -  (RE/  <RE4  Y  <3,  Kk)  >  >**2) 

IF  (OB.LT.(IA)  THEN 

C  OPTICAL  AIR  MASS  ALTITUDE  C  SEE  REF.  12,  EQ.  51  3 

I Z  =SGRT  ( 1 .  -(3A*QA )  •  ( RE-t-  Y  <  3 ,  KK )  )  -RE 


c 


ATMOSPHERIC  PRESSURE  FOR  QPTICAl,  AIR  MASS  ALTITUDE 
CAUL  V I RON ( Z  Z , P AM , DDO , DD2 , DD4 , ULO , LL2 ,  DD6 , DD8 , T  X ,  2 ) 


AM=35. l*PAM/PO 

C  SEE  REF.  12,  EOS.  40  52 

TRM2= . 5*  <  E  XP  < - . 65»AM) +EXP ( - . 095#AM) > 

FV-G5»CS»TRM2*TRM2/TRM1 

END  IF 
END  IF 

CC  END 

C  DIRECT  SOLAR  ENERGY  ABSORPTION 

G!2=C'a2*AU*FV«-HC 

C  ABSORPTION  OF  IR  ENERGY 

(53=C03*AR*BZ  »SR*T  I  **4 

C  NUSSELT  NUMBER  ROUTINE  C  SEE  THIS  REPORT  AND  REF.  IB  3 
IF  (ABS (SPEED) .LT.VL)  THEN 
CALL  PRGR  <  GF . . 67 , DM . SW, Y  < 4 . KK  > , TR . VS ) 

GNU=GN2* (2. +. 6#GP*«0.25> 

ELSE 

GNU= . 37*GN3*RN**0 . 6 
END  IF 

C  CONVECTIVE  HEAT  TRANSFER  BETWEEN  GAS  ENVELOPE  AND  AIR 
Q4=CI34*SA*  <TR-Y  <4,KK> ) *CA*6NU/DM 

C  IR  ENERGY  EMISSION 

Q5*C(J0*ER*BZ*SA*Y  x  4 ,  KK  > 

C  ROUTINE  FOR  DETERMINING  DIRECTIONAL  REFLECTIVITY  FACTOR  C  SEE 
C  REF.  IS,  FIG.  15  3 

ZS=57.29578*ATAN<S0RT( (l.-CA*aA>/(OA*OA> ) > 

FF^'l. 

IF  (ZS.GT.77.2)  FF=. 09375# ZS-5. 4375 

IF  <  (FF.EK.  1.  )  .AND.  (ZS.(3T.25.  ) FF=. Ol53»ZS+. 6169 

C  ABSORPTION  OF  REFLECTED  SOLAR  ENERGY 

G'6^CQ6*AV*(2,*HC>*GS*FF*RL»(l.-SC!RT(Y(3.KK>  /RE/2.  >  )*DA 

CC  07  THROUGH  Pit  ARE  BASED  ON  A  MODEL  PERMITTING  GAS  IMPURITIES  AND 

LL  THUS  ENERGY  ABSORPTION  BY  THE  INFLATANT  C  SEE  REFS.  18  !-  40  3 

IF  (EI.NE.O.)  THEN 

C  RADIATIVE  EXCHANGE  BETWEEN  INFLATANT  AND  ENVELOPE  FILM 

GI7=C07*E  I  *B  Z  <  V  ( 6  ,  F.K )  **4  - Y  ( 4 ,  FK )  *  *  4 )  tSA 

C  ABSORPTION  OF  L'lRECT  SOLAR  ENERGY 

GiB=COE)<>AB*FV*HC 
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•  .  *-  _■>  _*•  V 


c 


ABSORPTION  OF  REFLECTED  SOLAR  ENERGY 
Q9*CQ9*AG»!2.»HC>*6S#FF#RL*<l.-SQRT(Y(3,KK)/RE/2. ) )#G!A 


C  EMISSION  OF  IR  ENERGY  BY  BAS 

C 1 0-CG!  1 0#E6»B7  ttS  A»  Y  ( 6 ,  KK  >  **4 

C  ABSORPTION  OF  IR  ENERGY  BY  GAS 

011=C011»E6*B2»SR*TI*»4 

END  IF 


CC  END 

+  END  HEAT  TRANSFER  RATES  + 

4"4"f -f'i-'f  4-«"f +  -f-f  'f 

D(4,Ka)= (01+Q2+Q3+e4-D5+Q6+e7> /CF/WE 

;;  ROUTINE  FOR  EV-13  APEX  GAB  VALVE  OPERATION 
II-  <  <VV.6T.O. .AND.  (KK.EO.  D)  THEN 
IF  (P6.QE.0.93A)  VV==. 559-PG/8 
IF  <PG.LT. 0.936)  VV=. 72222»S0RT ( < 1 . 072-PG) *P5) 

VVt=6.  I27»VV*S0RT(3W/SWS-1.) 

END  IF 

IFIVD.LE.O.)  IVENT=0 
IF  nVENT.EO.O)  THEN 

C  EQUATIONS  FOR  FRACTIONAL.  VOLUME  LESS  THAN  1 

D<S,H,I)=-VV#SWG 

D<6.Ka)  =  (Q0+Q9~G!l-07-0t0+0U-SW*Y<7,KK)  »'D<3,KJ)  )  /Y(5.KK>  /(CV+R6) 

D  ( 7 ,  K  J )  =  Y  ( 7 .  tCK )  *  ( D  ( 5.,  I'vO )  /  Y  ( ir. ,  Kl : )  +D  ( 6 ,  K  J )  /  V  ( 6 ,  KL  )  t-D  ( 3 .  K  J )  /RA/  TR ) 
DO,K3)--SW»H)<3,KJ) 

ELSE 

t:  EQUATIONS  FOR  FRACTIONA  .  VOLUME  FQUAl.  TO  ) 

I'TaDTV 

VD=A1)*SQRT  <2.  *Gi(ABS  <  Y  (8,  Lf;)  -PE)  /SWG) 

D(;5!,KJ)=-<VV+VD'*SW6 

D(6,Kv1)  =  (Q0+09-Bl-07"Q10+G'l  1-PE*D(5.L-J)  /SWG)  /Y<5,LO  )  /CV 
D(7.);J)=i;i. 

D  ( 8 ,  l; J )  =  Y  <  s ,  lO.. ) «  ( n  <  r, .  t.  J )  / v  <  5 ,  )  *1)  <  6 ,  ko  )  / v  <  6 .  kk  > ) 

END  IP 

DO  )  1=1.8 

op=A().:j>-ti  <D<t  i,i.j)-P(i.J>*0''i  i.n.) ) 

Y(|.;isl.,]).-Y(KI.).K)rI)IIK?l.l 
1 HOO  Q  ( h  I ,  L  J  >  kO  a;  I , » ■). )  •»  3.  »-QQ~C  <  L  J )  *D  ( K I ,  L  J ) 

CONTI  MUF 

(  *  t  ■.  »  t  I  ♦+4+  +  4  «  +4 +4+ 44.4 +4+  + 4.,.,  4  4.,  +4.,  4.44.,  +44  ,.444444  +  444+  4  +  444,  444 

^  END  RIJNPE-I  l.'TTA  Mt'vCTiDN  + 

»  <  <  +  •  4  ♦  4  44  +44  4  +1  ...  -4  4  )  4*44  4+  +  -!  .r-4  *  +4444  4  <  4444444  4444++*  +  ++  +  +  +'.4+  +  +  +  +  +  4  +  +444 
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ROUTINE  TO  ADJUST  THEORETICAL  HAXIMUM  BALLOON  VOLUME  It  WARRANTED 
ON  THE  BASIS  OE  THE  ACCURACIES  OF  ACTUAL  CONSTRUCTION  P'RACTICES  I 
ANOTHER  OPTION  WOULD  BE  TO  SET  VT=Y<7,S),  IVENT=1  AND  TO  CONTINUE. 
IF  <Y(7,5).6T.VT)  THEN 
VT=YC>,  1) 

I VENT- i 
GOTO  900 
END  IF 

C  ROUTINE  TO  DETERMINE  IF  ANY  SELECTED  VARIABLE  GROWTH  RATE  IS 

C  EXCESSIVE  r  NOT  APPLICABLE  WHEN  BALLOON  DUCTS  ARE  VENTING  BAS  3 

IF  < I VENT. EG. 0)  THEN 

IF  (A'-i5(V(t,5)-Y<l,l)).GT.VHM>  GOTO  1700 
IF  (ABS(Y(3;,H)-Y(3,1').GT.AL1M)  GOTO  1700 
IF  (ABS(V(fl.S)-Y(4,l)>.GT.TLIM)  GOTO  1700 
IF  (ABS(Y(6,5)-Y<6, 1))  .LE.TLIM)  GOTO  1750 

C  ROUTINE  TO  ADJUST  INTEGRATION  INCREMENT  FOR  EXCESSIVE  GROWTH 

1700  DT=DT/2. 

IF  DT.LT.O.S)  THEN 

WRIiei*,!)’  BAD  EXIT  I  DT  <  O.S  I' 

STOP 
END  IF 
GOTO  GOO 

C  ROUTINE  TO  ADJUST  INTEGRATION  INCREMENT  TO  CONTROu  ALTITUDE 

r  DURING  AN  ITERATION  INTERVAL 

17S0  DTX=DTM 

DO  1800  1=1.7 

IF  (Y<l,5)»bTX.LT.AHM)  GOTO  1900 
1800  DT>:*UTX-2. 

1900  CONTINUE 

ENDIt' 

C  PRE -LAUNCH  STABILIZATION  CLOCK  TIMER 
IP  r^ALINCH.EO.O)  THEN 
TT=TT-*DT 
Dr=^DT>: 

IF  rn.GE.TS)  THEN 
I.AUI'ICH=  1 
DT  =  1. 

CD=('. 

ROTO  0200 
END  IF 
GOTO  2.;.0i;> 

FMDIF 

C  FLAFSED  FLIGHT  TIMER 

T -T+DT 


C,  ROUTINE  TO  GENERATE  OUTPUT 

2200  VLT=Y{-;j;.IO/Vr 

i;AL.L  M VBI-N  < VB,  CjL ,  gd  ,  gn ,  dm ,  sb.  dm ,  rh,  hc > 

CALL  U1R0N(Y(~,LK)  .DP.DDl.DDS.DDS.LLl.LL'i.DDr.DD-^.Tk,  1) 
l•t?~Y(^.b!/(ti^^DM) 

RN»-T)M»GW*ABS(Y(l,5)  )/VS 
X-60. *Y  <  I  .S) 

WR1TE<5,  l<b)T,Y!3,S)  ,X,RM*CD,PR,RN,VD, 

r  PARTIAL  FLIGHT  PROFILE  TERNINATiON  CHECk 

1F<(T.OE.TSTOP).OR.(Y(1<,5).GF.ASTOP>)  GOTO  2500 

C  SNTEGIWXTKIN  INCREMENT  CONTROL  DURING  flRST  10  SECONDS  OF  FLIGHT 
IF  (T.0T.1O.)  1)T<=:I)TX 

C  ROUTINE  TO  RE-INITIALI/E  RLINGE-KUTTA  ITERATION 
2300  DO  2^00  la 1,8 
Y(  I,  l)«Y(I,ti) 

24*Li  0(1,  l)siO(I,5) 

BO 10  BOO 


250«i  CLOSE  (5) 
STOP 

FNO 


i 

i 

■1 


ii'  r 


> 


N 

I 

f* 

s 

M 

V 


I 


i 


) 

1 
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BLIEiROUT  INF:  Pf<(3R  l  XO,  X  t ,  X2,  X3,  X4,  XS,  X6) 


»  » 

*  rOMPUTfLS  THK  PRODUCT  OF  THE  GRFBHOF  NO.  AND  1HE  PRANDTL  NO.  * 

*  * 

*»**#■«*»»*«■»*■*  im  »*■)(■»»*<  ***#*«**»»*##»*»*»)(  »***«#****it  **»»**»«■*#»»#*»»•«■## 


XO=:v.;.  174*X  1  »X2*  <  ( X2xX?./Xi)  »-*2>  -xABS  ( X4-X5)  /X5 
END 


UUDIVnUT  INF  NELSON  ( T ,  R,  X ) 


»  « 

♦  CONPUTEf)  FRACTIONAL  CAP ‘AREA  t  IN  FULL  OR  IN  PART  3  FOR  A  GIVEN  * 

*  SIGMA  AND  LENGTH  C  AS  A  FRACTION  OF  THE  ACTUAL  BDRELENGTH  ]  AS  ■» 

»■  MEASURED  FROM  IHE  THEORETICAL  APEX  POSITION  * 

»  « 


>»  Hit*  »«««««•*«»«•  »««««  »««(*«<(»-»««*««* 

DIMENSION  W(7.4) 

DATA  W/-. 190090090. -.9529(395, .019023809, .501211094, 

«  1 . 600860007.  1.11871  517:;.,  -  -  24:1042045, .  1 96645137,  4.  6444  2895, 

«  ■  2404 1  1 25 .  0. 472870627 .7.  20250734 1 .  4 . 25933462 ,  .  7859002 , 

»  3.  38  32fc12'’a2.  -6. 60.O2O4O:.'5,  4. 076000773,5.  56431556.-10.  053917, 

«5. 4<>4<  .9 12.  ■■ .  859152769,  -V.  3989509. 2.  9527 1 8 136,  - 1 . 645782348, 

■<  -2 . 706.3  ?  1B58 , 4.  499006036,-2.  28485046, .  3252352:r.5/ 


Y^O. 

XX.=-|.  R 
D>'  2<'0  ,.-1,4 
X  -  W  <  t , .) ) 

Oil  111.)  1-2,7 

lOO  2X-ZX‘WiI, J! »X*« (I -1 1 
2<m)  Y -y‘7X**X**  (J-l) 

Y-1.  Y 

END 


VG 


non 


SUBROUTINE  MYBLN  ( VBjBL,  GB.BN,  DH,SB,  DM,  RM,Hn) 


»  » 

*  COMPUTES  GEOMETRICAL  CHARACTERISTICS  OF  THE  SIMPLIFIED  PARACHUTE  * 

*  SHAPE  BALLOON  * 

*  * 


**#****»#**«*#****#*-»HHHI-**##fHHH(-************»*#**»***#***»#*#**«  «■*****#» 

COMPUTES  DIMENSIONLESS  HEIGHT  COORDINATE  (ZB)  OF  THE  GENERATOR 
SHAPE  t  AT  LINE  OF  INTERSECTION  WITH  THE  BABE  CONE  D  AS  A  FUNCTION 
OF  FRACTIQNAl.  VOLUME  <VB) 

CALL  BRAUNO.VB,  ZB) 

C  COMPUTES  HALF-ANGLE  (TH)  AT  NADIR  AS  A  FUNCTION  OF  ZB 
CALL  BRAUN (1, ZB, TH) 

C  COMPUTES  DIMENSIONLESS  HORIZONTAL  RADIUS  COORDINATE  (RED  OF  THE 
C  GENERATOR  SHAPE  C  AT  HEIGHT  ZB  I  AS  A  FUNCTION  OF  ZB 
CALL  BRAUN  (2,  ZB,  RB) 

C  COMPUTES  DIMENSIONLESS  GORELENGTH  COORDINATE  (00)  OF  THE  GENERATOR 
i:  SHAPE  C  AT  HEIGHT  ZB  .1  AS  A  FUNCTION  OF  ZB 

CALL  BRAUN  (3,  ZB,  00) 

C  CONVERTS  TH  FROM  DEGREES  10  RADIANS 
TH«0.  Ol74'j3292S2*TH 

C  computes  GB,  THE  RATIO  OF  THE  GENERATOR  SHAPE  GORELENGIH  TO  THE 
C  ACTUAL  BALLOON  GDRELEMGTH 

GB'I . /  (1. -OO+RD/SIMiTH) ) 

l:  aiMPtllH'-:  THE  GENERATOR  SHARP  GDRfi LENGTH 

BN-GL*GB 

C  CCJMPUTES  THE  MAXIMUM  HORIZONTAL  DIAMETER  OF  THE  GENERATOR  SHAPE 

i:  (  IHF  niAMEIER  OF  IMF  IDEALIZED,  PARTIALLY  FULL  BALLOON  i 

LV1=  ,  <5031_’*GN 

r,  CONVEPIS  ftTE  DIAMCIER  TO  THE  RADIUS 

RMrr|)M/;‘. 

C  COMPUTES  I  MF  VOL  I  IMF  01"  IMF  TANGENT  CONE 

vc='5.  /tamcth) 

C.  COMPUTLS  UIF.  HEIGHT  OF  THF  GENERATOR  SHARE 

ZM"' .  (ii ! '.TT’1‘11  f  >N 

Ci.M-ir-lHI  M-IF  .'GUI  HI  IMF  MF-LIUM  FU.LED  JME 
DH-.' M  ;B*(jN 

I.  I  lMl.N'.  l  >.■  li  ffc-RTNCF  RFIWEFN  THE  ACTUAL  BAl  LOON 

I  OTIRFl  I  NITU!  AN!;  TH!:  IT'AMEOUS  GENLiA'  p  SHAPE  OORCLENGTH 

0--1.  b 
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COMPUTES  FRACTIONAL  SURFACE  AREA  OF  ACTUAL  BALLOON  BELOW  THE 
LOCUS  OF  POINTS  CORRESPONDING  TO  DIMENSIONLESS  LOCATION  (0! 
CALL  BRAUN <4, 0,00! 

COMPUTES  FRACTIONAL  SURFACE  AREA  OF  THE  ACTUAL  BALLOON  SHELL 
ASSUMED  TO  APPROXIMATE  THE  AREA  OF  THE  ENVELOPE  SURROUNDING 
THE  GAS  BUBBLE 
SB=1.-00 


COMPUTES  AREA  OF  THE  MAXIMUM  HORIZONTAL  CROSSSECTION  OF  THF 
BALLOON 

HC=3. 1415927#RM»»2 
END 


SUBROUTINE  UPaQM(J.X,Y) 


-)»  # 

»  ROUTINE  TO  COMPUTE  CERTAIN  CHARACTERISTICS  OF  THE  UP!';ON  NATURAL *  * 

»  SHAPE!  BASED  ON  TABULAR  DATA  ORIGTNAI.LY  PREPARED  BY  0 .  H.  I5MAU-EY  * 

*  * 

»  FOR  J-1,  OUTPUTS  NON  DIMENSIONAL  SURFACE  AREA  ♦ 

»  A/S*»2  =  FCNiSIGMA),  R**2  =  .999799711  * 

*  « 

*  FOR  J-2  OUTPUTS  SIGMA  * 

«  SIGMA  ~  FC:N<G/(.  ).  R**2  =  .  99999999G  * 

«  * 

*  FOR  J-3,  OUTPUTS  NON  DIMENSIONAL  VOLUME  * 

*  V/S«#3  =  PCN<SIGMA>,  R*«2  =  .999999324  * 

*  * 


*****#**-***f#****-*#***»****»-**»***#*##***»'!<*#*«****«-*#****-»************tU<- 

DIMENSIUM  UF--<S),3) 

DATA  UP/ 1 . 23S90478fi,  .388987425-1 ,  -  .048794261 .  374649865, 

*  .  1451693335,  .  1961  1  38n:>,-.  1  1 701  7 1 7t.9,  6.  ,0.  37546901  4941 , 

*-.061  4H'l4-’:9632. .  f'l  793  '454(.>922,  -  1 . 839/2293864F-03, 

*2.  ■.•"•■9;.’.r.;  i7K;3E-04,  -1 .4.^41  7B71814E-05,  6.  .  .  1  260550EL 
*."<.1  >  -4  In'/flSf  ,“W.fj'.Vii/'>9/4*:.4I:-03.-  . i.>9i;)74n353, 

*.  OH'iuU  Um;  J33.  -6.  968t.>22fci::l  7E-‘.>3,  -  .  i.>!  i.'l  10745  7  963, 6.  / 

I  Ul-'  1  ,T,' 

I  1  U-^lNmipcH,  J)  ) 
oil  ! Hi •  1  !  .  1  J  ! 
t ■  V+UP ' 1  *  1 , J ' *' »* 1 


1.1'. 


.tw  v'n>w‘it.vHvu  jbif'fviM  ntb  /v)i' ru(f/W««»\rtt Wlv  vW'W‘nk’U;Vl'fVTrwwn*H’(/'W WWVH^W^Trjl 


SUBROUTINE  BRAUN(J,X,Y) 

***##»')ne*********»**#-i<-****#*****»***#in»-***»#*-K*»**»#*#***#)H(-»****innn(  >♦* 
■*  * 

*  MATHEMATICAL  MODELS  OF  FUNCTIONS  USED  IN  THE  SUBROUTINE  MYBLN  * 

*  * 

DIMENSION  BR(I0,15) 

DATA  BB/SA.  151674, 3:6.  9<.94I5, -8A3. 384153, 4340.  &007B. -23751 . 1652. 
*60663. 6084,-55452. 3401 , 0. .  2. ,  7. ,0. 00012929,  1 . 491 118, 

*0 .  89926756885  ,-12.491 6446722 . 15 . 759250787, 10.1738213142, 

*-26. 202606089,0..  2.  ,7. , .000128. 1 . 791098667,  0. 901 1 1501S, 

*-12. 302923389,  23.  701897323, -3. 8057025141,-8. 2295252039, 

*0. ,2. ,7. ,-0.00037994265363.0.019166414526.  1 . 07S93200657, 

*0. 10690678257,  1 .595f.4475729,-5. 19740094837, 2. 5981953799, 0.  ,  2. , 

*7.  , 1 . 0053899018856, -.2463831 77,0. ,0. ,0.  ,  0. , . 999695,8. ,0. , 

*13.  47532639,-13.453150437.0.  .0.  ,o!  ,0.  ,ol. ,  .998704865,  1.  ,2.  , 

*5. 70116145, -5. 668903443,0.  ,0. , 0. , 0. ,ol , .  996376. 1 . ,  2.  , 

*3. 17395699766,-3.  13250803076. 0. , 0. , 0. , 0. , 0. ,. 99216,  l.,2., 

*15. 735574 1 5, -67 . 659 1 64535, 90.679194639, -38.  706776393,0.  ,  0.  ,  0. , 
*.963, 2. , 4. , .  92472377876, 1 . 03075490637,-13. 8580229434,0. , 0. . 0. , 

*0.  .  .83035916,,  8.  ,0.  .  .  31091805108. .  27740682536,  -2.63309861357, 

*7. 406)32484  ;  -1 i  ,  064485699,8.  48105265334,  -2,69143362928. 

*. 498,?. .7. , .  340858423189.-.  18492521 608, 0, ,  0. , O. ,0.  ,0. , .  27218.  1 .  . 
*2.  ,  .  36290-' 275675,  - .  437916794366,  1 . 02361964267,-1 . 44464774489. 

♦0.  ,0.  .,0.  ,  .0508769,2.  ,4.  ,5.3392134074.-28.  1122B36054, 

*37 . 0 1 4 505 1 1 27 , 0 .  ,  0 .  ,  0 . , 0 . , .  00324869 , 8. , 0 . , .  383 1 79458027 , 

*-2. 5872040881 1,0.  .0.  ,0. .0.  ,0. ,  . 00015, !.,?./ 


(H  (J.l.  I.b)  8010  2<.Hi 
DO  100  J-5,15 

IF  <X.GT.9R(8. J) )  GOTO  200 

CONTINUE 

LL“lNT<E<Ff9,J)  > 

riti=iNT(eRno.  j) ) 

IF  'LL.E0.9)  THEN 

Y=  i'-BF.  (2.  J>-S0RT<BR(2,  J>*BR(2..J)  -4.  *BR(3,  J>*<BR(  1  ,..l) '  X) )  ) 
Y«V  /  ('2.*EiH(3,J)  1 
RETURN 
END  IF 
Y^BRO  .  J) 

DC)  3iin  I  =2, MM 
r=:-^friiR(  !  ,,I)*X*#  (l-l  : 


VMi  aval  9>/£M^tX-^ytliUm^\'0*XUAn\  r  tA»' AX<>^ib'i\»C^rN^AjC‘^VW!A»n9UQS>7Q>mV9!Rrni/?O0V)K1 


SUBROUTINE  VIRON< XA, PQ,EF,Ee, PX, JX, JV, RR, DY, TK,NN) 


*  « 

*  ATMOSPHERIC  MODEL  BASED  ON  1962  U.S.  STANDARD  ATMOSPHERE;  * 

»  WHEN  NN  =  2,  COMPUTES  ONLY  PRESSURE  FOR  AIR  MASS  ALTITUDE  * 

*  # 


COMMON  CAjDTI ,E(30,2),PE,P0,R0,SW,TI,TIRl,TIR0,TR,TRaP,VS 

DATA  (:C/34. 163195/ 

X=.3048037*XA 
X=6356.766»X/ (6356766. +X) 

IF  ((X.LT.E(JX,l)I.OR. (X.QE.EF))  THEN 
IF  (X.LT.E<JX,U  >  THEN 
JV=1 
PX--=PO 
GOTO  200 
END  IF 

100  IF  (ABS(DY) -GE.O.Ol)  PX=PX* ( (E (JX, 2) /EGf ** (CC/RR) > 

IF  (ABS(DY) .LT.O.Ol)  PX=PX*EXP (-CC»DX/E ( JX, 2) ) 

200  JX=JV 

JV=JV-»  1 
EP-E<.JV,1> 

ES'-E<JV,2) 

DY=ER-E(,1X,2) 
nX:=EF-E<JX  .1) 

RR=I)Y/DX 

IP  (X.BE.EF)  GOTO  100 
FNIUF 

DX=X-E (JX, 1) 

TK«E(JX,2i+RR*DX 

IF  (ABS5<DY).BE.0.0.)  P(3=PX#  <  (E  ( JX,  2> /TK)  *#  (CC/RR)  ) 

IF  (ABS<DY).  LT.O.Ol.)  FT3=PX#EXP  (-CC*DX/E  ( JX,  2)  ) 

IF  (NN.ED. n  THEN 

r:  TEMPERA  rURE,  RANiaNE 

Tfx'=  1 , 8*T); 

C  ATMOSPHERIC  PRESSURE,  L.BS/FT**2 

PE-2.0BS50*P(3 

C  SPECIFIC  WEIGHT  OF  AIR.  LBS/FT*+3 

GW:-.02t74fi4*PD/TK 

C  AIR  DEWSIIY,  SLU(3S/FT#*3 

Rr)^-SW/3:'.  1741 

C  VISCOHi  rV  OF  AIR,  LBS.'FT/SPX 

VS--7.  :r.02'3527h:-n7»  <TR**I .  5)  f  (TR+I98.  72) 


w<i*mrru>-vn3T’mr*T«jomjrrMir.if>nirtrpnffKftH'OTtrwjTr'njTTTir'^'.”J<-rjl>?\Jlfkwrvji-nyJ^ 


THERMAL  CONDUCTIVITY  OF  AIR,  LBS /SEC /DEGREE  RANKINE 
CA==3.  30S17E-04»  <TK#*1 .5t  /  <TK'*-24S.  4*  <10**  (-12.  /TK)  )  ) 

RADIATION  TEMPERATURE  OF  AIR,  RANLINE 
TI=TIR1 

iP  (X.LT.TROP)  TI=TIRi:hDTH*<X-E(1,  1)  > 

END  IF 


s'.=  .s  ess;:.-«:ssi5esiB:ss:sv.LJSi::e:i:istssS'Vts>^e£B4ae3s:nc=3:%3:ts?s:3SS:neatafes;es;e«9aB>sci9S:»:»  BS9ea:ese«;aeBsns3<9C*«ait9  9 

SUBROUTINE  INFORM 

*  «• 
*  ACCESSES  PROGRAM  NOTES  t  AUTHOR'S  STANDARD  ROUTINE  1  •* 


)i  •  T  I 

IS 

M 


DIMENS,.''M  -  INE(25) 

CIIARACTER»72  LINE 
CHARACTER*  15  FIN|vl,FNAME 
1  FORMAT (IT) 

?  FORMAT (A) 

T.  FORMAT  (4X,  A) 

4-  FORMAT (A, \) 

F I NN=' ELITE. FAX' 

L-0 

OFCN( 1, FILE «F INN, FORMa' FORMATTED’ ) 

READ(1,2)FNANE 
READ (1 .  n  I 

IF  (J.NE.67B)  GO  TO  9 

C  l:  NUMBER  OF  LINES  OF  TEXT  UNIT  TO  PE  DISPLAYED 

i,  READ  (1.1)1 

IF  (I.EO.O)  GO  TO  9 
N=0 

IF  fI.LE.20)  N=(24-I)/2 
CALL  CLEAR (25) 

DO  7  k=--i,i 

7  READ ( 1,2) LINE (K) 

DO  8  Kal.I 

B  WRITE!*, T)LINE(K.> 

IF  (N.6T.O)  CALL  CLEAR (N» 

PAUSE 

IF  (L.ECLl  )  GOTO  6 

NRITE(»,4>’  DISPLAY  PROGRAM  NOTES  ?  t  0/1  =  N/Y  3  ’ 

READ(*, 1)L 
IF  (L.EO, J )  GOTO  6 
9  CLOSE ( 1 > 

CALL  CLEAR (25) 

END 
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HUB5;!3!,iTlME  CLEAR <J) 


#»-it#»*##**#'**-!HH<-*lHHf*****#**»**»**»*»*»'lt*#*#*#*#*#»«»-X***»*##**#***1HHHI-» 


»  « 

*  CALLED  ONLY  BY  SUBROUTINE  INFORH  TO  CLEAR  SCREEN  OF  PREVIOUS 

»  NOTE  C  AUTHOR’S  STANDARD  ROUTINE  j  ■» 

*  * 


**-##*-S'#»Hnnnnnnni-********-»*-)m-»*»'t(«**#)HHH(***#tH(-**'»'im-»**#**#***«"»»*****#*-: 

1  FORMAT <A> 

DO  2  1  =  1,  J 

2  WRITE (#,1)’  ’ 

END 

===!==,==r==.===-.=r===,======:======  end  of  program  =================-.====-.-,^====.= 
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••V  -fW  *  -V  •  'I 


•  A  -.1  V^  H.'  M.1  ■•<.•.  •<  X V  1  MAX  Axnxx 


PROSRAM  CDMODEL 

#»)nnni«-)i-####»**#»innnnnr  ■»#*■«•##****■!*■»*•#*  «■*»#*##*»«■»*•■**)*•■»••»•**«**«■  **#**##*♦■» 


♦  « 

*  PROBRAM:  CDMOUEL  30  APRIL  * 

■«■  ■» 

*  THE  EXECUTABLE  VERSION  OF  THIS  PROGRAM  IS  DESIGNATED  AS  'CD’  * 

*  WHICH  HAS  BEEN  COMPILED  LINDER  MICROSOFT  FORTRAN  77  TO  BE  RUN  * 

ON  AN  IBM  PC.  * 

*■  « 

*  THIS  PROGRAM  FORMATS  DRAG  COEFFICIENT  MODELS  FOR  USE  WITH  THE  * 

*  EXECUTABLE  PROGRAM  'FLITE’ .  * 

«•  * 

*  THIS  PROGRAM  WAS  DEVELOPED  AT  AND  FOR  THE  AIR  l-ORCE  GEOPHYSICS  ♦ 

»  LABORATORY  AS  PART  OP  IN-HOUSE  WORK  UNIT  NO.  7i59llU.  * 

*  * 


*#-»»»**#**»»*##»*»*-»*-IHHHH<-***'lf#»»»**##*)HHHH******'tH>#*»'IHHr##*****-*'IHf»*'*** 

DIMENSION  PFV(4,20) 

CHARACTER* 15  FMOD 

2  FORMAT (A,'  ’,\) 

3  FORMAT (13) 

n  FORMAT (A) 

13  FORMAT (E15. 8) 

19  FORMAT  < 1 X . ’  CQEP .  FOR  TERM  ’ , 13. ’ ; ’ , \ ) 

20  FORMAT  (IX.’  REYNOLDS  NUMBER  EXP0(4EMT  ’,\) 

21  FORMAT (IX,'  FROUDE  NUMBER  EXPONENT  ’,\) 

22  FORMAT (IX,’  FRACTIONAL  VOLUME  EXPONENT 

CALL  INPORM 

!•;«.'  WRirE(*,2>'  ENTER  NUMBER  OF  TERMS  IN  MODEL:  t  .LE.20  3’ 

READ(*.  3)JSE(5 

IF  (JSEG.GT.20.0R. JSEG.Ll  . 1)  GOTO  100 
WRITE!*,  ID’  ’ 

DO  no  D  =  1,3SEG 
WRITE!*,  19),? 

READ!*, 13)RFV( 1,  J> 

WRITE!*,  ID’  ’ 

WIT  I TE  (*,20) 

READ!*,  1.3)RFV(:;,0) 

WRITE (*,21) 

READf*, 13)RFV \3, J) 

WRITE (*,22) 

READ!*, 13)RFV<4. J) 

Ui;i  WRITE!*, ID'  ’ 

WRITE!*,  ID’  INSERT  PROPER  DIS).  IN  Es-URIVE,  AND' 

WRITE!*.:;’  ENIER  MEW  MODEL  NAME:  B:  4  i  lef.per. .  RFV' 

READ!*,  IDFMDD 

OPEN  ( .  I- 1 LE-?'MOD,  STATUS -■  ’  NEW ' ,  FORM-  ’  FORMATTED ’  ) 

Wk!  IE  <3,  ID  FMOD 
WKr?E.::.,3)JSf;G 
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to 


:/j  Vi  v^,' v'j  vvir.; v;  ■»'/ vi%'v  wir.nrt’ v'jvj' ; v;v'/rvnrtvv'TK'’^vvrb» 

J 

J 


DO  120  J=1,JSEG 
DO  ir.'O  1  =  1.4 

120  WRnt(3,13)RF0<I.J> 

CLOSE  (3) 

END  -  ' 

SUBROUTINE  INFORM  = 

C  3  SEP  1984  • 

DIMENSION  LINE (25) 

C:HARfiCTER-!^72  L-IME  : 

CHARACTER  #15  FINN, FNAME  j 

FORMAT  (1.5) 

FORMAT (A)  ^ 

FORMAT (4X, A)  i 

FORMAT (A, \) 

FINN^'CD.FAX'  i 

L=0 

OPEN  < 1 , F I LE=F I NN . FORM= ’ FORMATTED  ’  >  i 

READ <1.2) FNAME  \ 

READ  <1,1)1  i 

IF  <1.NE.‘j78)  60  ro  9  j 

1:  NUMBER  OF  LINES  OF  TEXT  IJNiT  TO  BE  DISPLAYED  • 

Is  READ  <1,1)1  ] 

IF  (J.ED.O)  60  TO  9  ] 

N=0  1 

<F  (J.LE.20)  N=<24-I)/2  \ 

CALI.  CLEAR  <2S> 

DO  7  K=l,l  i 

7  READ'1,2)L.INE<(::> 

DO  8  F“1 . I  I 

a  WRITE  U.,3)HNE<l->  ] 

IF  IN.  61.0)  CAM.  niF'ARINj  ■ 

PAUSE  ; 

IF  <L.Et.I.  1  )  GOTO  6  j 

WRITE<*,4)’  DISPLAY  PROBRAM  NOTES  ?  I  0/1  =  N/Y  3  ’  I 

REAt)<»,l)L  i 

IF  <L.6Q. 1)  GOTO  A  • 

9  CLOSE  < 1 )  i 

CALI.  CLEAR  <  25)  1 

END  I 

SUBROUTINE  CLEAR <J)  | 

1  FURNA'f(A)  \ 

DO  2  1  =  1,  J  I 

I'  WRITE!#,  J)’  ’  5 

END  ' 


j 

■] 

i 

I 

j 


.v'\ 


1  V  1 » •?  *  X  « 


xV  '■*  < . 


^  ^  ^ LP.J v.iy,^v 


c 

c 


1 


4 

5 

6 
7 
B 

9 

10 
11 
12 


C 


SO 


PROGRAM  DfttftraRN 

ANdV-Jr  ^F[ue?EXE’?°°''''''  'OCD.eXE' 

DIMENSION  X  (2(1)  , E  (9, 2) ,  V  (30, 2) ,  Z  1 100,  4) 

CNARAC rER*64  FINN,FMAT 


«--a  A,-  '^r.  /^'-7  98. 7, 288.  15, 216. 65 

DATA  AL-N,AC4,bl/i*0./,CC/34. 163195/, OD, JIFF/S,  1/,H1/1013 

Cc^ilespec.FLV 

FORMAT (13)  ^ 

FORMAT (IX.FS.O, 1 X , F8. 0, 1 X, E15. 8) 

FORMAT (1X,FB.0, 1X,F8,0) 

FORMAT  (•'  NEW  FL  IGHT  DATA  FILE  tN/Y  =  0/13  ^  '  S) 

NEW  (:, AS/FILM  COEFFICIENT  FILE  CN/Y  =  0/ 1  D  '> 

NEW  HF.AT  XFER  COEFFICIENT  FILE  CN/Y  =  0/1]  '. 
ENTER  NAME  OR  OUTPUT  FILE;  B: f i lespec- OAF 
ENTER  NAME  OF  OUTPUT  FILE;  D: T i lespec. CMV 


FORMAT ( ’ 
FORMAT ( ’ 
FORMAT ( ' 
FORMAT C 
FORMAT  (A) 
FORMAT  <  ’ 


IS  INPUT  FOR  PROGRAM  “FROUDri"  CN/Y  =  0/1]  7 


,216.65, 

.25/ 

S\) 


!.V> 

%X) 


CALL  INFORM 


CRiEATES  D;  f  i  lespec .  FLY 
WRITE!#, 6) 

F^EAD(#,3>  I  AM 
IF  (lAM.EO. 1)  THEN 
WRITE  I#,!) 

READ(v,  IDFINN 
FMAT=' (Eia.G) ’ 

WRITE (*, 12) 

READ(#,3) JFP 

IF  (  (ilFD.  NE.  0) .  OR,  (.TFD.NE.D)  GOTO  50 
write:  (*,2' •'  SIGMA' 

READ(#,FMAT)X<1) 

WRITE!#, 2)'  GORELENGTH  CFT3’ 

READ(#,FMAT>  X (2) 

WR !  I  E  ( * ,  2 )  ’  CAI  •  LENGTH  t  FT  3  ’ 

READ{»i.FMAT)X<3) 

WRITE!*, 2)'  TOTAL  DUCT  AREA  ISO  FT 3’ 
READ<*,FMAT<X(4) 

WRl1E(#,2)'  MAXIMUM  BALLOON  VOLUME  rcu  FT3' 
READ(*,FMAT) X (5) 

WRITE!#, 2)'  FILM  THICKNESS  Cmils]' 

READ (*,FMAT) X (6) 

WRITE!*, 2)'  CAP  THICKNESS  Cmils]' 

REAI>!#,FMAT)  X  (7) 

WRITE!*,';)’  BALLOON  WEIGHT  CLBS]’ 

REAI)!#,FMAT)  X  (0) 

WRITE!#, 2)  ••  IRREDUCIBLE  PAYLOAD  WEIGHT  CLBS]' 


C'-' 
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<0^  wi'H/’  »wr^rviWT‘^-W'WT-*u^  "(/tjmi iviWT^ffS^wl'Rn 


Ri::ni)(*,rHAT)x  c?) 

wru  1 U  ( * ,  2 )  ’  BALLAST  WE  I  Gt  IT  C  LDS  3 ' 

REAn(*,FNAT)X (10) 

WRrTE(K,2)'  FRACTIONAL  FREE  LIFT’ 

REAIi()(,FMAT)X(U) 

WRnE(-«,2>’  BALLAST  RCiUR  RATE  LLBS/MINI’ 

READ(*,rMAT)X (12) 

WRlICi*,!’)-  GROUND  WIND  SPEED  l  lnoi-.s.r 
RE.M)(iiJ-MAT.)X  (t:.) 

WRITE  (K,  2)'  LAUNCH  SITE  ATH.  PRESSURE  fitibs]’ 

REAI)(  ).\FHAT>X  (14) 

X  (  ITj)-^!. 

X  (  I  b  >  - 1  . 

X(t7)-M. 

X(  18)-l.. 

Xi't'A  -1. 

rr  (JFD.ED.A)  THEN 

WRIIIKK.^)-  ESTIMATED  DRAG  COEFFICIENT  CO. 83’ 
READO.FI'IAT)  X(1S) 

WRnF;(*,L’>  ’  MAX.  REL.  ERR.  IN  ALTITUDE  CLOSURE  CO.  OIT’ 
READ(#,ENAf)X(l6) 

WRnF(*,2)'  DRAG  COEFFICIENI  ADJUSTMENT  CCD/AT’ 
REAIJ(*.FMfVI)X(18) 

WR1TE(*,2)’  ENABLE  GAS  ABSORPTION  CN/Y  =  0./1.3’ 
READ(*,FMAT)X(20) 

FNDIF 

WR  (  f  E  ( * ,  2 )  •  L  AUNCH  SHE  ELE VAT  ION  C ET 3  ’ 

KEAD(*,FI-IAT)X(21) 

WPirF(<t,2)’  GREENWICH  HOUR  ANBuE  rtlegrees3’ 

RtAl)(  »,FMAT  )X  (22) 

WRr.r;(>i,2)'  UECLINATION  tdec3rBes3’ 

Kr:AD<«,EMAr)X(2’.) 

Wl'-  n  I-  '  •  .  2 '  ■  1  ( )N(?  I TUITE  C  dsgr-  3  ’ 

RI.  MAI  )  X  (24) 

WRIT!;  ( 1,2)  •  lATlTUDE  Cd(?gr«es]' 

REAI)(*  .FMAT  )X  (2‘.-.) 

WRITC(t(,2>'  TEMFTRATURE  STABILIZATION  TIME  Cr>ec3’ 

READ ( 1 ,  FMAT ) X (26) 

WK!rE(«,2)'  NUMBER  OF  LOCAL  ATMOSPHERE  POINTS  C.'50  max.  3’ 
RI:  ADI  «,(  MAT  )  .<  (2T) 

WRt  ri:  (4,2)  ’  NUMBER  OF  FLIGHT  PROFILE  POINTS  Cl  00  maxJ  ’ 
RL:AI)<i,FHAr)X(2Et) 

IM^INTvX (27) ) 

INN-INT(X  (213)  ) 


DO  lOu  1=1,70 
Y<  1,  I)=^0. 
Y<I,2)--0. 

DO  200  I=r.l,lNN 
Z<I,  1)=:C/. 

Z ( (,2)=0. 
Z<1,B)=0. 
2(I,4)=:0. 


Li  VT'2^  *.i  *  '■.4  * 


2r..VS :  -V-  ’\'%.-'J’~\T'\ ' 


•/  'fTtf  ■ 


C  INPUT  LaCAL  ALTITUOa-tEMPERATURE  PROFILE 

WRITEU,!!)’  LOCAL  ATMOSPHERE:  ALT.  CSkml  -  TEMP.  CKelvinl' 
WRITE(*.2)’  IF  RANKINE  ENTER  l.SJ  IF  KELVIN  ENTER  l.O’ 
READ(*,FMAT)TCON 
FMATs’ (F9.3) ' 

DO  300  1  =  1,  IN 

WRITE(*,2>'  ALTITUDE  CQKm 3' 

READ(*.FMAT)Y(1, 1) 

WR I TE ( * ,  2  > '  ABSOLUTE  TEMPERATURE ’ 

READ  (*,FMAT) TABS 
300  Y(I,2>=TABS/TCaN 

WRITE(*,11)’  FLIGHT  PROFILE:  ALT.  tltl  !<  ELAPSED  TIME  Iser.3' 
350  WRITER,  2)’  ENTER  ALT.  SOURCE  NO.:  1  =  PRESSURE.  2  =  RADAR  ?’ 
READ t*. 3) ITEM 

IF  <  (ITEM.LT.D.QR.  <1TEN.GT.2>>  SOTO  350 

C  COMPUTES  STD.  ATM.  PRESSURE  FOR  INDICATED  FLIGHT  ALTITUDE 

FMAT=' (FB.O) • 

DO  <5.00  1  =  1,  INN 
WRITE(»,2>’  ALTITUDE  Elt3‘ 

READ(*.FNAT)I(I,2> 

WRITE(«,2) ’  TIME  CsecI’ 

REAl)(*,FMAT)Z(i,l) 

IF  IITEM.EQ.l)  THEN 
Q=.304a03V*Z (1,2) 

0=6354. 766«Q/ <6356766. *0) 

IF  ((0.6E.E(JD,n  ).AND.  (O.LV.ED)  GOTO  SOO 
400  JD=OIFF 

JIFF»JIFF+1 
E1=E<0IFF, 1) 

E2=E(JIFP,2) 

1D«E2-E;J!>,2> 
l.))(wLt'6.  (JD,  I ) 

R=TD/I)X 

IF  (O.I.r.Et)  SOTO  5(30 

1 F  ( APS  V  T  in .  oc:  .0.0'.'  H I  =H  I  *  ( \  F.  ( ,1  u ,  2 '  /  e :;  >  *  » ( cc  /  R )  > 

IF  ( AB  S  ( T  D ) .  L  T  .  0 . 0 ! )  H I  =  H I  * E  X  P  ( -CC  * !)  X  /  E  <  .1 1) ,  2 )  > 

GOTO  400 

500  DX--0-E(,ID,  I) 

TF=i; j)+R*nx 

I  F  (APH(  TO)  .(SE.o.ot )  Z  ( 1 »( 1E<  JD.2>  .'It  )  *«  (CL/R)  > 

I F  ( AB5  <  TI)> .  LT  .  0. on  Z  ( 1 . 3)  =H1  »EXP  ( -CC*DX/U  ( Jl),  2)  ) 

END  IF 

600  (4RI  TE  (»,  4)  Z  ( I.  1  )  .  Z  ( 1.2' ,  Z  ( 1 ,3> 
ir  (I  lEM.EO.;'.'  GOTO  210'' 

fi  1  *5  w  ii ! , 

JO  - IN 

C  Cnr.Ri-(  IS  ALTITUDE  FOR  LUCAI.  ATMOSPHERIC  PRESSURE 

DO  20>.'" 

AL-  Z  (I,/) 


.V 


ii 
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vV..'V.’V  V.V 


%  •- 


j 


DAL=BOO. 

DZ3a0.00S»2<l,3) 

700  a-.3048037»flL 

D-=63bi. .  766»tt7  <4356746.  +C> 

II"  (!«.G!K.Y<J|>,  l)).AND.  (O.UT.ED)  BUTO  900 
IF  (Q.GE.El)  GOTO  000 
aiFT=»l 
H1=X(14) 
aOO  00=^0  IFF 

olFI-=vHFF-tl 
EleYIJIFF, 1) 

ES=:Y(Jir-'F,2> 

TDsE2-Y<jb,2) 

DX=E1-Y<JD,1) 

9--=TD7DX 

!*•  (CJ.LT.El)  SOTO  900 

U  (AIJS(TD)  .GE.O.OU  H1=H1*  ( ( Y  (JD,  2) /E2)  *»  (CC/R>  ) 

IF  (ABS(TD).LT.O.Ol)  Hl=Hi«EXP (-CC«OX/Y ( JD, 2> ) 

GOTT  BOO 

900  l)X=tV-Y(JI).  1) 

TF«VvJD,2)+R»DX 

IF  <AOS<TD)  .GE.O.On  BP=HHK  <  Y ( JD.  D /TK)  «•* (CC/R)  ) 

IF  (AB3(TD> -LT.O.Ol)  BF=Hl*EXP(-CC»0X/Y(aD,2) > 

C  PRESSURE  CONVERGENCE  CHECK 

GI=AIQ«KBP-  2  (1,3)  > 

C  CORRECTS  MONOrONlC  /ASSUMPTION  IF,  REQUIRED 

IF  (IQ.LT.2)  THEN 
IQbIQ+I 

IF  (IQ. EG. 2)  THEN 

ir-  <(  (0*000)  .6T.O.  )  .AND.  (APS(PPQ)  .GT.ABSIQ)  )  )  THFN 
AUA=-l. 

AC 3-0. 

AC4oO. 

AL=Z (I .2) 

DAL-QOO. 

ROTO  7<.*0 
FNDIF 
GOTO  950 
FNDIF 
OUO-Q 
HMD  IF 

C  LIMIT  CHECF 

950  IF  (ADS(e).GT.DZ;.)  THEN 

C-  NOT  WITHIN  LIMITS 

IF  (D.IT.O.)  GOTO  tOOO 
AC 3“ I . 

S(5N^l. 

HO  10  1 ! OO 
lOOo  0C4=I. 

SGN=-I. 

I  to*.'  DAL^^DAI./  ( ACT-* At;/) ) 


I 

i 


\ 
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AL=AL-^S(iN*DAL 
GOTO  700 
END  IF 

C  WITHIN  LIMITS 

AC3-0. 
ftC4=0, 

Z  <I.4>sZ<t.2) 

Z  <1,2)=AL 

■SOOO*  WRITE  (♦,  S)  2  ( ! ,  2) ,  2  ( 1 , 4) 

2100  0FEN<,5,FlLe«FINN,STATUS=*NEW'  ,FORtt>=* FORMATTED* ) 
WR1TE{3,  lUFlNN 
FMAT*’ (EIS.O) * 

DO  2200  l“l,2e 

2200  WRnE(3,FNAT))(<l) 

FNAT«*  r^F'P.S)  ' 

DO  2300  I » 1,30 

2300  WRITE (3, FNAT) Y ( 1 , 1 ) , Y 1 1 . 2) 

•  fMAT«’ (2F8.0) ’ 

1)0  2400  1^=1,  INN 

2400  WRITE(3,FMAT)2(l,l),ZfK2> 

CLOSE (3) 

END  IF 

C  CREATl-S  BJ  TUor.pec.GAF 
WRITE (♦,?) 

READ<».3) 1AM 
IF  <!AM.ER.1)  THEN 
WRITE!*, 9) 

REAIX'MUFINN 
FNAT** (ElS.a) * 

WRITE (  U  2)*  ALF* 

REAl)(ii,FMAT)XU) 

WRITE (♦.2)’  BET* 

RE  AD  (It.  FNAT)  X  (2) 

WRITE!*, 2)’  GAN* 

READ!*,FMAT) X!3) 
wRrn;!»,2)’  c:v’ 

READ!*.  FNAT) X (4) 

WMn:(*,2)'  CF’ 

RKAD(«,FMAr)X(5) 

WRITE!*, 2)'  WOW' 

RFAO(*,FMAT>  XvT,) 

WRIIE(*.2>*  A/V 
RF.AI)«*,FMA1 )  X  <  :*' 

WR  Id.' <  k , ■  im'’ 

Rf.'AI><«.FMm)x<H) 

WRiril(*.2)'  IW 
READ!*,  FNAT 'X!9> 

WRUT!*.!*''  AXir 
RK.’iD**  .  •  X  t  !».' 

WRl  IE  ■  X  .  r>  ■  E'.F' 

RF:a1i'»  ,  FNAT'  X  .'111 


WF?ITt:(*,2)  RYR' 

F!t-;AI){*.FMflT>X(12) 

WRITfi(*,2)’  TYR' 

REftD(*,FMAT>  X (13> 

WRITE!  AYRG’ 

re;ad(*.fmat>x(14) 

aPEN(3,FIL.E=FINN.STATUS='NEW'  ,FOftM-’f-ORMATTED’  ) 
WRITE(3,  IDFINN 
DO  2600  1=1,14 
2600  WR I TE ( 3 , FMAT ) X ( I ) 

C:i-0SE(3> 

END  IF 

C  CREATES  B:f ilespec.CMV 

WRITE(»,a) 

READ!*, 3) I AM 
IF  (lAN.EQ.t)  THEN 
WRITE!* . 10) 

READ!*,  1  DMNN 
FMAT=- (EtS.9> ' 

WRITE!*, 2>'  CRl’ 

READ!*,FMAT)X!1) 

WRITE!*, 2)'  C02' 

READ!*, FMAT) X (2) 

WRITE!*, 2)’  CQ3' 

READ!*, FMAT) X  <3) 

WRITE (*,2)’  C04' 

READ!*. FMAT) X<4) 

WRITE!*,2)‘  CG5’ 

READ<*,FMAT>X(S) 

WRITE!*, 2)’  C06' 

READ!*, FMAT) X< 6) 

WRITE!*,  2)-  Ct!7' 

READ!*, FMAT) X!7) 

WRITE!*, 2)'  COB’ 

READ!*, FMAT) X !0) 

WRrrE!i,2)’  CQ9’ 
re;ad!*,fmat)X(V) 

WRITE <*,2>'  COIO’ 

READ<*.FMAT)X(10) 

WRITE!*, 2)’  coir 
READ!»,rMAT)X(ll) 

WRITE!*,  2)’  ewr 

READ !*, FMAT) X( 12) 

WRITE!*.  2)'  GNl” 

F!EAD!*.F'MAT)X!13) 

WRITE!*. 2)’  GN3* 

READ ! * . FMAT ) X ! 1 4 ) 

WRITE!*, 2)'  VIRTUAL  MASS  COEFFICIENT’ 

READ !*, FMAT) X 115) 

WRITE!*, 2)’  EFFECTIVELY-ZERO  VELOCITY' 
READ1«,FMAT)X116) 

0PEN!5,FILE=FINN,STATUS='NEW’  , FORM=’ FORMATTED’ ) 
WRITE13,  IDFINN 
DO  2000  1=1,16 
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t'lM  (T  I  ’.,1  HAT)  X '.1> 

ri  usT.  ( ■'.) 

IMOU 

tMc:^ 


tUFiAOm  INE  INFORM 
:■  SEE  ISEI^ 

DIMENSION  LINE<25> 

i;haracier«vi'  lime 

CHARACTER*  1 S  t-T  NN ,  FNAME 
I UFVMAT  (13) 

IOr:MAI(A) 

FORMAT  (.4X,A) 

4  FORMAT (A, \) 

FlNN=’t!Cl)r)ATA.FAX’ 

L.  -11 

OHEN<  I  ,F'  Il.E=FINN,EORM=’EORMArTEL''  > 

ri:AD(T,;')FMAME 

FlADCl,  1 )  1 

IF  (I.,NE.67B)  (30  TO  9 

C  1:  NIJMB(-:R  of  lines  of  text  unit  to  be  DISFLAYED 

r,  READ  (1,1)1 

U  (l.EO.O)  BCJ  TO  7 
I'l-C.i 

IF  <I.LE.20)  Nr--(24-I>/2 
r,'Al  1.  a.EAR(2S) 

Dll  7  F.  =  1,I 

/  READ  (1.2)  LINE  O  ') 

1)0  a  1-1,1 

0  WRI  TE  («,  3)  L  INE  ((■■) 

IF  •'N.OI  .  0)  CALL  CLEARTN) 

PAUSE 

IF  (I...ED.  t  )  OOTO  h 

WryITE(*,4)  niBPLAY  PROGRAM  NOTES  ?  t  0/1  =  N/Y  3 
READU',  1)L 
IF  (L.EO. 1 )  OOTU  6 
7  CI..OSE  <  1 ) 

■:AL(.  clear  (25) 

END 


SUBROUTINE  CLEARIJ) 
1  FOFI-MATIA) 

DO  2  1:^1,. T 
7.  WRITE(*,1)’  ' 

END 
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GLOSSARY  OF  PROGRAMS  'FINDCD'  &  'FRDUDE 


NAME  REFERENCE  DIMENSION  DESCRIPTION 

A<N)  constant  -  RUNGE-KUTTA  AND  OTHER  CONSTANTS 

ABL  UPSON  -  BALLOON  SURFACE  AREA  i  6L*#2 

AC3  default  -  CONVERGENCE  ROUTINE  PARAMETER  <  0  or  1  ) 

AC4  default  -  CONVERGENCE  ROUTINE  PARAMETER  (  0  or  1  > 

AD  input  ft**2  DUCT  AREA,  TOTAL 

AG  Model  -  EFFECTIVE  UV  ABSQRPTANCE  OF  GAS 

AGAS  input  -  IR  RESPONSE  CONTROL,  GAS  C  0  or  1  3 

ALF  def/inp  Ibrn/ft/ser  VISCOSITY  COEFFICIENT,  GAS 

ALIM  input  ft  ALTITUDE  GROWTH  LIMIT  IN  TIME  DT 

AM  Model  mass  units  OPTICAL  AIR  MASS 

AR  -  EFFECTIVE  IR  ABSDRPTANCE,  FILM  (  =  ER  > 

ARMS  Model  mass  units  OPTICAL  AIR  MASS 

ASTOP  input  ft  ALTITUDE,  COMPUTATION  ABORT  (see  notes) 

AV  Model  -  EFFECTIVE  UV  ABSDRPTANCE  OF  FILM 

AYR  def/inp  -  COEFF.  OF  ABSORPTIVITY,  IR 

AYRG  def/inp  -  COEFF.  OF  ABSORPTIVITY,  SOLAR  (gas) 

AVR61  formula  -  <1-AYRB) 

AYV  def/inp  -  COEFF.  OF  ABSORPTIVITY,  SOLAR 

B(N)  constant  -  RUNGE-KUTTA  AND  OTHER  CONSTANTS 

BE  formula  lb/ft**3  SPECIFIC  LIFT 

BET  def/inp  -  VISCOSITY  EXPONENT,  GAS 

BP  '/IRON  mb  PRESSURE,  ATMOSPHEKTC  (not  u^jed) 

BLKnV  lb  BUOYANCY,  TOIAL 

BZ  constant  Ib/f t/sGr:/R*#4  STEFAN- BOLT Z MAN  CONST.  (3.6995E-10> 

C(N).  constant  -  RUNGE-KUTTA  AND  OTHER  CONSTANTS 

CA  VIRON  lli/sec/R  COEFFICIENT.  THERMAL  CONDUCTIVITY  OF  AIR 

r:c  constant  GK«i/Kelvin  VALUE  (3A.  lf,3195) 

CD  -  COEFFICIENT,  DRAG 

(::F  def/inp  ft/Rankine  SPECIFIC  HEAI  OF  BALL(30N  FILM 

C6  formula  lb/«ec/R  COEFFICIENT,  THERMAL  CONDUCTIVITY  OF  GAS 

CLEN  input  ft  LENGTH,  CAP 

CM  def/inp  -  COEFFICIENT,  VIRTUAL  MASS 

C(i!l  clef/inp  — . -  CORRECTION  FACTOR  FOR  01 

CO  10  def/inp  -  CORRECTION  FACTOR  FOR  DIO 

tXUl  def/inp  -  CORRECTION  FACTOR  FDR  D1  1 

CD?  def/inp  -  CORRECTION  FACTOR  FOR  02 

CG!3  def/inp  -  CORRECTION  FACTOR  FOR  03 

C04  def/inp  -  CORRECTION  FACTOR  FOR  04 

COS  def/inp  -  CORRECTION  FACTOR  FOR  05 

CG'o  def/inp  -  CORRECTION  FACTOR  FOR  06 

CD'/  def/inp  -  CORRECTION  FACTOR  FOR  C!7 

CDS  def/inp  -  CORRECTION  FACTOR  FOR  OG 

CD9  def/inp  -  CORRECTION  FACTOR  F(DR  09 

CS  default  -  SOLAR  RADIATION  FACTOR 

CV  def/inp  ft/Rankine  SP.  HEAT  AT  CONST.  VOL.,  GAS 

C)  formula  -  COS  ( XL)  *CO.S  ( Xl>) 

DiM.M)  formulas  -  TINE  DERIVATIVE  OF  Y(M,N) 

PAl  ft  ALLOWABLE  ALTITUDE  CLOSURE 


DB 

default 

lb /sec 

DBB 

Ib/sec 

DBB 

input 

Ib/min 

DCD 

input 

DDO 

VI  RON 

GKffl 

DDl 

VI  RON 

GK(t» 

DD2 

VI  RON 

Kelvin 

DD3 

VI  RON 

Kel' in 

DD4 

VIRON 

mb 

DD5 

VI  RON 

mb 

DD6 

VIRON 

Kel vin/GKm 

DD7 

VIRON 

Kel  vin/GKd) 

DDB 

VIRON 

Kelvin 

DD9 

VIRON 

Kel vi n 

DM 

MYBLN 

ft 

CRAG 

formula 

lb 

DT 

dcf aul t 

sec 

DTI 

Model 

Rankine/GKm 

DTM 

def /i np 

sec; 

DTV 

clef  Zinp 

sec 

DTX 

seu 

DUl 

input 

DU2 

i  nput 

— 

DU3 

i  nput 

CU-1 

i  nput 

DUS 

input 

E(N, 1) 

input 

GKm 

E(N,2) 

1  nput 

Kel vin 

EG 

Model 

EGAS 

— 

El 

Model 

- - 

ELL 

1  nput 

ft 

ER 

Model 

ERR 

input 

EYR 

def /i np 

EYRG 

Model 

EYR(31 

f  cirmul  a 

- - - 

F(N) 

f ormul as 

(van  able) 

FE 

Model 

Ft 

i  nput 

FLYiN, 1 ) 

input 

sec 

FLY (N.2) 

j  nput 

ft 

FORCE 

lb 

FR 

— 

KRD 

— 

FROM,  1) 

out  put 

— 

FR(N.2) 

LIU  (.put 

• -  -  ~ 

FR(N.:.) 

C'UipiM 

— 

FR  (N,'T) 

OLltput 

FV 

Model 

Ib/it/sec 

0 

constant 

ft  .'S.e'C»«2 

GAM 

dfc'f  / 1  np 

ft  Ih/lbm/R 

GD 

const  ant 

— 

6B 

MYBLN 

— 

I3B1 

input 

( var  i  dbl  n) 

ACTUAL  DEBALLASTINB  RATE  (default  is  0) 
MAXIMUM  DEDALLASTING  RATE  (converted) 
MAXIMUM  DEBALLASTING  RATE 
CD  ADJUSTMENT 

STORAGE  C  EF  ]  for  altitude  ZZ 

STORAGE  C  EF  3 

STORAGE  C  EG  3  for  altitude  ZZ 

STORAGE  [  EG  3 

STORAGE  C  RX  3  for  altitude  ZZ 

STORAGE  C  PX  3 

STORAGE  C  RR  3  for  altit>  dc  ZZ 

STORAGE  C  RR  3 

STORAGE  C  DY  3  for  altitude  ZZ 

STORAGE  C  DY  3 

DIAMETER  OF  BALLOON 

DRAG,  AERODYNAMIC 

TIME,  INTEGRATION  INCREMENT 

TEMPERATURE  GRADIENT,  IR 

TIME,  MAX.  INTEGRATION  INCREMENT  <  20.) 

TIME.  VENTING  INTEGRATION  INCREMENT  <  O. 

iNTERt^EDIATE  VALUE  OF  DT 

DUMMY  (unused) 

DUMMY  (unused) 

DUMMY  (unused) 

DUMMY  (unuBCJd) 

DUMMY  (unused) 

ALTITUDE 

TEMPERATURE,  ATMOSPHERIC  PROFILE 
EFFECTIVE  IR  EM I SSI VI TY  OF  BAS 
IR  RESPONSE  CONTROL,  GAS  t  =  ABAS  3 
RADIATIVE  EXCHANGE  COEFF.  ,  GAB  S<  WALL 
fc  l.t VA  T 1  ON ,  LAUNCH  S 1  I E 
EFFECTIVE  IR  EMISSIVITY  OF  WALL 
ALLOWABLE  R.E.  IN  COMPUTED  ALTllUDE 
COEFF.  OF  EMISSIVITY,  IR 
COEFF.  OF  EMISSIVITY  (yaa) 

(  1  ••  EYRG  ) 

(locally  defined) 

DIRECTIONAL  HEMISPHERICAL  REFLECTIVITY 

FREE  LIFT.  FRACTIONAL 

TIME 

ALTITUDE 

LIFT,  NET 

FROUPE  NUMBER 

FROUDE  NUMBER 

REYNOLDS  NUMBER 

FROUDE  NUMBER 

FRACTIONAL  VOLUME 

DRAG  COEFFICIENT 

UV  FLUX 

DRAVITATIONAL  CONSTANT  (32.17.11) 

THFRMAI.  rONUUCTTVITV  COFFF.  ,  GAS 
(  CLEN/GL  ) :  used  once  as  DUMMY 
(  GN/6L  ) 

DLIMMY  (  locally  defined  ) 
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>  V  \ 


i 


m 

cr 

* 

h 


BB3 

input 

<var i able) 

DUMMY  (  locally  defined  ) 

GH 

input 

decireejs 

GREENWICH  HOUR  ANGLE 

Gl. 

i  nput 

ft 

(50RELENGTH,  GALLOON 

GLNIO 

constcint 

— 

NATURAL  LOG  OF  10 

GN 

NYCd-N 

ft 

GORELENBTH,  GENERATOR  SHAPE 

GNl 

def /inp 

CORRECTION  FACTOR 

BN2 

de-f  /  i  np 

CORRECTION  FACTOR 

GN3 

def / i np 

CORRECTION  FACTOR 

GNU 

Model 

NUSSELT  NUMBER 

BP 

PRBR 

6RASH0F  NO.  *  PRANDTL  NO. 

constant 

Ib/f t/sec 

SOLAR  CONSTANT  (96) 

HC 

MYBL.N 

f  t'lm^ 

AREA.UV  ABSORPTION 

IN 

— 

IKTIXIN) 

INN 

— 

INTtXlNNI,  Kubiect  to  option 

INN 

input 

— 

NO.  OF  FLIGHT  DATA  POINTS  TO  BE  ANALYZED 

IVPNI 

— 

STATUS  OF  DUCT  VENTING  <  Y/N  1/0  ! 

JSEG 

input 

— 

NUMBER  OF  DRAG  MODEL  SEGMENTS 

KI 

— 

INDEX,  RUNBE-KIJTTA 

!■.  .*f 

INDEX,  RUNGE-l'UTTA 

[■X 

— 

INDEX,  RLINGE-KUTTA 

LAUNCH 

de-f  aul  t 

— 

STATUS  OF  LAUNCH  (  Y/N  =  1/0  ) 

LEAP 

default 

— 

STATUS  <  TT.GT.DT  /  TT.LE.DT  =0/1  1 

LL.0 

VI  RON 

— 

STORAGE  C  JX  I  for  altitude  12 

LH 

VI  RON 

STORAGE  t  JX  J 

LL2 

VI  RUN 

— 

STORAGE  1  JV  ;l  for  altitude  ZZ 

ll::. 

VI  RON 

— 

STORAGE  I  JV  I 

PO 

1  nput 

mb 

PRESSURE,  ATMOSPHERIC  AT  LAUNCH  BITE 

PAM 

VIRON 

mh 

PRESSURE.  ATMOSPHERIC  (at  alt.  Z2) 

f’E 

VIRQN 

lh/ft*iF2 

PRESSURE.  ATMOSPHERIC 

PCs 

formula 

111-  f  t**2 

PRESSURE  DIFFERENTIAL,  BALLOON  APEX 

PMM,N) 

f  orniulaT. 

'■  var  i  abl  e'j 

VARIABLES.  RUNGE-KUTTA 

l.'l 

Model 

ft  lb 'Set. 

CONVECTIVE  HEAT  TRANSFER,  GAS  «!  WALL 

010 

Model 

ft  lb -•  set 

INFRARED  EMISSION.  GAS 

Q1 1 

Model 

ft  lb /sec. 

INFRARED  ABSORPTION,  GAG 

02 

Model 

ft  lb/ sec 

DIRECT  SOLAR  ENERGY  ABSORPTION 

03 

Model 

ft  lb /sec 

I NFRARED  ADSORPT 1 ON 

04 

Model 

ft  lb /see 

CONVECTIVE  HEAT  TRANSFER,  WALL  !•  AIR 

OS 

Model 

ft  lb /sec 

INFRARED  EMISSION 

06 

Model 

ft  Ib/sec 

REFLECTED  SOLAR  ENERGY  ABSORPTION 

cr; 

Model 

ft  lb ’sec 

RADIATIVE  EXCHANGE,  GAS  !<  WALL 

08 

Model 

ft  lb/ sec 

DIRECI  SOLAR  ENERGY  ABSORPTION.,  GAS 

ov 

Model 

ft  1  b  /  ?iec 

REFLECTED  SOLAR  ENERGY  ABSORPTION,  GAS 

QA 

f ormul a 

— 

CUS (SOLAR  ZENITH  ANGLE) 

OG 

ormula 

COS  ‘.f  cntal  t  itudel ) 

00 

formal  a 

<  var  i  abl  e> 

DUMMY,  RUNGF-kllTTA 

FVA 

iionetant 

ft /Rani i ne 

GAS  CONSTANT,  AIR  (SS.ZSS) 

PE 

c  cinst  ant 

ft 

RADIUS  OF  EARTH  (20,0^.5,278) 

PPV-M, 

N)  input 

— 

EXPONENTS  AND  COEFFICIENTS.  DRAG  MODEL 

const  ant 

ft /Rant  1 ne 

GAS  CONSTANT,  HELIUM  (3B6.076) 

PI. 

Model 

— 

REFLECTANCE 

MYULN 

It 

RADIUS,  GAS  BUBBLE 

l;N 

formula 

REYNOLDS  NUMBER 

RD 

VIRON 

slug/f t**3 

density,  ftIR 

PYR 

clef  / 1  'ip 

— 

COEFF.  OF  REFLECTIVITY,  IR 

PVV 

def  •■'i  np 

- . - 

COEFF.  OF  REFLECTIVITY.  SOLAR 
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SA 

Model 

SB 

I'tYBLN 

SIG 

input 

SILTX 

UPSON 

- - 

SPD 

1  nput 

Knots 

SPDSO 

(-ft/sec>«-*2 

SPEED 

Tt/sec 

SR 

Model 

f  t#*2 

SW 

VI  RON 

lb/tt**3 

SWG 

far  mill  a 

U)/fttnf3 

SX 

T  ormul a 

T 

!5ec 

TCf^ 

i  nput 

111  i  1  u 

THI 

1  iiput 

mils. 

TI 

VI  RON 

Rankine 

TIPI 

Model 

Rank  me 

TIRO 

Madel 

Rank  j  lie 

Tt.' 

VIRON 

Kel vin 

TLIM 

i  nput 

Rank;  me 

TR 

VIRON 

Rankine 

TRMl 

Model 

trm:? 

Model 

TIVOP 

Model 

6K(ifi 

T5 

input 

T'lT  OP 

1  nput 

sec: 

1 1 

nt5c 

TT 

ITT 

TX 

VIRON 

(-.el  vm 

TYR 

lieF  / 1  np 

■ - 

TYV 

dot / j  np 

U 

NELSON 

Uu 

NELSON 

Ft*-*3 

vn 

— 

VI) 

■formula 

ft**:). 'set: 

VL 

def / i np 

•ft. 'sec 

VLIN 

i  nput 

f  t /sec 

VS 

VIRON 

lb. 'ft. 'sec: 

VBiS 

f  ar  mu  1  a 

1  b.'f  t  ,'sec 

VT 

f  l*»3 

VTM 

1  nput 

f  t**3 

VU 

UPSON 

— 

VV 

Model 

— 

VV 

formula 

f  t*if3/set: 

W(M,N) 

NELSON 

— 

WB 

i  nput 

lb 

wc 

f  cirmul  a 

lb 

WE 

■f  ormui  a 

lb 

WOW 

de-f /inp 

Ib.'f  tif*2/mil 

WP 

i  nput 

lb 

ws 

lb 

WT 

1  nput 

Ih 

WTX 

lb 

WV 

5et:/f  t**2 

EFFECTIVE  SURFACE  AF^EA  OF  C3AS  BUBBLE 
AREA  OF  BALLOON  BUBBLE,  FRACTIONAL 
SIGMA,  BALLOON  SHAPE  FACTOR 
SIGMA,  IN  FLIGHT 

LAUNCH  WIND  SPEED  (c:onverL«ci  to  Tt/aec) 
(  not  used  ) 

ASCENT  RATE  (  or  launch  wind  speed  > 
EFFECTIVE  IR  SURFACE  AREA 
SPECIFIC  WEIGHT,  AIR 
SPECIFIC  WEIGHT,  GAS 
'!IN(XL>  *  SINIXB) 

TIME,  TOTAL  ELAPSED 
THItXNESS  OF  CAPS  <suni> 

THICKNESS,  BALLOON  WALL 
TEMPERATURE,  EOUILIBRIUM  RADIATION 
TEMPERATURE,  IR  <  at  trapopause  ) 
TEMPERATURE,  IR  !  at  insl  ) 

TEMPERATURE,  ATMOSPHERIC 
TEMPERATURE  GROW I II  LIMIT 
TEMPERATURE ,  AT MORPHER I C 
TRANSMITTANCF,  ATMOSPHERIC  (solar) 
TRANSMITTANCE,  ATMOSPHERIC  (solar) 
HEIGHT,  TROPOPAUSE 
TINE.  TEMPERATURE  STABILISATION 
TIME.  CDMPIJTATION  ABORT  (goc)  notes) 
TIME,  ELAPSED  PRIOR  TO  STABILI 2A1 ION 
TINE,  REMAINING  IN  INTERVAL 
S  TOR-AGE  C  T  1 

TEMPERATURE  AT  ALTITUDE  ll  (not  used) 
COEFF.  OF  TRANSMISSIVITY,  IR 
COEFF.  OF  TRANSMISSIVITY,  SOLAR 
FRACTIONAL  PART  01-  00 

ratio:  cap  area  /  balloon  surface  area 

VOLUME  OF  GAS,  PRELAUNCH  (no  superheat) 
VOLUME  OF  BAl.LOON,  FRACT  lONAt. 

VOLUME  FLOW  RATE  THROUGH  DUCTS 

EFFECTIVE  2ER0  VELOCITY 

VELOCITY  GROWTH  LIMIT 

VISCOSITY,  AIR 

VISCOSITY.  GAS 

VOLUME.  AT  CEILING  ALTITUDE 

VOLUME  OF  BALLOON.  THEORETICAL  MAXIMUM 

RATIO:  INSTANTANEOUS  VOLUME  /  6L*»- 

DISCHAR6E  COEFFICIENT.  EV-13 

VOLUME  FLOW  RATE  THROUGH  EV-13 

COEFFICIENTS.  CAP  AREA 

WEIGHT.  BALLOON 

EFFECTIVE  WFIIGHT.  CAPS 

EFFECTIVE  WEIGHT.  BUBBLE  ENVELOPE 

UNIT  WEIGHT  (  POLYETHYLENE  =  0.0046  ) 

WEIGHT.  IRREDUCIBLE  PAYLOAD 

WEIGHT,  SYSTEM 

WEIGHT,  BALLAST 

STORAGE  r  WT  3 

STORAGE  r  SW/VS  3 


X 

1  .-ar  1  ahl  e?) 

niJMftX  <  1  oc  a  1 1 V  dt'f  1  titsd ) 

XI) 

input. 

dour GPS 

URCtlNAI  KJN,  LAUNCH 

XB 

1  nput 

degr  ees 

LONGITLIOe,  lAUNOI 
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